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Abstract 

The field of astroparticle physics is currently developing rapidly, since new exper- 
iments challenge our understanding of the investigated processes. Three messen- 
gers can be used to extract information on the properties of astrophysical sources: 
photons, charged Cosmic Rays and neutrinos. This review focuses on high-energy 
neutrinos {E^, > 100 GeV) with the main topics as follows. 

• The production mechanism of high-energy neutrinos in astrophysical shocks. The 
connection between the observed photon spectra and charged Cosmic Rays is 
described and the source properties as they are known from photon observations 
and from charged Cosmic Rays are presented. 

• High-energy neutrino detection. Current detection methods are described and 
the status of the next generation neutrino telescopes are reviewed. In particular, 
water and ice Cherenkov detectors as well as radio measurements in ice and with 
balloon experiments are presented. In addition, future perspectives for optical, 
radio and acoustic detection of neutrinos are reviewed. 

• Sources of neutrino emission. The main source classes are reviewed, i.e. galactic 
sources. Active Galactic Nuclei, starburst galaxies and Gamma Ray Bursts. The 
interaction of high energy protons with the cosmic microwave background implies 
the production of neutrinos, referred to as GZK neutrinos. 

• Implications of neutrino flux limits. Recent limits given by the AMANDA exper- 
iment and their implications regarding the physics of the sources are presented. 
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1 Introduction 



The multimessenger connection between Cosmic Rays, photons and neutrinos 
of different wavelengths is crucial to comprehend in the pursuit of a deeper 
understanding of the fundamental processes driving non-thermal astrophysi- 
cal sources. None of these messengers alone is able to give a complete picture. 
While photons reveal the surface of those objects due to the optical thick- 
ness of the sources, charged Cosmic Rays give a direct sight into the inner 
acceleration processes. However, they do not point back to their origin, since 
scrambled by intergalactic magnetic fields. Neutrinos are produced in inter- 
actions of protons with a photon field or other protons. The detection of an 
extraterrestrial neutrino signal at the highest energies [E^, > 100 GeV) renders 
possible the exploration of the particle acceleration region itself. Neutrinos es- 
cape the acceleration region and propagate through space basically untouched 
and this advantage is directly connected to the drawback of the very low de- 
tection probability of neutrinos. In this review, the connection of the three 
messengers will be reviewed with the focus on the resulting neutrino fluxes 
and their detection probabilities. 

Astrophysical neutrinos are produced over the broad energy range, covering 15 
orders of magnitude. An overview of the neutrino spectrum is shown in Fig. 1. 
From the lowest energies of E,y ~ meV to the highest energies of E^y ~ EeV, 
the intensity of the signal decreases by 42 orders of magnitude, making it 
necessary to introduce new methods of neutrino detection and analysis in 
order to increase the sensitivity to the neutrino fluxes especially at the highest 
energies. 

The Cosmic Neutrino Background (Cz/B) is an isotropic neutrino flux having 



decoupled in the early Universe, only 1 s after the Big Bang, see |Wei72 
and references therein. It is the neutrino equivalent of the cosmic microwave 
background (CMB). The temperature of the blackbody spectrum has dropped 
to ~ 1.9 K today due to the expansion of the Universe, and the flux peaks 
at milli-eV energies. While this background is essentially predicted in the 
standard model of cosmology, it was not possible to test it experimentally yet. 
Recently, the measurement of the Ci^B was proposed by u sing the i nteraction 



of cosmic neutrinos with nuclei undergoing beta decay [CMMOTl ]. with 
as an optimal candidate. For the case of a 100 g ^H-detector, an event rate 
of a few tens to a few hundreds is expected. The exact number depends on 
the distribution of cosmic neutrinos and the neutrino mass. This is currently 
not within experimental reach, but can probably explored in the future with 
improved experimental techniques. 

The sun emits neutrinos in different fusion processes in the MeV range. In the 
figure, neutrinos from pp interactions and the ''B spectrum is shown. Massive 
stars (M ~ 10^ ° • M^ti) are expec ted to emit neutrinos at ~ 1 MeV due to 



Silicon burning |OMK04l . IMOK06 
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At slightly higher energies lies the neutrino spectrum from SN 1987A. The 
expected flux from a SN at a distance of 10 kpc is indicated. Such a close 
supernova will be observed by todays neutrino telescopes. The total diffuse 
flux from SNRs ("relic"), however, is about four orders of magnitude lower 
and could not be tested yet. 

At energies of -Ej^ > 0.1 GeV, the measured spectrum of a tmospheric ne utri- 



nos is indicated, squares are measurements from Frejus B"'"87l. |D"^95| an d 



dots are AMANDA measurements from the years 2000-2003 [M+OTI . iMiinO? 



At the highest energies, a generic spectrum from GRBs is indica ted [WB97 



WB99|, as well as the maximum contribution from AGN [MPROll ] and the ex - 
pected neutrino flux from the absorption of protons by the GZK effect YT93 . 
These sources have not been observed yet due to the high atmospheric back- 
ground and hmited sensitivities of the detectors. During the past ten years, 
high-energy neutrino detection still succeeded in improving the instruments in 
a way that made it possible to cover five decades of energy and to improve the 
sensitivity by 17 (seventeen!) orders of magnitude. We have today reached the 
point where the most optimistic models could already be used to constrain 
the physics of different sources classes and where the more realistic models 
are being challenged by the upcoming generation of neutrino detectors. This 
review discusses their options and possibilities. 

In particular, the different source types being able to accelerate neutrinos to 
the highest energies are examined in more detail. The hypothesis of neutrino 
emission from these objects is reviewed quantitatively in the context of multi- 
messenger physics. Section 3 focuses on what is known about the non-thermal 
Universe from photons and charged Cosmic Rays. The connection from these 
two messengers to the third, complementary particle, the neutrino, is drawn. 
In Section 4, the focus lies on the methods used to detect neutrinos at the 
highest energies. What neutrino signal to expect from extragalactic sources is 
considered in the following sections. Possible neutrino emission from galactic 
sources is presented in Section 5. Active Galactic Nuclei (AGN) are reviewed 
in Section 6, the neutrino flux models from Gamma Ray Bursts (GRBs) are 
discussed in Section 7 and starburst galaxies are presented in Section 8. Fi- 
nally, production scenarios of neutrinos generated by the interaction of high- 
energy protons with the CMB are presented in Section 9 before closing with 
a summary in Section 10. 
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Fig. 1. The astrophysical neutrino spectrum including different source predictions 
ranging from meV up to EeV energies. Point source fluxes have been scaled b y 
1/(4 vr) in order to be comparable to diffuse spectra. Figure after Kos92l . RouOC | . 
The individual spectra are explained and referenced i n the text. The atmospheric 
prediction, averaged over the solid angle, i s taken from the atmospheric data 

are from the Frejus ex perimen t D"'"95ll (red squares) and from the AMANDA 
experiment (blue circles) [m^07 . MiinO?! ] . The fluxes based on mere predictions are 
shown as dashed lines. The solid lines represent those fluxes already measured. 
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2 Introductory notes 



2.1 On the representation of flux models in figures 

Figure 1 shows the neutrino flux models in a double-logarithmic representa- 
tion, implying that the dependence of the differential flux dN/dEy with the 
energy is shown as \og\dNy / dE ij]{\og By). In general, non-thermal particle 
spectra can usually be approximated by powerlaws, dN/dE = A ■ E'"". In a 
double-logarithmic representation, this leads to a straight line, 

log [dN/dE] (log E) = log [Ay ■ E""] = log [A] -a- log [E] . 

The slope of the line is given by the spectral index a and the y-axis intercept 
represents the normalization A. Such powerlaw spectra are observed in the 
case of photons and charged Cosmic Rays, and are expected in the case of 
neutrinos. In the following, all spectra will be shown in the double-logarithmic 
representation. Since many spectra are very steep [dN/dE ~ E'"^ — E~^) for 
all three messengers, it is useful to weight the y-axis with a power of the 
energy, E^, straight line, 

log [E"" ■ dN/dE] (log E) = log[Ay E-^^""] ( 1 ) 

= log [A] - (a-w)- log [E] . 

For w = 2, which is often used in the case of neutrino spectra, an E~^- 
spectrum is then represented by a flat line, which is convenient, since neutrino 
spectra at E^, > 100 GeV are often expected to be close to Ey~^'^ — Ey~'^'^. 
The double-logarithmic representation with a weight is chosen to simplify to 
read off the spectral index of the models and measured fluxes. 

2.2 Flux and limit notations 

Throughout this review, different fluxes and flux limits will appear. Usually, 
the flux at Earth is given per energy interval, dN/dE, in units of [dN/dE] = 
GeV~^ cm~^ s~^ sr~^. Point source fluxes have the same notation, but the unit 
[dN/dE] = GeV~^ cm~^ s~^. The convention for the different particle and 
electromagnetic fluxes is as follows: 

• Charged Cosmic Rays 

For the total spectrum of charged particles, the spectrum is written as 

^ = AcR- Ec%-- . (2) 

If only protons or electrons are considered, the following notation is used: 
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Photons 

Photon spectra are usually given as: 

^^^.E-. (4) 

The total power Ptotiyj), depending on the angular frequency u;, is written 
as 



Ptot oc (5) 



with s = —aj — 1. 



For the special case of TeV photons, the notation 



dNxeV f ETeV\ 



dEr.v '^'^ VlTeVy 



(6) 



is used. The spectrum is normalized at 1 TeV. 

Neutrinos 

The neutrino spectrum is usually multiplied by the energy squared, E^"^, 

E.' ^ = A. ■ E.-- . (7) 



Alternatively, the spectrum itself is given as 
with a,y = a' — 2. 



V 



V 

Neutrino flux limits are usually given in the form of $ = EJ' ■ dNy/dE^ and 
are denoted as follows: 

• Diffuse Limit (DL) given in units of GeV cm ^ s ^ sr ^ . 

• <|)^-'^: Stacking Limit (SL) in units of GeV cm^^ s~^, obtained for the point 
source flux from a certain class of AGN. The stacking method is explained 
in Section 4. 

• Stacking Diffuse Limit (SDL), derived from the stacking limit in 
the same units as the diffuse limit, GeV cm~^ s^^ sr^^. It is determined by 
taking into account the contribution from weaker sources as well as yet 
unidentified sources, present in a diffuse background. 

A similar convention is used to denote the corresponding sensitivities: 

• Diffuse Sensitivity (DS) in units of GeV cm~^s~^sr~^. 

• <|>P^: the sensitivity to a single point source in units of GeV cm~^ s^^. 
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3 The multimessenger connection 



Neutrino flux predictions are built on the direct connection between the ob- 
served Cosmic Ray spectrum and the non-thermal emission from astronomical 
sources. Here, it is reviewed how the different parts of the Cosmic Ray spec- 
trum can be connected to the various objects and source classes and how this 
in turn leads to the production of neutrinos. 



3.1 Cosmic Rays 



Charged Cosmic Rays (CRs) have been observed from energies of Ecr ~ 
10^ eV up to Ecr ~ 10^° eV. The spectrum is altered by the Solar wind for 
energies below about Ecr ~ 20 ■ Z GeV, where Z is the charge of the nucleus. 
At higher energies, the energy spectrum follows a powerlaw with two, possibly 
three breaks. The powerlaw structure can in general be explained by shock 
acceleration in astrophysical sources. In this section, the observed spectrum 
with its features is discussed as well as the acceleration mechanism built on 
stochastic acceleration of test particles at magnetic field inhomogeneities in 
astrophysical shocks. 



3.1.1 Observation of charged Cosmic Rays 

Already in the early 20th century, it was discovered that the Earth is e xposed 



to a co ntinuous flux of charged particles from outer space, see e.g. |Hesl2 
Kohial . Viktor Hess and others performed balloon flights proving that the 
ionization of the atmosphere increases with height. This contradicted the hy- 
pothesis that the flux of ionizing particles arises from radioactive matter in 
the Earth's rocks exclusively. 

Today, the spectrum of charged CRs, dNcR/dEcR, has been examined over 
a wide range of energies Ecr, using balloon experiments and satellites for 
low energies {Ecr < 10^^ eV) and Earth-bound experiments for high energies 
{Ecr > 10^"^ eV). The all particle energy spectrum of CRs is shown in Fig. 2. 
The spectrum is weighted by E^j^ in order to have a flatter representation of 
the very steep spectrum. The powerlaw behavior of the spectrum is clearly 
visible, 

^ oc EaT"" ■ (9) 
dEcR '''' ^ ' 

Two kinks can also be seen, referred to as knee and ankle. The spectral indices 
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for the different parts of the spectrum are WBM98I . |V"'"99 



O^CR 



2.67 for \og{EcR/eV) < 15.4 

3.10 for 15.4 < log{EcR/eV) < 18.5 
2.75 for 18.5 < log{EcR/eV) . 



(10) 



A second knee around Ecr ~ 4 ■ 10^^ eV i s discussed toda y, with an even 
steeper behavior up to the ankle, see e.g. Hor03l . lHKT07l | and references 
therein. The general spectral powerlaw-like behavior can be explained by 
stochastic particle acceleration in collision-less plasmas. 



3.1.1.1 Cosmic Rays & directional information Charged Cosmic Rays 
below Ecr < Z ■ 10^^ eV do not point back to their origin, since they are 
scrambled by interstellar magnetic fields. The strong connection between non- 
thermal emission from astrophysical sources and particle acceleration can, 
however, be used to establish a model for different sources and source types 
to explain the Cosmic Ray spectrum. Sources within our Galaxy can produce 
the Cosmic Ray spectrum up to the ankle. Events at higher energies have to 
be extragalactic: 

(a) no galactic source class is energetic enough for the production of particles 
at such high energies as discussed later in Section 3.1.2 and 

(b) the particles' gyro-radius becomes too large and they escape from the galaxy 
already at lower energies. 

(c) In addition, at energies as high as Ecr ~ 10^° eV, the parti cle dif fusion 



is low compared to the traveling length through the Galaxy Hil84l ]. The 
observed particles point in this case back to their original source. The ob- 
served events are isotropically distributed, which is only possible for travel- 
ing lengths longer than the diameter of the galaxy. 

This leaves extragalactic or exotic sources as the origin of the highest energy 
events, for which acceleration up to Ecr ~ 10^^ eV is possible. At the highest 
observed energies, E > 5- 10^^ eV, nearby galaxies with distances smaller than 
~ 50 Mpc can still carry directional information. At these distances, most 
galaxies are located in the supergalactic plane. Most recently, the Auger 
experiment observed a first evidence for the correlation between the Veron- 
Cetty (V -C) catalo g of Active Galactic Nuclei and Cosmic Rays of the highest 
energies jAA+07a|. In the analysis, three parameters have been optimized in 



order to look for a correlation between the sources in the catalog and the 
events at the highest energies: The maximum angular separation ip between 
Cosmic Ray event and source, the maximum redshift ^max for which AGN 
in the catalog are still considered and the threshold energy Eth, giving the 
lowest energy at which Cosmic Ray events are still considered. These three 
parameters ensure that 
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Fig. 2. All particle Cosmic Ray spectru m. Data p oints coi ne from the experi ments a s 
listed in the bottom left corn er: Auger IyA+OtIi . HiRes |Hig02^ . AGAS A iY+95al |. 
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(a) magnetic deflection is considered to a certain amount [ip is larger than 
the point spread function of the detector), 

(b) only the local Universe is considered (the choice of ^^ax leaves out distant 
galaxies) , 

(c) only events with the highest energies are considered. Events of energies 
below Eth are discarded, since lower-energy events are more sensitive to 
magnetic fleld deflections. 
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Fig. 3. Skyplot of 27 Auger events at the highest energies {E > 5.7 • 10^^ eV, 
circles with radius ip = 3.1°) and the sources in the V-C catalog up to a redshift 
^max = 0.018 (red stars). Within this distance, 318 of 472 AGN in the V-C catalog 
are in the field of view of Auger. The blue outlines show the sensitivity of Auger. 
Darker color indicates larger relative exposure. The black, dashed line rep resents th e 



supergalactic plane. The closest AGN, Gen A, is marked in white. From AA"'"07a| |. 
Reprinted with permission from A A AS. 

The final parameter set is {ip , z^^ax, Eth) = (3.1°, 0.018, 57 EeV), and the sky 
plot is shown in Fig. 3. The maximum redshift of z^^x = 0.018 corresponds to 
a distance of 75 Mpc. The correlation is a first evidence that sources following 
the structure of the supergalacitc plane are the sources of the highest energy 
Cosmic Rays. The analysis does not allow for the direct identification of the 
sources for different reasons. Firstly, the catalog of AGN used in the analy- 
sis cannot considered to be complete and is partly inhomogeneous, including 
different types of AGN. Secondly, sources other than AGN, which also follow 
the local structure of the Universe, may be responsible for the emission of 
UHECRs rather than AGN themselves. Nonetheless, this result is a first step 
towards the identification of the sources of extragalactic Cosmic Rays. 



3.1.1.2 The Greisen Zatsepin Kuzmin cutoff Protons effectively lose 
energy by interactions with the CMB on their way to Earth: 



P ICMB 



A+ 

(ii; 



The energy loss length for pair production at ~ 10^^ eV is about ~ 500 Mpc, 
while it is only ~ 50 Mpc for the production of a A— resonance. The latter 
process is therefore responsible for a rapid decrease of the particle spectrum 
above Eqr > 5 • 10^^ eV, if the sources of Cosmic Rays are located at distances 
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larger than ~ 50 Mpc Gre66l . lKZ68l ]. This effect is named "Greisen Zatsepin 
Kuzmin (GZK) cutoff" . 

At the highest energies of Ecr > 5 ■ 10^^ eV, a discrepancy between the CR 
fluxes as observed by two experiments, AGASAH and HiRe^ was found. 
AGASA data are represented by stars and HiRes is shown as open boxes 
in Fig. 2. While AGASA detected several events above 5 ■ 10^^ eV, HiRes 
observed a decay of the spectrum. Given the prediction of the GZK cutoff, the 
result from HiRes is expected, while the AGASA events need to be explained 
by exotic phenomena. A surface array for the measur ement of the charged 
component of the showers was used by AGASA Y"'"95bl ] while HiRes HiR04 
consisted of telescopes measuring the emission of the showers' fluorescence 
light. At these energies, there are quite large uncertainties in the calibration 
of the spectra due to the low statistics and large systematic errors. Taking 
this into account, it is possible to interpr et the two results in a way which 
would still fit a single theory, see C"'"06bl ]. The Auger expe riment is being 
built to resolve the issue of the highest energy events, see e.g. [DA+Od l. With 
Auger, the two different techniques as used by AGASA and HiRes are 
combined in a hybrid array. This allows for the investigation of systematic 
uncertainties. The size of the array ensures sufficient statistics. Auger has 
now instrumented more than 85% of a 3000 km^ surface array and all four 
telescopes for fluorescence measurements are operating since M arch 200 7. The 
hybrid array is expected to be completed by the end of 2007 DA"*" 061 ]. First 



results from the AuGER array show that a component beyond the GZK cutoff 
can be excluded: a power law beh avior above E = 10^^'^ eV can be excluded 
with a significance of 6 a YA"'"07 . 



3.1.2 Production of ultra high-energy Cosmic Rays 

There are two basic scenarios for the production of charged Cosmic Rays at the 
highest energies {Ecr > 3 ■ 10^^ eV), referred to as the top-down and bottom- 
up scenarios. In the top-down scenario, UHECRs come from the decays of 
superheavy particles with masses ranging from 10^^ GeV up to the GUT scale, 
Mgut ~ 10^'' eV. Such decays result in protons at the highest energies up to 
Ep ~ 10^2 gv, and in this scenario, the GZK cutoff can be avoided, since the 
protons are produced in the Earth's vicinity. Furthermore, the proton signal 
is expected to be accompanied by a large flux of neutrinos and high-energy 
photons, see e.g. iBSOd . The absence of such signatures in combination with 
the confirmation of the GZK cutoff favors the production of protons in particle 
acceleration processes in distant sources, referred to as the bottom-up scenario. 



Akeno Giant Air Shower Array 
^ High Resolution Fly's Eye Detector 
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3.1.2.1 Bottom-up Astrophysical environments are often characterized 
by the colhsion of different plasmas. As an example, a remnants of supernova 
explosions are observed for typically more than 1000 years, as a consequence of 
the supernova shell being accelerated into the interstellar medium. A second 
example is the collision of two galaxies due to gravitational interaction. A 
shock front is produced, when a gas encounters other gas or a wall, with a 
velocity faster than any signal velocity. Such phenomena are not only observed 
in astrophysical environments, but also in other media, for instance in the 
atmosphere, such as supersonic movement of planes or bullets in air. The plane 
or bullet moves faster than the characteristic speed o f the raedium, the speed 
of sound, and produces a shock wave, the Mach cone MW84I . lMW85l . lMac98 . 
In astrophysical shocks, the characteristic speed of the plasma is the speed of 
magnetic waves. 

The concept of stochastic particle acceleration has for the first time been 



presented by Fermi [Fer49l . lFer54| . In astrophysical shocks, inhomogeneous 
magnetic fields are responsible for the acceleration of the charged particles. A 
single particle crosses the shock-front back and forth, gaining a constant frac- 
tion of energy per encounter with a magnetic-field inhomogeneity. The par- 
ticle leaves the shock-front region when carrying sufficient energy to escape. 
Stochastically, the acceleration of a whole population of particle results in a 
powerlaw-behavior as observed for the spectrum of ch arged Cosmic Ra ys. The 



theo ry of ac celeration was refined i n the 1 970s by Bell [BelTSal . lBel78bl | , Krym- 
skii Kry77l ]. Blandford fc Ostriker B078| and Axford, Leer & Skadron ALS78 
While Bell worked out a microscopic approach, in which individual particles 
are traced, Krymskii, Blandford & Ostriker and Axford worked on the macro- 
scopic description of astrophysical shocks, neglecting any individual movement 
of particles. One important consequence of acceleration theory is that the max- 
imum energy depends on t he ma gnetic field and the size of the acceleration 



region as derived by Hillas |Hil84 
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Ps ■ Z ■ B^G ■ Rkpc ■ 



(12) 



Here, -E^ax '■— -E'max/(10^^ eV) is the maximum energy which can be achieved, 
Ps = Vs/c is the shock velocity in terms of the speed of light and Z is the 
charge of the accelerated particle in units of the charge of the electron, e. 
Furthermore, B^^g '■= B/{1 fiG) is the magnetic field of the acceleration region 
in units of 1 fiG and Rkpc '■= kpc) is the size of the acceleration region 

in units of 1 kpc. The spectral index of the particle spectrum depends on 
the conditions of the astrophysical shock, i.e. the magnetic field strength, the 
orientation of the shock-front towards the magnetic field, the shock-front's 
velocity and the extension of the shock itself. 
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3.1.3 Primary spectra and radiation fields 

It is essential to study the correlation between accelerated primaries and sec- 
ondary radiation effects in order to connect photon observations to Cosmic 
Rays and neutrinos. In this section, the focus lies on the synchrotron emis- 
sion, since its observation reveals the spectral behavior of the primaries as 
discussed in the following. 



3.1.3.1 Total synchrotron power The total radiated power PiQt is pro- 
portional to with m as the mass of the particle. Since the electron-proton 
mass ratio is rrie/mp ~ 5 ■ 10~^, the radiated energy from electrons is a factor 
~ 10^^ higher than for protons. Protons only lose energy to synchrotron ra- 
diation for extremely high energies and large magnetic fields, since the power 
increases with the squared product of the external magnetic field B and par- 
ticle's energy Ep, Ptot oc {B • Ep)"^. Electrons undergo synchrotron losses at 
moderate energies already. 



3.1.3.2 Synchrotron radiation and the spectral shape For non-thermal 
spectra, the spectral index of the primary shock-accelerated particles can be 
expressed in terms of the synchrotron spectral index of a source. Electrons and 
protons follow the same distribution, i.e. the spectral index of the electrons is 
the same as for the protons, Up. The power P{uj) per unit frequency duj for a 
particle accelerated by an external magnetic field B can be written in terms of 
a fu nction F{x), which only depends on the dimensionless variable x := uo/ujc 



see 



RL79| for details): 

PMcxF(-). (13) 

The critical frequency of the synchrotron spectrum is given as 

3(7-8 sin a r, , ^, 

= 7' , (14) 

2mc 

representing a measure for the maximum frequency of acceleration for the 
particle spectrum. In this expression, q is the charge, m the mass and 7 is the 
boost factor of the accelerated particles. The latter can be expressed in terms 
of the energy, 

Ep = -fmc^. (15) 
Since shock-accelerated primaries follow a powerlaw distribution. 



dK 



dEp 



dEp oc Ep~"^ dEp , (16) 



the total radiated power Ptotioj) can be expressed as 



PtoMo: r[ P{u)E;"-dEp^ F{—)E;''''dEp. (17) 

J Er, J E-n (jJf 
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With Ep (X •y (X cj^^ (see Equ. (14) and Equ. (15)), 7 can be substituted for 
X := u/uci'y) and the total power can be written as 

Ptot OC cU-("^-^)/2 . n . 2;K-3)/2 _ (^g) 

J Xl 

Since the integral does not depend on a;, the frequency dependency is given 

as 

Pto, OC ^-("--1)/^ (19) 

The total synchrotron spectrum therefore follows a powerlaw with a spectral 
index s 

Ptot cx . (20) 

This leads to a linear correlation between synchrotron and particle spectral 
index, 

' = ^- ^21) 
The differential spectral index which will be used in the following is = —s—1 
with dN^/dE^ oc E^"^ . The flatness of the synchrotron spectra is limited 
by the theory of synchrotron radiation. Since a single electron produces a 
spectrum with = —2/3, the total synchrotron spectrum of an electron 
population cannot be flatter than < —2/3. This phenomenon is referred to 
as the line of deatl^ 

The connection between synchrotron and particle spectral index becomes im- 
portant, for instance, in the case of Gamma Ray Burst (GRB) spectra, which 
are typically explained by synchrotron emission of electrons. In some cases, 
the burst spectra are flatter than the maximum values, which indicates that 
other phenomena like absorption due to high optical depth need to play a role 
in the radiative processes in GRBs as well. 



3.1.3.3 Electron cooling The calculations above assume that the dy- 
namical timescale of the system is much shorter than the cooling time ofelec- 



trons due to radiation losses. This is called the slow cooling regime, see |Kar62 
In the fast cooling regime of long dynamical timescales compared to radiative 
cooling, the photon spectrum is flatter by 1/2, 

s = l. (22) 

For the prompt emission in GRBs, for example, the dynamical time scale is 
short compared to the cooling time and slow cooling has to be considered. 



^ In the integral representation, S oc the spectrum behaves as S oc u^^^. While 
in GRB physics, it is more common to use the differential representation, the integral 
form is more common to use in the case of AGN spectra. 
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3.1.3.4 Further radiation effects Note that the synchrotron spectrum 
is in many cases altered by further radiation effects. The synchrotron field can 
interact with the electron population, leading to the Inverse Compton (IC) 
effect, which boost photons to higher energies. This scenario is called Syn- 
chrotron Self Compton (SSC). Processes like optical depth effects, extinction 
by dust, pair production or bremsstrahlung can additionally alter the observed 
spectrum. When external photon fields are present. Inverse Compton scatter- 
ing does not need to be induced by the synchrotron field. Such a process is 
referred to External Inverse Compton (EC) scattering. At the highest ener- 
gies {E^ ~ TeV), the decay of 7r° particles resulting from proton-proton and 
proton-photon interactions can also dominate the spectrum. The last process 
competes with the SSC model. TeV emission in such h adronic models are re- 
ferred to as Proton- Induced Cascades (PIC), see RacOd ] for a discussion of the 
emission features. The question whether high-energy photon signals originate 
from hadronic (vr^-decays) or leptonic (SSC/EC) processes is one of the most 
striking these days. In some cases, TeV photon emission can also be explained 
by proton synchrotron radiation. 



3.2 Sources of high- energy photons 



From the radiation processes described above, it is clear that the emission of 
high-energy photons is usually connected to the acceleration of electrons or 
protons in astrophysical sources. In this section, the most energetic sources in 
the sky are discussed with respect to their observation in photons, and the pos- 
sible contribution to the spectrum of charged Cosmic Rays which is displayed 
in Fig. 2. The observed photon spectra are essential for the prediction of neu- 
trino fluxes, since photons are the only messengers giving direct evidence on 
the properties of the sources. Galactic sources are supernova remnants (SNRs) 
as likely sources for the production of Cosmic Rays up to the knee as well as 
X-Ray Binaries (XRBs), in particular microquasars, and pulsa rs whic h are 
candidates for the production of Cosmic Rays above 100 TeV |Gai90| . The 
most energetic, extragalactic sources are Active Galactic Nuclei (AGN) as 
permanent sources in the sky as well as Gamma Ray Bursts (GRBs) as tran- 
sient eruptions. Considering the abundance of the different source types and 
their individual electromagnetic output, leads to the expectation that galactic 
sources can produce the Cosmic Ray spectrum up to the ankle while extra- 
galactic sources are responsible for the CR flux above the ankle. The power of 
electromagnetic output mirrors the power in Cosmic Rays, since electromag- 
netic radiation originates from the charged particles in the source. Table 1 
lists source classes with their intrinsic luminosity and possible contribution to 
the Cosmic Ray spectrum. 
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Source class 


typical eni. 




life 


energy range 


Rcf 




outDnt 




time 






AL AO i lU 


bOUROEb 










lu erg/s 




iuuu yr 


10^° eV< Ep < lO^^ eV 


[GS64] 


SNR-wind 






innn vr 

J- \J\J\J jL 


10^° eV< Ep < lO^s eV 


[VB88] 


X-ray binaries 


1038 erg/s 




105 - lO'^ yr 


10^^ eV< Ep < lO^s eV 


fGai9nl 


Pulsars 


1037 gj-g/g 




10^ vr 


IQi^ eV< Ep < 1018 eV 


[Gai90l 


rL/X 1 KAGAL AO i lO 


oUUKOto 








Galaxy clusters 


~ 10^^ erg/s 




10^ yr 


3 • 10^8 eV< Ep < 10^1 eV 


fKRB97l 


AGN 




?/s 


lO'^ yr 


3 • 10^8 eV< Ep < 10^1 eV 


[BS87] 


GRBs 


1049 _ io5i erf 


?/s 


> 1 - 100 s 


3 • lO^s eV< Ep < 10^1 eV 


fVie95. WB971 



Table 1 

Source classes possibly contributing to the Cosmic Ray spectrum. The values for the 
typical electromagnetic emission and the life time are estimates based on the objects' 
observations. It should be noted that the luminosity and lifetime distributions can 
scatter, and there are objects within the given classes which have values deviating 
from the given typical ones In the case of GRBs, a lower limit is given, the actual 
value depends on whether or not GRB afterglow emission contributes to UHECRs. 
The energy ranges are estimates from the given references. 



3.2.1 Active Galactic Nuclei 

A class of galaxies with a particularly bright core was detected for the first 
time in 1962. An object, appearing star-like in the sky, showed extreme radio- 
emission features and could therefore not be classified as a star. The interpre- 
tation that this object, today known as 3C 273, was indeed a distant galaxy 
with a very bright core, was su ggeste d for the first time one year after the 



Sch63l ]. This class of objects was referred to 



detection by Maarten Schmidt 
as Quasi Stellar Objects (QSOs). 

Today, it is known that QSOs fit into the general classification scheme of Ac- 
tive Galactic Nuclei (AGN), objects which are believed to be powered by a 
rotating supermassive black hole in the center of the galaxy. A schematic view 
of the general picture of AGN is shown in Fig. 4. The core is "active" due 
to the accretion disk which forms around the central black hole and radiates 
strongly at optical frequencies. The disk is fed by matter from a dust torus. 
Perpendicular to the accretion disk, two relativistic jets are emitted, trans- 
porting matter in form of lobes. Knots and hot spots along the jets emit radio 
emission, leading to the strong observed radio signal of AGN. It is expected 
that these knots and hot spots represent shock environments in which particles 
are accelerated to hi gh ene rgies, in the case of hadrons up to proton energies 
of Ep ~ 10^1 eV, see |BS87] . In this section, a general classification scheme for 
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Fig. 4. Scheme of a cylindrically symmetric AGN shown in the r-z-plane, both axes 
logarithmicahy scaled to 1 pc. It is indica ted which objects are beheved to be seen 
from which direction. Figure after 

|zBoi]- 

AGN is presented as well as spectral and temporal properties of the sources. 



3.2.1.1 AGN unification scheme Three main criteria can be used for 
the unification scheme of Active Galactic Nuclei which is indicated schemati- 
cally in Fig. 5: 

(1) The activity of the source at radio wavelengths yields a division into radio 
loud and radio weak objects. About 90% of all AGN are radio weak and 
are usually hosted in spiral galaxies, while radio loud nuclei are located 
in the centers of elliptic galaxies. 

(2) The luminosity of the object is a further classification criterion. Radio 
weak sources are subdivided into optically strong and optically weak 
sources, which can be distinguished by considering the features of the 
emission lines. Optically strong sources usually lack narrow emission 
lines which are present in the optically weak case. Both source types ap- 
pear to have broad emission lines. Radio loud sources with extended jets 
(~ 100 kpc) are subdivided at radio wavelengths into low luminosity and 
high luminosity objects at a critical luminosity of = 2.5 ■ 10^^ W/Hz. 
The jets of compact objects such as GHz-Peaked Sources (GPS) and 
Compact Steep Sources (CSS) are believed to get stuck in matter. 

(3) The third classification criterion is the orientation of the AGN towards 
the observer. AGN are axisymmetric along the jet axis. In the branch 
of radio loud AGN, an object is classified as a blazar if one of the jets 
is pointed directly towards the observer. Flat Spectrum Radio Quasars 
(FSRQ) are the high luminosity population of the blazars while BL Lacs 
form the corresponding low luminosity population. So called Faranoff Ri- 
ley (FR) galaxies are being looked at from the side, so that jets and torus 
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Radio loud, elliptical C~10%) 



D 
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disturbed 
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logL >26.5 
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logL <26.5 
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Optically strong 
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broad 
lines 
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lines 
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FR-i 




Radio quiet 
Quasar 


[ssrq] 


Radio galaxy 




Radio gai. 





GPS/CSS 



[ hbl] ( TbT ) 



{ Blazars 



HBL High energy cutoff BL Lac 
LBL: Low energy cutoff BL Lac 
GPS: GHz peaked Sources/ CSS: Compact steep spectrum 
FSRQ: Flat Spectrum Radio Quasar 
SSRQ: Steep Spectrum Radio Quasar 
RIQ: Radio Intermediate Quasar 



Inclination angle 



luminosity 



tiost galaxy 



Fig. 5. AGN classification scheme, figure after AI"'"06bl |. 



are usually clearly visible. The high luminosity FR-II galaxies show a very 
strong radio emission at the outermost end of the jets, while the radio 
emission of the low luminosity class FR-I happens in knots throughout 
the jet. 

For radio weak AGN, the objects are called radio weak quasars in the op- 
tically strong case and Seyfert-I galaxies in the optically weak case when 
looked at the gap between jet and AGN torus. The radio weak equiva- 
lent to FR galaxies are Radio Intermediate Quasars (RIQ) and Seyfert-II 
galaxies, where the observer's view is directed towards the torus. 



3.2.1.2 Multiwavelength observations of AGN AGN have been ob- 
served in all frequency bands, ranging from radio observations up to TeV mea- 
surements. As an example, the lightcurves of the BL Lac object lES 1959+650 
have been investi gated i n a multiwavelength campaign between May 18, and 



August 14, 2002 KjlOJ]. The results are shown in Fig. 6. The bandpasses 



Lnopa 

shown are TeV emission as detected by Whipple (stars) and HEGRAqJ (cir- 
cles). X-ray emission as measured by RXTE0, optical emission in the Violett, 
Red and Infrared band as well as radio measurements from UMRAO0 at fre- 
quencies of 14.5 GHz and 4.8 GHz. As can be seen from the example, AGN 
are highly variable objects at all wavelengths. In multiwavelength campaigns 
as the presented one, the correlation between the temporal behavior in the 
different bandpasses is examined. Such campaigns are relevant to determine 
the origin of the radiation. In the case of the observation of lES 1959+650, a 



^ High Energy Gamma-Ray Astronomy 

^ Rossi X-ray Timing Explorer 
6 



University of Michigan Radio Astrophysical Observatory 
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so-called orphan flare was observed, representing a rapid increase in intensity 
only at TeV energies. Since SSC models necessarily predict the correlated emis- 
sion of TeV photons and X-rays, this scenario can be excluded for the observed 
flare, while an EC scenario is still possible. A hadronic scenario in which the 
TeV photons come from the decay of 7r° particles produced in proton-photon 
interactions does not require the connection between X-rays and high-energy 
photons. Orphan flares are therefore very interesting in the context of neutrino 
emission. 



Recent results from |Aha"'"06a| could restrict the temporal variability at TeV 
energies to less than 2 days. This implies scales of the order of the Schwarzschild 
radius of M 87, indicating that the TeV signal originates from the core of the 
object and not from the jets. 

Most recent ly, the d etection of 3C 279 by MAGIC at energies of > was 



announced [PMiQ7|. The detection this distant source {z = 0.538) was un- 
expected, since very high-energy photons are believed to be absorbed by the 
interaction with the infrared (IR) background, see also Section 3.2.1.3. The 
flux at energies of above > 100 GeV is therefore expected to be absorbed 
completely for sources with distances above z > 0.2. The spectral energy dis- 
tribution (SED) of AGN typically shows two main bumps, see Fig. 7, apart 
from dust radiation, which can also lead to a bump in the spectrum, see 



e.g. fc+si- 

The lower energy hump is believed to arise from synchrotron ra- 



diation of electrons. Different radiation processes can be responsible for the 
second hump at higher energies (GeV- TeV): in the Synchrotron Self Comp- 
ton scenario, the synchrotron photons are up-scattered to high energies by 
the primary electrons by the Inverse Compton effect. This implies the direct 
correlation of the two humps in the case of an intensity variation. A second 
possible scenario is the production of TeV photons in vr^-decays. vr" particles 
are produced in hadronic interactions with photon fields or with each other 
in the source. Such a hadronic scenario leads to the coincident production of 
high-energy neutrinos. It does not necessarily imply the coincident variation 
of the lower and higher energy hump, since the two emission signatures are 
not directly linked. A third component which can contribute at TeV energies 
is the synchrotron radiation of protons. The latter radiate in the case of high 
magnetic fields in the acceleration region. Depending on the energy range of 
the second hump in the SED, BL Lac objects can be divided into a further 
sub-class of High-peaked BL Lacs (HBLs), Low-peaked BL Lacs (LBLs) and 
Flat Spectrum Radio Quasars (FSRQs). If the peak occurs at TeV energies, 
sources are called HBLs, while they are referred to as LBLs at peak energies in 
the GeV range. The source is called FSRQ for even lower peak energies. The 
source classes can be identified by taking their radio emission at 5 G Hz as a 



measure. Figure 7 shows the SED for a sample of blazars presented in F"'"98 
The sample is divided into five sub-samples, selected via their radio luminos- 
ity at 5 GHz, \og{Lradio/eTg) = {< 41, 42 - 43,43 - 44,44 - 45, > 45}. The 



overlaid curves are analytic approximations as described in F"'"98 
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Fig. 6. Ligh tcurve of lES 1959+650 between May 18, and August 14, 2002. Figure 
from K"'"04l |. with permission from the author; Data (from the top): (a) TeV detec- 
tion by Whipple (stars) and HEGRA (circles); X-ray measurements from RXTE, 
(b): 10 keV, (c): spectral photon index between (3 — 25) keV; Optical results, (d), 
(e) and (f): V and / magnitudes from BOLTWOOD and R magnitude from BOLT- 
WOOD (Swiss crosses) and Abastunami (circles); (g) and (h): radio measurements 
at frequencies of 14.5 GHz and 4.8 GHz from UMRAO. 



3.2.1.3 TeV photons from AGN and the extragalactic background 
hght Since 2003, the second generation of Imaging Air Cherenkov Telescopes 
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log(i/) 



Fig. 7. Averaged SED for FSRQs, LBLs and HBLs. Reprint from |f+98I |. with 
permission from MNRAS. 



(lACTs) is taking data, with leading results from H.E.S.sJZl and MAGK 

and most recently also from VERITAS0. As of December 2007, 17 HBLs have 
been detected. Most recently, the first LBL, BL Lacer tae, has been identified 



at TeV energies by the MAGIC experiment AMjlOTaJ. In addition, M 87 as a 
misaligned blazar, or a FR-1 galaxy, is identified at energies of > 100 GeV. 
With 3C 279, the first FSRQ has now been identified by MAGIC at > 
100 GeV. This source is by far the most distant object detected at these high 
energies. The AGN identified at TeV energies and their spectral properties 
are listed in table 10. It is not possible to give a time-independent spectral 
parameterization of these variable TeV sources, since both the spectral index 
and the normalization vary with the intensity of the signal. In the table, intense 
flare events have been chosen. The values represent the maximum output per 



High Energy Stereoscopic System 
^ Major Atmospheric Gamma Imaging Cherenkov Telescope 
^ Very Energetic Radiation Imaging Telescope Array System 
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second observed so far. The parameterization for the spectra is taken as 



TeV 



TeV 



^^^^ ■ [rfev) 



(23) 



with the spectral parameters as measured at Earth. The signal of TeV pho- 
tons is (partially) absorbed at the highest energies by the Extragalactic Back- 
ground Light (EBL), comprised of the infrared background as the trace of 
star formatio n, the c osmic microwave background and other photon fields in 
the Universe |SdS92t . For a lower energy threshold of E^^'' ~ 100 GeV, it was 
expected that sources up to 2; ~ 0.2 can be detected. Thus, the detection of 
3C 279 aX z = 0.538 was unexpected and makes it necessary to reconsider 
the basic physics. Two assumptions go into the conclusion that sources far- 
ther away than z ~ 0.2 should be completely absorbed at > 100 GeV: 
(a) the presence of the IR background from star formation; (b) the intrin- 
sic particle spectra are not flatter than dNxev I dExev oc E^ly . With the 
detection of 3C 279, it is most likely that one of the assumptions does not 
hold. The intensity of the EBL at IR wavelengths is difficult to determine 
directly. However, star formation unavoidably leads to the production of at 
least a minimum component of IR radiation. The observation of such a distant 
source is difficult to explain by altering the EBL model within the standard 
models of cosmology and particle physics. Beyond the standard mo del of par - 



ticle physics, it is possible to introduce a new boson as done in AMR07 
Photons are in oscillation with a light boson of a mass m < 10^^'^ eV, lead- 
ing _to_tlie significant reduction of the ELB. An further approach is presented 
explaining the long attenuation length via the violation of the 



PMOO 



m 

Lorentz invariance at the highest energies. An explanation within the stan- 
dard model is that the photon spectra are indeed extremely flat intrinsically, 
which can be th e result of very flat primary particle spectra. The statement 
from Aha"'"06bl | that the p hoton s pectra can intrinsically not become flatter 
than E^ly is questioned in SBS07|. The basic assumption which produces the 
E~^-^ limit is that the primary particle spectra producing the Inverse Compton 
component at GeV- TeV photon energies is not flatter than ~ E~'^ as it was 



shown by |B098l . |K"'"00| . |OB02| | for relativistic, parallel shocks in the small 
angle scattering regime. It is pointed out in different papers, however, that 
the assu med shock conflguration is very specific and does not apply to most 



sources 



BarOd lMBQ07al . lMBQ07bj . In |SBS07 



it is shown that using larger 
scattering angles, very fiat spectra can be produced. This leads in turn to 
fiatter IC spectra, making it possible to extend the absorption cutoff to higher 
redshifts. To resolve this issue, further observations of more distant object are 
necessary, and the lowering of the threshold energy for photon detection to a 
few tens of GeV will also help. 

While lACTs are designed for the detection of point sources, surface water 
Cherenkov detectors instrument a large area with photomultipliers, being able 
to observe both photons, electrons and muons from electromagnetic showers. 
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This technique allows for the observation of 27r sr at once, however, with a 
relatively low sensitivity and a high-energy threshold. MILAGRO succeeded 
with the observation of Mkn 421 and the Crab nebula. In addition, diffuse 
emission regions in the Cygnus region were found, which is not possible with 
lACTs due to the small field of view. The next generatio n water Chere nkov 
telescopes, including the planned HAWCH^ experiment SSM05 . S"'"06 



are 



expected to give more information on diffuse emission due to a higher sensi- 
tivity. 



High Altitude Water Cherenkov 
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Source 


type z OiTeV FxeV first TeV 

detection 
10- (TeV s cm2) Ref 


M 87 


FR-I 0.004 -2.22 11.7 HEGRA 
fH.E.S.S. 2005^ [Aha+06al fAha+03al 


Mkn 421 


HBL 0.031 -2.14 311 Whipple 
(Whipple 2001) K+Ol] fPW+921 


Mkn 501 


HBL 0.034 -1.92 1001 Whipple 
fHEGRA 1997) [Aha+99l fQW+96l 


lES 2344+514 


HBL 0.044 -2.54 510 Whipple 
(Whipple 1995) [S+05bl fCW+981 


Mkn 180 


HBL 0.045 -3.3 12.0 MAGIC 
fMAGIC 2006) fAM+06al fAM+06al 


lES 1959+650 


HBL 0.047 -2.83 740 Tel. Array 
fHEGRA 2002) fAha+03bl fNis99l 


BL Lacertae 


LBL 0.069 -3.6 3.6 MAGIC 
fMAGIC 2006) fAM+07al fAM+07al 


PKS 0548-322 


HBL 0.069 n/a n/a H.E.S.S. 

fDH+071 


PKS 2005-489 


HBL 0.071 -4.0 2 H.E.S.S. 

fH.E.S.S. 03/04) [Aha+05al fAha+05al 


PKS 2005-489 


HBL 0.071 -4.0 2 H.E.S.S. 

fH.E.S.S. 03/04) [Aha+05al fAha+05al 


RGB J0152+017 


HBL 0.080 n/a n/a H.E.S.S. 

fNC+071 


H 1426+428 


HBL 0.129 -3.5 50.0 Whipple 
(Whipple 2001) fP+02l fHV+Oll 


lES 0229+200 


HBL 0.139 -2.5 6 H.E.S.S. 

fH.E.S.S. 05/06) [Aha+07al [DH+071 


H 2356-309 


HBL 0.165 -3.0 3.00 H.E.S.S. 

fH.E.S.S. 2004) [Aha+06cl fAha+06cl 
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Source 


type 


z 


OLTeV 


FTeV 

10-i3/(XeV s cm2) 


first TeV 
detection 
Ref 


lES 1218+304 


HBL 


0.182 


-3.0 

(MAGIC 2005) 


12.7 

fAM+06bl 


MAGIC 
fAM+06bl 


lES 1101-232 


HBL 


0.186 


-2.88 

(H.E.S.S. 2006) 


4.44 

fAha+07bl 


H.E.S.S. 

fAha+07bl 


lES 0347-121 


HBL 


0.188 


3.10 


4.52 

fAha+07cl 


H.E.S.S. 

fAha+07cl 


lES 1011-^496 


HBL 


0.212 


-4.0 

(MAGIG 2006) 


3.2 

fAM+07bl 


MAGIC 

fAM+07bl 


3C 279 


FSRQ 


0.538 


n/a 


n/a 


MAGIC 
fPM+071 


PG 1553-M13 


HBL 


? 


-4.2 

(MAGIG 05/06) 


2.1 

[AM+07cl 


H.E.S.S. 

fAha+06dl 

MAGIC 

fAM+07cl 



Table 2 

AGN with TeV photon emission. The spectrum is dNxeV / dE^ = FtcV • 
(E'.y/lTeV)"^"^ Sources are sorted by redshift. The redshift of PG 1553-1-113 is 
still unknown. The sources with no spectral index announced (n/a) are preliminary. 
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3.2.2 Gamma Ray Bursts 



Photon eruptions of unknown origin were detected in the 1960th by both 
American and Soviet mihtary satelhtes. While it was immediately clear that 
these events were not man made, but originated from outer space, the publi- 
cation of the first observat ion in 1 967 did not happen before 1973 in the case 
of the Vela Satellites KS073l | and only a few months later in the case of 
the Soviet Kqsmq s -461, a nd the American 0S0-7f"1 and IMP-d^ satel- 
lites ICD73l . IMGI74I ] . Systematic studies of these Gamma Ray Bursts 



(GRBs) were done with B ATSEHH on board of the CGRGjfJwhich was tak- 



ing data for 9 years, between April 1991 and June 2000 [Pjl99|. During that 
time, 2704 GRBs in the energy range of (20, 2000) keV were detected I. 
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Fig. 8. Distribution of tgo for 1234 GRBs in the BATSE 4B Catalog |P+99l . lBATn8l |. 
Events with tgo < 2 s are classified as short bursts, GRBs with tgo > 2 s are called 
long bursts. Figure courtesy NASA/BATSE and NASA Gamma- Ray As tronomy 
group at National Space Science and Technology Center (NSSTC) |NSS08l ]. 



Figure 8 shows the distribution of the duration of the bursts tgo BAT08l . lP"'"99 



defined in the way that 90% of the signal was received during that time. Two 
populations of bursts can be identified, classified as "short" (tgo < 2 s) and 
"long" (tgo > 2 s) bursts. The spatial distribution of GRBs in galactic coor- 
dinates as observed by the BATSE reveals an isotropic distribution with no 



Orbiting Solar Observatory-7 

Interplanetary Monitoring Platform-6 

Burst and Transient Source Experiment 
-"^^ Compton Gamma Ray Observatory 

In the following, the detection energy range of all quoted instruments is given in 
the notation (E'min) -E'max) heV for simplicity. 
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visible clustering in the galactic plane or anywhere else. This indicates an ex- 
tragalactic origin of the events. However, scenarios of a g alactic halo with so fa r 



unknown sources of GRBs were also proposed, see e.g. [H+941 . |PRR95l . IFS96 
A first indication of a cosmological origin was given by the spatially non- 
Euclidean distribution of the source luminosity. The final proof of the cosmo- 
logical distance of GRBs was possible i n 1997 by the first afterglow observation 



by the BeppoS Axl I satellite, see e.g. |WRM97l |. While the prompt emission is 



mainly detected in the keV-MeV band, the so-called afterglow continues until 
long after the prompt emission and is seen in basically all wavelength-bands, 
from the radio band up to GeV-energies. From the afterglow-observation, host 
galaxies can be identified, or absorption and emission lines can be measured to 
determine the redshift at which the GRB occurred. These redshifts are cosmo- 
logical, so that GRBs are known to happen outside of our Galaxy. The reason 
for the intense discussion of a galactic origin was that the photon fiuence 
in the (20, 2000) keV energy band as measured by BATSE scatters around 

~ 10"^ erg/cm^ . (24) 

The term "fiuence" is used here as opposed to flux, since the units are erg/cm^, 
while a flux is typically measured per area and time interval. For extragalactic 
distances, the total luminosity Lgrb of a GRB event lies around 

Lgrb ~ 10'' erg/s (25) 

for isotropic emission, but a little bit lower in the case of beamed emission 
favored currently. This tremendous output lies more than four orders of mag- 
nitude higher than the typical output of AGN, the most luminous permanent 
source class in the sky, 

Lagn 10'' - lO^^erg/s . (26) 

While GRBs emit only for a short time of a few seconds, AGN are active over 
long periods, so that the time integrated output is comparable, ~ 10^^ erg. In 
the past decades, different models have been developed to explain the output 
from GRBs and AGN. Many of the classical arguments for the description of 



GRB physics ar e borro wed from supernova remnants |Wol72l . lCox72l . ICS74 



or AGN models ReeTOj], since the phenomena are similar, although on differ- 
ent spacetime-scales. Both AGN and GRBs show a variable time structure, 
supernova remnants as well as long GRBs are produced in the explosion of 
stars. In all three objects, shock fronts are responsible for particle accelera- 
tion and therefore, the non-thermal electromagnetic spectrum can explained 
by synchrotron radiation of electrons. Inverse Compton scattering and also 
by proton-photon interactions. The favored model which is able to explain 



Beppo stands for Guiseppe Occhialini, and "SAX" is an acronym for Satellite 
per Astronomia X 
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most of the phenomena connected to a GRB is the fireball model, explaining 
the huge elect romagnetic emission by shock formation of relativistic plasma 
shells, 



see e.g. 



Pir99 



PirOSl . IZM04j . Alternatively, th e Cannonball model tries 



DR04. Dar06 . However 



to explain GRBs by colliding plasma balls, see e.g. 
the latter meets observations difficult to match the predictions in this model, 
e.g. the prediction of apparent motion for the radio emission, which is not 
observed for GRBs. Thus, only the fireball model will be discussed in more 
detail in the following paragraph. 



3.2.2.1 Fireball model A schematic view of the fireball model is shown 



in Fig. 9. It is based on the model of stellar outbursts as described in [Sed58 



and was enhanced including relativistic effects in order to match GRB obser- 
vations. The fireball model does not give any constraint on the progenitor. It 
yields a phenomenological description of the actual burst observations. The 
basic idea is that a large amount of mass is ejected within a short time inter- 
val by a central engine. The plasma is ejected successively in shells. At some 
point, the outer shells slow down and are caught by inner shells and a shock 
front is built up, accelerating electrons and baryons in the plasma up to high 
energies. While proton s can be accelerated basically loss-free up to energies as 
high as lO^i eV 



Vie95. WaxOO 



electrons lose their energy to synchrotron ra- 
diation, escaping from the shocks as soon as the region becomes optically thin. 
This is observed as prompt emission from GRBs. Those shocks resulting from 
collisions of shells are called internal shocks. So-called external shocks result 
from collisions of the shells with the interstellar medium leading to afterglow 
emission as described in the following paragraph. While the prompt emission 
occurs mainly at energies of > 100 keV, afterglow emission is observed in 
almost all wavele ngth bands. Reviews on the details of the underlying physics 



are given m e.g. 



Pir99, Pir05. ZM04 



3.2.2.2 GRB experiments after BATSE After the BATSE era, many 
GRB satellites were taking data, each covering a much smaller field of view 
and thus providing much less statistics than BATSE. Additionally, many 
of the satellites were not able to give directional information. For the de- 
termination of the GRB spatial origin, it is necessary to have at least three 
detectors. BATSE had four energy channels and could localize GRBs. The 
precision was between approximately 10° and a few degrees. Other satellites 
have only one or two instruments on board. To improve the localization of 
GRBs, the Interplanetary Network was created already in the 1970s as 
an interconnection of all GRB satellites. The currently active third Inter- 
planetary Network, IPN3, was formed with the launch of Ulysses in 



1990 H"'"92l . lUlyOSi ] . Since the beginning, more than 25 spacecraft missions 



have participated. The CGRO joined in 1991 when it was launched. Today, 
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Fig. 9. Schematic view of the firebah scenario. 

HETE- lil!] |S+05al. lHET08|. INT EGRAlH IWinOl Im+OS^ . RH ESSf . 
RHEOSl ^^iARsQDYSSEY |h+06I . [Mar08| . Ulysses |h+91 lUlvOSt . Konus 



Wind |M+05al . lKon08| and Swift form the IPN3. The most important GRB 



experiments with some of their individual properties are hsted in table 3. With 
information of more than two instruments, GRB positions of a accuracies up 
to several square-arcminutes can be reconstructed. Before the launch of Bep- 
poSax in 1996, this was the only possibility of arcminute precision measure- 
ments. HETE-II and INTEGRAL are able to localize GRBs without addi- 
tional information from IPN3. KoNUS is very sensitive to short GRBs and a 
catalog of 130 short GRBs is examined in the context of neutrino emission, 
see Section 7. 



High Energy Transient Explorer-II 

INTErnational Gamma-Ray Astrophysics Laboratory 

Ramaty High Energy Solar Spectroscopic Imager 
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GRB Sat. 


Launch-Demise 


FoV*** 


E-range 
[keV] 


Localiz. 
Precision 


Reference 
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r\K /I CidCi /I nvn 

uti/ iyby-ub/ ly /y 


b.i X b.i 
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1 r^alU / 1 , 
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~ 27r sr 


/on onnn^ 
[ZD, zUUUj 


degree 
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U4 / iyyD-U4 / zuuz 


U.D — i.o 


{ r\ 1 Q^^^ 
(U.i, oUUj 


arcmin 


IB^97l. 
loeDUol 


iViAHo (JU. 


U4 / zuui-now 


bz 


{oU, iUUUUj 




\ti UbI, 
iiViaruo| 


NRAR 


n9/19Qfi-09/90ni 




(^ inooo") 




fNF, A 081 


J-VL^iN U O V V . 


11/1 QQzL nnw 


27r sr 


I XUUUU 1 




1 iVUliUO 1 




no /onno n,^TTT 
Uz / zUUz-now 




(^o, zUUUUj 




IT} wpnsl 
1 xtrLii/Uo 1 


INTEGRAL* 


lO/2002-now 








fM+05bl 


CT3T 






(i5,oUUUJ 












/l^ 10000"! 


(XL (jlLLLLL 




"LJTPT^P TT* 

rlrj 1 Hz-ii 


1 n /9nnn nmir 


1.5 — 3 sr 


(^U.O, 4UU J 


arcinm 


FQ+O^qI 
lo uoai 


bWIFT 


ii/zUU4-now 








iuniub|. 


BAT 




2 sr 


(15, 150) 


arcmin 


fSwi08l 


UVOT** 




17' X 17' 


(170, 650) 


arcsec 




XRT 




23.6' X 23.6' 


(0.2, 10) 


arcsec 





Table 3 

Some GRB Satellites and the basic properties of the main instruments for the de- 
tection of prompt emission. FoV= Field of View. /* : Current member of IPN3. 
/** : the energy range for UVOT is given in terms of wavelength, in units of nano- 
meters (nm). /***: FoV for BATSE and KoSMOS-641 is actually Att sr, but here 
it is assumed that ~ half a hemisphere is occulted by the Earth. 
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The Swift satelhte was launched in November 2004 - for a review see 
e.g. 



[Chi06[ | and references therein. Swift is a dedicated GRB satellite with 
four instruments on board. The main purpose of the BAT0 detector is the 
discovery of prompt emission from GRBs. The main sensitivity is in the en- 
ergy range of (15, 150) keV and the field of view is about 2 sr. About 100 
GRBs per year are detected with BAT. The XRT0 covers an energy range 
of (0.2, 10) keV and serves afterglow observations. The UVOt[^ is an instru- 
ment for the detection of the optical afterglow at (170, 650) nm wavelengths. 

The advantage of a satellite carrying both prompt emission and afterglow 
instruments is that the afterglow can be followed almost starting from the 
prompt emission phase. This has already lead to the discovery of unexp ected 
temporal behavior directly after the prompt emission, see e.g. MesOq for a 
review. The basic features of the early afterglow are shown in Fig. 10. While the 
afterglow appearance at t ~ 10'* s after the prompt emission had been known 
before, the temporal behavior at earlier times was unexplored until the launch 
of Swift. Two main features are found. Firstly, a break in the temporal decay 
structure is observed. The decay index changes from at ~ —3 to ~ —0.5 
at around 100 s to 1000 s after the prompt emission and it changes back 
to the previously observed at ~ —1.3 behavior at 10"^ s to 10^ s. Secondly, 
X-ray flares are detected during the early afterglow, indicated by the dashed 
triangle in the curve. X-ray flares occur only in a fraction of the observed 
bursts. Different models to explain both the changes in the deca y index and 



the X-ray flares have been developed which are reviewed in e. g. |Mes06 



3.2.2.3 Observed Redshifts from GRBs and GRB progenitors The 

first afterglow observation by BeppoSax for GRB970228 also implied the first 
measurement of the redshift of a GRB. Between 1997 and November 2004 - the 
time of the launch of SwiFT - redshifts of 44 GRBs have been detected. Since 
the launch of Swift, 58 long and 6 short GRB redshifts were measured as of 
March 12 2007. Figure 11 shows the redshift distribution of those GRBs. The 
solid line represents Swift bursts while the dashed line shows pre-SwiFT 
measurements of redshifts. The mean of the distributions is shifted. Swift 
observations show a higher contribution of very distant bursts. One reason for 
this is the better sensitivity of Swift compared to pre-SwiFT instruments. 
Bursts at high redshifts have typically a weaker fluence and may have been 
missed by instruments of lower sensitivity. The higher statistics of bursts with 
measured redshifts relies on the possibility of early afterglow observation of 
the Swift instruments XRT and UVOT. 

It is known since 2003 that long GRBs are connected to super nova explosions 



of type Ic, which follow the death of Wolf-Rayet stars |M"'"03al | . Two scenarios 



^'^ Burst Alert Telescope 

X-Ray Telescope 
22 UV/Optical Telescope 
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of producing jets in exploding stars have been discussed: 



1) the cor e collapse from the center of a flattened, fast rotating star, see |Sak71 



M5176 



(2) 



the magneto-rotational core c ollapse with a n initial dipole magnetic field 
configuration as suggested in iBisTOl . iBisTli 



Short bursts have been proven in 2005 to originate from the merging of two 



neutro n stars o r a neutron star and a black hole in a binary system H^05 



V"'"05d . lG"'"06bl |. The differences in the observed events lie not only in the 
duration of the bursts, but are also seen in the redshift distribution. While 
long GRBs are most likely to follow the star formation rate and are located in 
starforming regions, short bursts happen in regions of rather low star formation 
rate and at small redshifts {z ~ 0.1). 

New results from the Swift satellite show, however, that this scheme is still 




Fig. 10. The figure shows the early afterglow flux decay with time. The scale is double 
logarithmic and the temporal behavior follows t"^* with differe nt deca y indices at for 
different time periods as indicated in the figure. Figure after [MesOd ] . After a steep 
decay of t~^, a shallow decay phase follows, ~ t~^'^. In this second phase, afterglow 
flares can occur in the burst. The observation of this early afterglow up to about 
10^ s after the burst has for the first time been observed by the Swift experiment. 
The phase at later times, t > 10"* s was known since BATSE times. A spectral 
behavior of t~^'^ is observed at first, with a break in the slope at ~ 10^ — 10^ s to 
a steeper behavior of ~ 
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Fig. 11. Observed bursts with measured redshifts. The sohd hue represents long 
bursts with redshifts from the SwiFT time, the dashed Une shows long GRBs in the 
pre-SwiFT era. There is a clear shift of the distributions. Bursts in the pre-SwiFT 
era (mostly BATSE measurements) have lower redshifts on average. The dotted 
line shows the distribution of short bursts with measured redshifts in the Swift 
era. There were no redshift detections of short bursts before the launch of Swift. 

too simple: there are exceptio nal bursts whic h do not fit into this scheme 
(e.g. GRB060218, GRB060614) |C+n6al . lG+n6j . Also, Swift does not see the 
strong di stinctio n between short hard and long soft bursts as it was observed by 
BATSE fSaSol. This indicates that the classification scheme is more complex 
than it can be determined yet and needs to be refined in the future. 



3.2.2.4 Classification of GRBs In order to explain the observation of 
the prompt emission in GRBs, a boost factor of F > 100 for the shock fronts 
is necessary. For lower boost factors, the shock region is optically thick to pair 
production processes. On the other hand, the boost factor must not exceed 
F ~ 1000, since proto ns woul d lose most of their energy due to synchrotron 
radiation in that case HaH02|. The peak energy of the GRB, Ep^ak is directly 
correlated to the boost factor and lies around Ep^ak ~ (100, 1000) keV for 
regular GRBs. Typically, a boost factor of F = 300 is assumed. The two main 
sub-classes of GRBs are long and short ones. The differences do not only lie 
in the duration of the events, but also in the hardness of the spectra. Short 
bursts typically have much harder spectra than long GRBs which is why they 
are usually referred to as Short Hard Bursts (SHBs). 
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A measure for the hardness is the ratio of soft to hard emission, 



S^(hard) 

hardness ratio := H := — — — . (27) 

D(sort) 

with S as the photon flux. In the case of BATSE, which had four energy 
channels, channel 3 with an energy range of (100, 300) keV and channel 4, 
(50, 100) keV were used to determine the hardness ratio. Bursts with tgo < 2 s 
are generally harder than long bursts (tgo > 2 s). It should be noted, though, 
that follow-up experiments like HETE-II, Ko NUS a nd Swift do not see 



equally hard short bursts. For more details, see |SaS06 



Detailed studies by HETE-II have shown that apart from the regul ar, lon g 



GRBs, there are bursts having peak energies in the X-ray regime |S"'"05a 
For the classification of events in terms of the energy band of emission, the 
hardness ratio for the flux at (30, 400) keV and (2, 30) keV was examined: 

_ ^(30, 400) 
^^^^^ - 5(2, 30) • ^^^^ 

Regular GRBs have H > 1. About 2/3 of the 45 HETE-II bursts have, 
however, H < 1. This class of GRBs peaking in the X-ray regime, has further 
been subdivided into X-Ray Rich bursts (XRRs) with X-ray and soft 7-ray 
emission, 0.3 < H < 1, and X-Ray Flashes (XRFs) with only X-ray emission, 
H < 0.3. The ratio between these three burst classes as observed by HETE-II 
is 

(GRB:XRR:XRF) = (1:1:1). (29) 

A schematic view of the classification of GRBs into long and short events 
and the further subdivision of the long events into regular GRBs, XRRs and 
XRFs is shown in Fig. 12. The r atio of long to short GRBs is indicated of 2 : 1 



as detected by BATSE P"'"99[ |. However, HETE-II rather detects a ratio 



of 3 : 1. The detection ratio should be considered as dependent on the instru- 



ment propertied I. For the analysis of the GRB-XRR-XRF relation, however. 



it is stated that not many GRBs are miss ed relative to XRR and XRF events 



due to any observational effects |S"'"05al |. so that a 1:1:1 ratio seems to be 
reasonable. 

One interpretation of the strong variation of the peak energy is based on 
variations in the baryonic load of the sh ocks, wh ich is an important parameter 



for the development of the boost factor |DCB99| . With a high baryonic load in 



the evolving jet, the system cannot be accelerated to high energies due to the 
high mass of the baryons. Thus, a this dirty fireballhas boost factors of an order 
of magnitude less than a regular GRB, F ~ 10. Since the peak energy evolves 



Swift for instance detects only one short GRB for 18 long ones. This is also due 
to the fact that Swift detects at relatively low energies E < 150 keV, while the 
emission for short GRB rather happens at higher energies as discussed before. 
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Fig. 12. Classification scheme of GRBs. 

with r, the typical values for dirty fireballs are Epeak ~ (I7 100) keV and XRRs 
or XRFs are observed. GRBs as typically observed by BATSE are classified 
as regular fireballs. Clean fireballs have a very low baryonic load and thus, a 
high gamma of around F ~ 1000, leading to peak energies of Ep^ak > MeV. 
The duration of such short high-energy bursts (SHE) would be small, since 
the lack of heavy baryons enables the jet to evolve more rapidly, so that 
tgo ~ 0.1 s. Such events have, however, not been observed yet. With the lau nch 
of glastS, the detection of such phenomena will be possible |GM99[ |. A 
summary of the different burst types and their basic parameters is given in 
table 4. The fact that most bursts occur as regular fireballs shows that a certain 
fraction of baryons needs to be present in the jet which inevitably leads to the 
production of neutrinos in proton-photon interactions. The question about the 
intensity of such a neutrino signal has still to be solved. 



parameter 


clean 


dirty 


regular 






(XRR/XRF) 


(GRB) 


r 


~ 3000 


~ 30 


~ 300 


Epeak [keV] 


> 1000 


1 - 100 


100 - 1000 


^90 [s] 


0.1 


~ 10 


~ 10 



Table 4 

Basic parameters of clean/dirty and regular fireballs. 



The classification scheme as presented here needs to be refined including a 
more detailed view on the matter in the future: today, there are several 
bursts which fall out of the scheme, e.g. GRB060218 and GRB060614, and 

Gamma-Ray Large Area Space Telescope 
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Swift se ems to detect a different sub-class of short bursts as compared to 
BATSE jSaS06 |. A summary of the obser vation of t he dif ferent sub-classes 



(XRF/GRB/short) with Swift is given in |Z+07j . In |Z+07| . the detection of 
the shallow decay of the early afterglow as recently observed in many Swift 
bursts as energy injection features. Such a scenario disfavors the interpretation 
of XRFs as low— r GRBs. Swift will help to improve the current classifica- 
tion scheme. For now, the scheme as presented above is still useful, since it 
can describe the majority of bursts. 



3.2.3 Galactic sources 

The electromagnetic output of galactic sources of non-thermal emission can 
be used to estimate their contribution to the Cosmic Ray spectrum. 

The Cosmic Ray luminosity can be written as 

LcR = jE{E„nn) (30) 

assuming the production of Cosmic Rays in the Milky Way with a residence 
time Tij ~ 6 ■ 10^ yr in the volume of the galactic disk Vgd ~ 10^^ cm^. Here, 

jEiE^in) = [ ^EpdEp (31) 

is the Cosmic Ray energy flux. The Cosmic Ray luminosity depends on the 
minimum energy which is produced. The calculation can only be valid at ener- 
gies below the ankle, Ecr < 3 ■ 10^^ eV: at the highest energies 
{EcR > 3 ■ 10^^ eV), the observed spectrum is too isotropic to be of galactic 
origin. Figure 13 shows the Cosmic Ray luminosity versus minimum energy. 
Three potential source classes are indicated as filled/hatched areas below the 
curve. In the following, the different source classes will be reviewed with their 
potential contribution to the Cosmic Ray luminosity. Here, it is assumed that 
the luminosity in Cosmic Rays produced by a source class must be less than 
the total electromagnetic output of the same class, L^"^^ < Lf^^'^. 



3.2.3.1 Supernova remnants The CR spectrum at energies below the 
knee is commonly believed to be produced by the shock fronts in expanding 
shells of supernova remnants (blue filled area). When a star perishes in a 
supernova (SN) explosion, the emitted material encounters the interstellar 
medium (ISM), building a shock front of typical velocities Vs ~ 10^ m/s. 
With the typical rate of SN explosions in a galaxy, dNsN/di ~ 3/(100 yr), 
and a mean ejected mass of Msn ~ 10 Mq per SN, the supernova remnant's 
shock front is active for about AIsnr ~ 1000 yr. The luminosity of a single 
SNR, 

LsNR ~ 10^^ erg/s, can then be converted into the total luminosity of 



38 
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Fig. 13. Cosmic Ray luminosity from the Milky Way and possible sources given their 
electromagnetic output. The upper bound of the curve indicates the Cosmic Ray lu- 
minosity for energies Ep > Emin- Cosmic Rays above the ankle, i.e. Ep > 3 - 10^^ eV, 
cannot have a galactic origin, since the observed signal is too isotropic. Since su- 
pernova remnants in our Galaxy can produce a luminosity of Lsnb. ~ 10^^ erg/s, 
they can be responsible for Cosmic Ray production starting at the lowest energies, 
i.e. Ep ^ 1 GeV (blue filled area). X-ray binaries (Lxrb ~ 10"^^ erg/s) and pul- 
sars (^Lp2ilsar ^ 

10^7 erg/s) can only be responsible for the Cosmic Ray flux above 
10^'^ GeV (XRBs, red hatched area, lines) resp. lO^-S GeV (pulsars, green hatched 
area, circles). Alternatively, supernova winds can be responsible for the contribu- 
tion above the knee, since higher energies can be reached in winds as opposed to 
SN explosions into the interstellar medium. 



SNRs in the Milky Way, 



tSNR ^ 
^tot ~ 



SN 



dt 



^SNR 



■ AtsNR ~ 3 ■ 10^" erg/s . 



(32) 



Given the integral luminosity required for the production of Cosmic Rays at 



Ep> 10^ eV, 



LcRiEp > 1 GeV) ^ 10^^ erg/s 



(33) 



SNRs are good candidates for the production of Cosmic Rays. It is, however, 
difficult to explain the break in the spectrum at Ep ~ 10^^ eV. One possibility 
is that leakage of particles out of the Milky Way becomes important, leaving 
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only heavy elements at the higher energies. This can lead to a steepening of the 
spectrum. Another possibility is that SN explosions into their own winds can 
be able to accelerate particles to higher energies, since higher mass numbers 
(Helium up to iron) are produced. SN explosions of type lb and Ic lose their 
hydrogen (for Ic also the helium) envelope before collapsing. This leads to a 
higher density of particles when the shock forms and thus to different shock 
conditions. Regular SNRs can in this scenario produce CRs up to the knee 
and SNR- Winds are responsible for the spectrum between the knee and the 
ankle. 



3.2.3.2 Pulsars, X-ray binaries and microquasars As an alternative 
explanation for the contribution above the knee, systems including neutron 
stars or black holes are considered. 

Neutron stars can be observed due to their emission of electromagnetic radi- 
ation along the magnetic field axis. Since the rotational axis of the objects 
does not align with the magnetic field axis, they are observed as pulsars: the 
emission is only seen when the particle jet points towards Earth. Pulsars have 
periodic signals ranging from several seconds down to milliseconds. The Crab 
as the most prominent, since most luminous, example is a millisecond pulsar. 
It is a neutron star which was presumably produced in a SN explosion ob- 



served on July 04, 1054 [fiJ. The pulsar wind nebula has been observed at all 
wavelengths. From radio up to X-ray energies, the supernova remnant is seen, 
while the pulsar itself is visible at X-ray and higher energies. The observed 
TeV signal from the Crab is non-the rmal, an i ndication for particle accelera- 



tion in a shock environment, see e.g. Aha^06el |. The main reason why pulsars 



are good candidates for particle acceleration are the very high magnetic fields 
of around B ~ 10^^ G. Pulsars have spin-down luminosities around 10^^ erg/s 
and can therefore not be responsible for the CR flux below the knee. They 
can, however, contribute to the region between knee and ankle (see Fig. 13, 
green hatched area, circles). 

Magnetars, which are observed as Anomalous X-ray Pulsars (AXPs) or Soft 
Gamma Repeaters (SGRs), represent pulsars with even higher magnetic fields 
of 5 ~ 10^^ G. For a summary on AXPs and SGRs see e.g. |WT04 |. Five 



SGR candidates have been observed in the Milky Way so far. These reveal 
themselves by randomly emitting 7 radiation from time to time. The smaller 
eruptions have usually thermal spectra, while giant bursts of non-thermal 
emission are observed from time to time. Famous events are the outbursts of 
SGR 1806-20 on January 7, 1979 and on December 27, 2004 as well as the 
giant emission of SGR 1900+14 on August 27, 1998. Flare luminosities of SGRs 
range from ~ 10'^'^ erg/s up to ~ 10^^ erg/s. AXPs show similar phenomena, 
but with emission at X-ray energies rather than 7 rays. The original definition 



For a summary of historical SNe including SN 1054, see 



40 



of AXPs was the steady emission of X-rays. About half of the AXPs today 
are known to be variable in X-rays, though. Eight AXP candidates have been 
observed in the Milky Way so far. The eruptions from magnetars are believed 
to come from star quakes, exciting the surface of the magnetar which leads to 
the emission of high-energy radiation. The energy output from AXPs lies at 
around 10^^ — 10^^ erg/s. 

Binary systems including a neutron star or a black hole are good candidates 
for shock acceleration as well. X-ray binaries consist of a compact object 
and companion star. In the case of Low-Mass X-ray Binaries (LMXBs), the 
non-degenerate star has about a few solar masses. High Mass X-ray binaries 
(HMXBs) are fed by a blue star (0/B star). Once the companion exceeds the 
Roche volume of the binary system, it starts to feed the compact object with 
matter. For HMXBs, accretion can also occur through accretion from stellar 
winds. The neutron star or black hole in turn emits the gained energy in form 
of X-rays and sometimes in a jet along the magnetic axis. Such systems can 
lead to particle acceleration up to the ankle at most. The typical electromag- 
netic energy of such X-ray binary systems is ~ 10'^^ erg/s. This is much lower 
than the total energy required to explain the CR luminosity for Ep > 10^^ eV. 
But considering only events at higher energies, i.e. Ep > 10^^'^ eV leads to a 
much lower CR luminosity due to the steeply falling spectrum (see Fig. 13, 
red hatched area, lines). For a summary of the XRB and CR connection, see 
e.g. 



Gai90 



A special case of an X-ray binary is the so-called microquasar, which are 
X-ray binary systems with photon emission along two co-linear jets. S uch a 



phenomenon has for instance been observed in GRS 1915+105 [MR94I ] . It is 
believed to be caused by instabilities in the accretion disk of the system. A 
burst is caused every time a particularly large amount of matter is accreted 
from the accompanying star. The development of the burst can be traced by 
the observation of the electromagnetic emission. Enhanced X-ray emission is 
seen close to the accretion disk, and radio to optical emission can be obse rved 



along the jet. For a recent summary of microquasar physics, see |Mir06a . 

For binary pulsar systems, a heavy Be-star or O star is accompanied by a neu- 
tron star. Periodic emission of high-energy photons can be observed in this 
case: Be-stars are massive stars with a circumstellar disk. The path of the neu- 
tron star around the Be-star is highly elliptic. This can lead to the accretion 
of matter by the neutron star from the massive star only near th e perias - 



tron, the point of closest distance of the two stars as discussed in |Mir06b 
A periodic emission of TeV photons was observed from three systems in the 
Milky Way, IS 5039 and P SR B1259-63 were detected by the H.E.S .S. experi - 



ment 



Aha+06ilAha+05b| . while LS I 61+303 was seen by MAGIC |AM+06d . 
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Fig. 14. Cross section for proton photon interactions in the center of mass system. 
Dotted hnes represent baryonic resonances. Direct pion production processes are 
shown as the dashed hne, while non-diffractive and diffractive multipion production 
are shown by the dash-dotted, upper and the dash-three-dotted, lower line, respec- 
tiv ely. Mo re information on the com parison between the data and the model is given 



m 



M+0d |. Reprinted from [m+00| |. with permission from Elsevier. 



3.3 High-energy neutrinos from astrophysical sources 

Neutrinos are produced in astrophysical shock fronts in proton-photon and/or 
proton-proton interactions via pion-production. The dominant channels are 



pn^ , fraction 2/3 
P7^A+^<( ' (34) 

r2 7r+ , fraction 1/3 

pp7r° , fraction 2/3 
pp — ^{ (35) 
pmr^ , fraction 1/3. 

The same processes occur for incident neutrons instead of protons, leading to 
the production of tt" par ticles. At higher energies, kaons can also contribute 
to the spectrum RacOd ] . Higher order processes are usually referred to as 



multipion production processes. The total cross section for proton photon 
interactions in the center of mass system is shown in Fig. 14. 
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While the resulting neutrons are likely to interact before decaying, charged 
pions decay and produce neutrinos, 



IT ^ /i U^^e Z/e Z/^ Z/^ . 



(36) 
(37) 



3.3.1 Neutrino oscillations 

Assuming that pions of negative and positive charge occur equally, the pro- 
duction flavor ratio of neutrinos at the source is 



[V^ : V,. 



;i:2:0). 



(3J 



This implies that tau neutrinos are not produced in astrophysical sources. 
However, since neutrinos have non-vanishing mass eigenvalues, they oscillate 
on their way to Earth. Neutrino fiavor-eigenstates >, a = e,fi,T, are 
connected to neutrino mass-eigenstates {uj >, j = 1,2,3 via 



i=l,2,3 



(39) 



Here, Ej is the eigenvalue for the Hamiltonian operator in vacuum, and thus 
the energy of the mass-eigenstate. The mixing matrix {Uaj) depends on three 
angles, 9i2, 6*13, and 6*23, and a phase 6. The matrix is given as 



U 



C12 Ci3 
~Si2 C23 

S12 S23 - C12 C23 sise 



■S12 Ci3 

iS 

C12 S23 Sl3 e C12 C23 
<5 



■5l2S23'Sl3e 



Sl3e 
■§23 Ci3 



i5 



-C12 S23 - S12 C23 Si3e'" C23 Ci3 



(40) 



Here, Cjk := cos{9jk) and Sjk := sm{9jiS), j = 1,2 and k = 1,2,3. The proba- 
bility for a neutrino to oscillate from a flavor state a to a flavor state /5 in a 
time t starting from the emission of the neutrino at the source, t = 0, is given 
as 



p 



\< i^p{t)W{t = Q) >| 



- 4 ^ Uai UfSi Uaj Upj SW? ( 



4^. , 



(41) 
(42) 



for units of c = h = 1. Here, Smfj is the difference of the squared masses of 
two flavors i and j. While solar neutrino detection basically provides mea- 
surements for 6*12, i.e. the oscillation from electron neutrinos to muon neutri- 
nos, atmospheric neutrinos provide 623, where muon neutrinos oscillate into 
tau neutrinos. The third angle ^13 appears to be the smallest and limits are 
pr ovided by reacto r experiments. The parameter set is approximately given 



as 



AS+041 . lAK+04 
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~ vr/6 
6*23 ~7r/4 

^13^0. 



(43) 
(44) 
(45) 



Thus, the mixing matrix can be written as 



U 







/ V3 1 

2 2 

]_ J_ 

2^2 2^2 V2 

1 73 1 

V 2^2 2^2 V2 / 



(46) 



For a non-monochromatic neutrino beam, the probabihty has to be averaged 
over the energy spectrum. Therefore, the sin terr a in th e oscillation proba- 
bility is averaged to ~ 0.5 for large distances L LP95I ]. Therefore, the os- 
cillation probability is independent of time for distances larger than the size 
of the solar system. The probability matrix for a neutrino flavor vector of 

^^souTce^ ^source ^^source^ changing tO a flaVOr vector ^lyEarth^^Earth^^Earth-^ -g 



given as 



V 



10 

4 
4 



At the source, the flavor ratio is 



4 
7 
7 



/ ^Earth \ 



V, 



V. 



Earth 
Earth 



(47) 



source 
source 



I 



Applying Equ. (47), the flavor vector approaches 



V' 



.Earth 

e 

. .Earth 
.Earth 



J 



(4J 



(49) 



for path lengths exceeding the size of the solar system. 

Thus, equal numbers of electron, muon and tau neutrinos are expected to be 
observed at Earth. 



3.3.2 Normalization of the spectrum 

Given that the existence of high-energetic protons in astrophysical environ- 
ments is bound to neutrino production, the observed flux of charged Cosmic 



44 



Rays can be used to esti mate t he neutrino intensity expected from the sources 
of Cosmic Rays, see e.g. |Hal06l ] and references therein. The Cosmic Ray energy 
density above the ankle j£;(-E'min = 3 ■ 10^^ eV) = JE^in=3-w^^'=v dNp/ dEp Ep dEp, 
is expected to be proportional to the neutrino energy density, 

jE{Emm) = Xu/p ■ J E^ dEy . (50) 

The constant of proportionality x^jp depends on the optical thickness of the 
source to proton-photon interactions and on the fraction of energy transferred 
to the charged pion. The fraction of energy going into the neutron does not 
contribute to neutrino production. The average energy which is transfered 
from the proton to the neutrino derives from the fraction of energy going 
into the charged pion, (xp^j^) ~ 0.2. The four leptons resulting from the 
pion decay are assumed to carry an equal amount of energy. Consequently, 
each neutrino which is produced carries 1/4 of the pion's original energy. The 
relation between neutrino and proton energy is therefore Ei^ = Ep/20. 

Due to the second channel in which 7r° particles are produced, the neutrino 
flux is expected to be accompanied by a high-energy photon flux: 

7r°^77. (51) 

The resulting photons are produced at > TeV energies. Thus, optically thin 
sources emit TeV photons in coincidence with high-energy neutrinos and the 
energy densities are proportional: 

/ dE^ = x,/^ ■ [ ^E,dE,. (52) 

The constant of proportionality x^/-y depends on the fraction of energy going 
into pion production. For optically thin sources, in the case of pp interac- 
tions, 1/3 of the proton energy goes into each pion flavor and the energy in 
^fi +'i^ti corresponds to the energy in photons, Xu/^y ~ 1. For p'-f interactions, 
Xi^/y ~ 1/4. If the source of neutrino emission is optically thick to TeV pho- 
tons, the TeV signal avalanches to lower energies until it reaches a level at 
which it can escape the source. Thus, sources of GeV to sub-MeV emission 
can be sources of neutrino production as well, assuming an optically thick 
environment. Equation (52) holds with different values for x^f-y and modified 
integration limits. 

In conclusion, the same sources which have been described as the potential 
origin of Cosmic Rays and which show signs of non-thermal photon emission 
are selected when discussing neutrino emission. In the following sections, the 
possibility of neutrino emission from galactic sources like supernova remnants, 
pulsars and X-ray binaries will be mentioned, as well as the predictions of neu- 
trino production in the decay or interaction of exotic particles like WIMPs or 
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monopoles. Furthermore, extragalactic neutrino emission, i.e. neutrinos from 
AGN, starbursts and GRBs is discussed. 
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4 High-energy neutrino detection methods 

When a neutrino ui = u^, Ufj,, Vr interacts with a nucleon via charged current 
interactions, a lepton / = e, r is produced, 

Nvi^Xl. (53) 

Here, X indicates the hadronic product of the interaction, leading to a hadronic 
cascade. Charged particles in a medium emit Cherenkov hght if traveling faster 
than the spe ed of l ight in the same medium. The Cherenkov effect is described 



in detail in | JacOGi ] . The blue Cherenkov light can be detected by photomul- 
tiplier tubes and both the incidental direction and the neutrino energy are 
reconstructible. The signatures for electrons, muons and tauons are very dif- 
ferent. Electrons and tauons produce electromagnetic and hadronic cascades 
within the detector. About 20% of the total energy go into the hadronic cas- 
cade, which arises from the nuclear recoil. The remaining 80% of the energy 
are carried by the electromagnetic cascade, produced by the charged lepton 



interacting with the electrons of the medium |Kle04j ] . In contrast, muons only 
undergo radiation losses and leave track-like signatures. The track of the muon 
can be reconstructed, since the Cherenkov signal is emitted in form of a di- 
rected cone. Tauons have very distinct signatures. A first cascade is produced 
in the first Vr N interaction, where a r and a hadronic cascade is produced. The 
r in turn decays, producing a second cascade in which a further r is generated. 
Thus, tau neutrinos have the unique property of regenerating themselves. 

In the case of neutral current interactions, 

Nui^Xui (54) 

and a hadronic cascade X can be observed. 
4 ■ 1 Signatures 

The primary detection technique for Cherenkov arrays is the measurement of 
neutrino-induced muon tracks in water or ice. To guarantee that the observed 
muons are neutrino-induced, the Earth is used as a filter. While muons pro- 
duced in the atmosphere are absorbed by the Earth, neutrinos traverse the 
Earth and the signature is unique. After the filtering of atmospheric muons, 
the remaining signal mainly consists of neutrinos produced in hadronic showers 
in the atmosphere. Various analysis methods have been developed to separate 
this background from a potential signal from extraterrestrial sources. In ad- 
dition to the neutrino-induced lepton, a cascade produced at the interaction 
vertex contributes to the signal. 

While the muon leaves an optical track, the cascade events produce a ra- 
dio and an acoustic signal in addition to the optical emission. The Landau- 
Pomeranchuk-Migdal (LPM) effect becomes important at the highest energies. 
El, > 10^^ eV, where the cross sections for electron bremsstrahung and photon 
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electron-positron pair production are reduced. As deduced in Kle04l ] . the elec- 
tromagnetic shower lengthens to about 100 m at 10^° eV. Most of the signal 
is radiated at the end of the cascade for each of the three detection meth- 
ods - optical, radio and acoustic. This is why a two-component signal can be 
expected at the highest energies. 

In the case of tau neutrinos, instead of a single cascade, multiple cascades are 
produced due to the regeneration of the tau. With IceCube, so-called double- 
bang signatures are expected, where two cascades are observed within the 
instrumented volume. If one of the cascades happens outside the instrumented 
volume, the signature is called lollipop. A single cascade is observed together 
with a muon bundle, so that the cascade finally looks like a lollipop. The lower 
energy detection threshold is given by the fact that the cascade-like events 
need to be larger than the spacing of the IceCube strings, which is about 
125 m. The assumption that the extension of the cascade event needs to be 
larg er than about ~ 200 m leads to a lower detection threshold of -Emin ~ PeV, 



see 



B+03 



The reproduction of tauons renders possible the signal propagation through 
the Earth. Therefore, tau neutrino cascades are suitable for the investigation 
of horizontal or Earth skimming neutrino air showers. The cascades are re- 
produced within the Earth until the shower comes out of the Earth and can 
be detected by either an Imaging Air Cherenkov Telescope or an air-shower 
surface array. Regular hadronic or electromagnetic air showers do not traverse 
the high column density and the neutrino sign ature is u nique. This technique 



is for instance used by the Auger experiment AA+07bl ] in th e case of s urface 



arrays and by the MAGIC experiment in the case of lACTs BM+07al . 



4.2 Atmospheric neutrino background 

The main background in the search for extraterrestrial neutrinos are atmo- 
spheric muons and neutrinos. When primary CRs interact in the upper atmo- 
sphere, a shower develops. If the primary is a hadron, e.g. a proton, secondary 
mesons and muons are produced while the shower propagates through the at- 
mosphere. Secondary mesons like charged pions and kaons produce neutrinos 
when decaying. The neutrino spectrum is close to an behavior in the 

energy range of 10^ GeV< < 10^. This is about one power steeper than 
the primary spectrum, ~ i?"^ '', since a considerable fraction of pions and 
kaons (r ~ 10~^ s) interact again before decaying. Neutrinos from pion and 
kaon decays are usually referred to as conventional atmospheric neutrinos as 
opposed to prompt atmospheric neutrinos from heavier hadrons. 
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4-2.1 Conventional atmospheric neutrinos 

An anal ytic des cription of the muon neutrino spectrum at > 100 GeV is 
given in 



VolSd 



as 



dN 



dEudn 



0.0285 • ^-2.69. 



1 



+ 



l+GE^/E^{e) 
0.213 



1+IA4E^/E ±(9) 



100 < < 5.4 • 10^ 



(55) 



0.48 • ^^4.04. 



> 5.4 • 10^ GeV 



and the electron neutrino spectrum is given as 

0.0024 • E-'^-^^- 



dN 



dEydn 



0.05 

l+1.5E^/E^± (61) 



+ 



+ 



0.185 



l+lA4E^/Ej^o{e) 
UAEiW 



100 < 



< 3.7 • 10^ 



(56) 



0.0071 • 

[EK±{e)+?,.7EKo{e)] E^ > 3.7-105 GeV. 



Here, E^,, Ek± and Eko are angle distribution parameters, depending on the 
zenith angle 9. The factor is energy dependent as 



m = <o)^h{d)-\og{E,). 



(57) 



The values of i?,r, Ek±i E^o, a{6) and b{6) are listed in table 5. At lower 
energies, i^j^ < 100 GeV, the description is given numerically by H"'"95 . 



cos(^) 


1 


0.6 


0.4 


0.3 


0.2 


0.1 


0.05 





EAO) [GeV] 


121 


202 


298 


392 


572 


886 


1060 


1190 


EjMO) [GeV] 


897 


1500 


2190 


2900 


4220 


6540 


7820 


8760 


EKo{e) [GeV] 


194 


324 


473 


628 


915 


1410 


1690 


1890 


a{e) 


-1 


-0.355 


-0.687 


-0.619 


-0.384 


-0.095 


0.0 


0.083 


m 


0.0 


-0.23 


-0.01 


-0.007 


-0.09 


-0.165 


-0.186 


-0.215 



Table 5 



Angle distribution parameters for the atmospheric neutrino flux V0I80I ] 
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The ratio of neutrinos to anti- neutrinos is energy dependent. In |Wei93l ] it is 
given as 



RatmiK) = {dNj{dE,dn))J {dNj{dE,dn))- = 

1 E^<1 GeV 

1 + 0.3 ■ log {EjGeV) 1 GeV < E^ < 100 GeV (5J 

0.6 + 0.5 ■ log {EjGeV) 100 GeV < E^ < 1000 GeV 

1.2 + 0.3 ■ log {E^/GeV) E^ > 1000 GeV . 



The spectrum is one power steeper than the spectral behavior of the primaries, 
since the probabi lity tha t pions and kaons interact before decaying increases 
with the energy GaiQOj . The intensity of the spectrum increases with the 
angle of the neutrino towards the zenith angle. The larger the angle, the lower 
the density gradient of the atmosphere and therefore, more pion decays are 
possible for horizontal events as opposed to vertical events. 



4.2.2 Prompt atmospheric neutrinos 

When heavy hadrons with charm and beauty content are produced in the 
atmosphere, the resulting neutrino spectrum appears flatter than the one of 
conventional atmospheric neutrinos. Due to short lifetimes ( ~ 10~^^ s) these 
hadrons do not interact before decaying for energies up to E^, < 10^ GeV. 
Therefore, the spectrum of the prompt neutrinos follows the primary index, 
oc Ey~'^''^ . The exact normalization and spectral behavior is difficult to predict. 
In order to estimate the energy spectrum of neutrinos from hadrons with charm 
contents, it is necessary to simulate interaction conditions at extremely high 
energies {Ei, > 10^ GeV). Interaction models are based on accelerator data 
reaching particle energies up to i?^ ~ 10^ GeV at most. The second challenge 
is the behavior of interactions at small Bjorken— x. The Bjorken scale invariant 
X is given as 

^ = • (59) 

Here, is the negative four-momentum transfer, M is the nucleon mass 
of t he interaction and v is the energy transfer in the laboratory system, see 



e.g. |Sch97l . IHM84I ]. The scale invariant is physically constrained to < x < 1. 
Interactions with small values of x imply forward-scattered processes. Acceler- 
ator experiments only yields measurements for x > 10~^, since the beam pipe 
does not allow the detection of particles at smaller x. Therefore, extrapolations 
to high energies and small x are needed in order to calculate the atmospheric 
charm flux. Figure 15 shows the a tmospheric prediction including conven- 



tional neutrinos according to H+04| together with different predictions of the 
prompt atmospheric neutrino flux. At energies of around E^ ~ 10^ — 10^ GeV, 
the prompt neutrino flux becomes stronger than the signal from the conven- 
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tional neutrinos. The uncertainties in the models are larger than one order of 
magnitude. Large volume neutrino detectors can help to determine the flux of 
prompt atmospheric neutrinos and thereby improve the knowledge of particle 
interactions at the highest energies. 



4-3 General detection method 



An estimate of the detection rate of neutrino-induced muons/cascades R can 
be given by folding the expected neutrino flux at Earth with the probability 
of the detection of the neutrino, 

Pshadow 

{e,E,)—^dE,. (60) 

Here, P^^i{E^, -Emin) is the probability that a neutrino interacts with a nucleus 
to produce either a muon - / = /i - or an electromagnetic or hadronic cascade 
- / = cascade - which is detectable in a large volume neutrino detector. It can 
be written as 

r^-' da 

P,^i = Na / dEi—ri{EuE^,,), (61) 

where Na is Avogadro's constant. The range of the produced lepton within 
detection range ri is determined by the geometric detection volume. For the 
case of muon neutrinos, the resulting muon can be observed even if produced 
outside t he detec tor itself. The average muon energy loss rate in a medium is 



given as 



GHS95 



l^-^) = -<E,)-^{E,).E,. (62) 

The flrst term occurs due to ionization losses of the muon. Monte Carlo sim- 
ulations show that a can approximately be regarded as a constant with a 
value of a ~ 2 MeV/(g cm^) in media like rock. The second term is due to 
bremsstrahlung, pair production and nuclear interactions. The parameter (3 
is determined via Monte Carlo methods, with a value which is approximately 
constant in a given medium. For rock-like media, /? ~ 4 ■ 10~®/(g cm^). For 
a discussion of a and (3 in the co ntext of muon tracking in a medium using 



Monte Carlo methods, see RCOl 



The muon range is then given by the integral of the inverse of muon loss rate 
over the muon energy: 

1 /-^M 1 



'log^(^. (64) 



/Ej^in {dE/dX) Je^^'^ a + b- E 

1 a/b + 1 
b a/b + E\ 

The differential cross section da/dEi includes charged and neutral current 
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Fig. 15. Atmospheric neutrino flux models including the contribution from hadron 
decays with char m con tents. The flux of conventional atmospheric neutrinos is 
given by Honda [H"'"04l. The prorn pt neutrino flux predictions are models given 
in: Fiorentini_etal. [FNVoH, Ib+89I | {Naumov RQPM and Naunwv QGS Mh Mar- 
tin et al. |MRS03I | {Martin MRS, Martin KMS and Martin GBW); Costa |cos01 | 
{pQCD opt, pQCD pes, QGSM opt, QGSM pes, RQPM opt, RQPM pes). pQCD 
models are determined using pertubative QCD in next to leading order. QGSM 
models are half-empiric, using the quark gluon string model pased on non-per- 
tubative QCD calculations. The recombination quark parton model (RQPM) is 
phenomenology-based and non-pertubative. Models marked opt are optimistic pa- 
rameterizations, pes stands for pessimistic. Figure courtesy Kirsten Munich, see 
also jMunOTj l. 



52 



Fig. 16. Traveling distance of a neutrino through Earth. 

interactions for cascade events, da/dEcasc = da/dEcc + cr/dEj^c- In the 
case of a muon signal, only the charged current cross section contributes, 
da/dEcasc = da/dEcc- The energy threshold of the detector is given as -Emm- 
The cross section is determined by using the parton distributio n functio ns 



given by jCER| . For the calculations in this review, the model of lGRV92f 
chosen 



IS 



The probability that a neutrino is absorbed by the Earth due to neutrino- 
nucleon charged and neutral current interactions is given by the shadow factor, 

PshadUX) = exp(-X(^)/A) , (65) 

in which X is the column depth in units of cmwe and A is the mean path 
that the particle can survive without interacting. The mean free path can be 
written in terms of the total cross section atot as follows 

Na CTtot 

Subsequently, the shadow factor can be expressed as 

Pshadow{X) = exp{-NA (Ttot X) . (67) 

The column depth depends on the traveling distance through Earth. It is also 
dependent on the Earth's density layers it has to pass. 

In order to calculate the absorption probability the neutrino's path through 
the detector needs to be calculated. Assuming a neutrino is entering the de- 
tector under a zenith angle 6 (see figure 16) and applying the sine theorem to 
the triangle including 7 := 180° — 0, f3 and a, the relation 

sin7 ^ i^e /^^x 
sin/3 R(B-d ^ ^ 
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is found, i?® is the radius of the Earth while d is the depth at which the 
detector is located under the Earth's surface. The angle a is determined as a = 
180° —7-/3. Applying the sine theorem once more and using the equation (68) 
with the approximation d « R^, the distance d^, which the neutrino has to 
travel through Earth is given ae 



26 



du 



sm a 



sm7 



(69) 



The column depth X{6) is the product of the density of the Earth, p^, with 
the traveling distance of the neutrino through Earth: 



Xi9) = ■ dM ■ 



(70) 



Depending on the zenith angle, the neutrino has to travel through different 
Earth layers which h ave diff erent densities. A model of the density profile of 



the Earth is given in [DA81 



p{r)[g/cm? 



13.0885 - 8.8381 -x x < 0.192 

12.5815 - 1.2638 • x - 3.6426 • x^- 

-5.5281 • x^ 0.192 < x < 0.546 

7.9565 - 6.4761 • x + 5.5283 • x^- 
-3.0807 • x^ 



5.3197 - 1.4836 • x 

11.2494 - 8.0298 • x 

7.1089 - 3.8045 -x 

2.691 + 0.6924 • x 

2.9 

2.6 

1.02 



0.546 < X < 0.895 
0.895 < X < 0.906 
0.906 < X < 0.937 
0.937 < X < 0.965 
0.965 < X < 0.996 
0.996 < X < 0.998 
0.998 < X < 0.999 
X < 1 . 



(71) 



Here, x is given as x := r/R^ with R^ = 6378 km as the radius of the Earth. 

The energy dependence of the shadow factor Pshadow and interaction probabil- 
ity is shown in Fig. 17 for the case of electron neutrinos. Solid lines represent 



Note that this is only the case for sin 7 7^ 0. For sin 7 = and 6 = 0°, the traveling 
distance of the neutrino through Earth is equal to the depth of the detector under 
the Earth's surface (d): di,{0 = 0) = d. If, on the other hand, sin 7 = and 9 = 180°, 
the traveling distance is dy{6) = 2-R^—d since the neutrino arrives from the opposite 
direction. 
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neutrinos and dashed lines show anti-neutrinos. The shadow factor is close to 
unity for low energies, while the interaction probability increases with energy. 
The product of both determines folded with the neutrino spectrum gives the 
event rate for a given flux. 




Fig. 17. Probabilities as a function of energy for the case of electron neutrinos. The 
shadow factor is shown for a nadir angle of 9 = 110°. The nadir angle is measured 
from the South Pole towards northern directions. Higher inclination angles result in 
the shift of the cutoff towards lower energies. The interaction probability increases 
with energy. For anti-electron neutrinos, the Glashow resonance produces a peak at 
10^'^ GeV. The product of the two probabilities (thick lines) is a measure for the 
detection probabilities (see text). 

Three main signatures can be detected from neutrino-nucleon interactions. 
The main detection method is the optical detection of the Cherenkov light. 
Detectors in water and ice have been and are being built. The main differences 
between the detector performance in ice and water lie in the scattering and 
absorption lengths in water and ice. Ice absorbs only little of the light, water 
has better scattering properties. In water, there are two additional background 
components, bioluminescence from bacteria populations in the sea water, and 
light from decay. 

Neutrino-induced cascades can also be detected at radio wavelengths and as 
acoustic pulses. These cascades are detectable over long distances and there- 
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fore, radio and acoustic arrays do not need to be instrumented as densely as 
optical detectors. On the other hand, larger volumes need to be deployed, since 
the radio and acoustic detection methods become effective at high threshold 
energies, -Emin ~ 10^ — 10^ GeV, at which the expected neutrino flux is already 
very low. 



4-4 Optical detection in ice and water 



4.4.1 Ice 

The largest currently operating, completed high-energy neutrino detector is 



AM AND A-Ill I. The low neutrino-nucleon cross section demands the instru- 
mentation of large volumes to achieve a significant rate of events. AMANDA 
is located in the Antarctic ice, between 1500 m and 2000 m below the surface 
and covers a geometric volume of ~ 0.016 km^. The detector consists of 19 
strings with a total of 677 PhotoMultiplier Tubes (PMTs). AMANDA can 
detect neutrinos with energies of a lower threshold energy of -Emin > 50 GeV. 

AMANDA'S successor IceCube is currently being built around the 
AMANDA experiment in 1500 m to 2500 m depth. The final configuration 
will include a geometric volume of 1 km^ with 70 — 80 strings, carrying 60 
PMTs each. During the past four polar seasons (each year between October 
and February the following year), a total of 40 strings has already been de- 
ployed (1 string in 04/05, 8 strings in 05/06 and 13 strings in 06/07, 18 strings 
in 07/08). As half a cubic- kilometer is already instrumented, completion of the 
full cubic-kilometer is expected in 2011. Since the spacing between the PMTs 
and also between the strings is larger than in AMANDA, the threshold en- 
ergy is slightly higher, -Emin ~ 100 GeV. With the increased volume, higher 
energy events are easier to identify (the light deposition and the path lengths 
are increased with increasing energy) and better sensitivities can be achieved. 



4.4.2 Water 

The optical detection of neutrinos in water via the Cherenkov effect has al- 
ready been probed by lower energy detectors like ( Super) Kamiokande and 
SNO as discussed above. In the case of high-energy neutrinos, natural water 
reservoirs are used for detection. The first generation arra y Baikal detect s 



neutrino-induced cascades and also muons in Lake Baikal AB"'"06al . IWis06 
The advantage of this lake is that it is deep enough to deploy ~ 80 m long 
strings in depths of ~ 1100 m and that the detector can be reached in winter 
due to the thick ice layer covering the lake. Baikal has been taking data with 
the NT200 detector since 1998. In April 2005, the detector was enlarged to 
NT200-I-, enclosing 0.005 km'^ of water. The sensitivity to a diffuse flux for all 
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three neutrino flavors lies around ~ 2 ■ 10 ^ GeV cm ^ s ^ sr ^ AB+OGb 



Currently, different experim ents are under c onstruction in the Mediterranean 
sea. In the ANTARESH |AA+97l . lMon03[ | project , 300 m Ions strings are 
being deployed up to a depth of ~ 2500 m |AA+97I . IANTOS^ . As of Decem- 
ber 2007, 10 lines were deployed and the first muon tracks could already be 
identified. The final detector configuration will consist of 12 strings with an 
instrumented volume of ~ 0.03 km'^. A second array is being constructed 
near the coast of Italy, the NEMC@ detector |Sap06l . iNEMOsj . at a depth 
of ~ 3500 m. Such a high depth has two main advantages: the flux of at- 
mospheric muons from above is lower in the detector, and also the biolumi- 
nescence background is low at high depths. The first tower was deployed in 
2006. In a second phase, the deployment of one or more additional towers are 
planned. NESTORl^^l is a th ird experiment, being deployed near the coast of 
Greece at a depth of 3800 m [AN+Oa . iNESOSl j. The first phase of the detector 
was deployed in March 2003, and first results are pre sented in form of the 



measurement of horizontal atmospheric muons AN"'"05 



The primary aim of the three arrays described above is to explore different sites 
and techniques for neutrino detection in sea water in order to build a 1 km^ de - 



tector, choosing the optimal location and design, namely KM3NeT |KM308 



This detector will be complementary to IceCube, detecting mainly muon neu- 
trinos from the southern hemisphere. With both arrays, a view of the complete 
sky is possible. The view on the southern hemisphere enables KM3NeT the 
observation of the galactic center, which is only restrictedly visible from the 
South Pole. 



4.4-3 Analysis strategies 

In order to keep the analysis blinded to avoid experimenters bias, analysis 
cuts are optimized using off-source samples created by scrambling the right 
ascension of events or excluding the time window of transient emissions under 
investigation. Diffuse flux analyses are optimized on a low energy sample, 
where the signal is expected to be negligible. 

The search for steady point sources is possible in two ways: the selection of data 
according to specific source candidates and the search for spatial clustering 
in the data independent of any source location. The mean sensitivity in the 
northern hemisphere to an E^^'^ neutrino flux from a poi nt sourc e, = 



EJ^ ■ dN/dE\n^^^, given by the AMANDA experiment is |AI+07b 



$PS ^5g.2o-8GeVcm-2s-^ (72) 



Astronomy with a Neutrino Telescope and Abyss environmental RESearch 
NEutrino Mediterranean Observatory 

Neutrino Extended Submarine Telescope with Oceanographic Research 
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for 5 years of data taking AI"'"07b . 



The search for single point sources can be complemented by stacking classes of 
sources according to the direct correlation between the photon output and the 
potential neutrino signal. This was done with AMANDA data for 11 differen t 
AGN samples that were selected at different wavelength bands, see AI"'"06bl . 
The optimum sensitivity was typically achieved by the stacking of the ~ 10 
strongest sources in the catalog. The principle of source stacking is indicated 
schematically in Fig. 18. As an example, three source classes, FR-I and FR-II 
galaxies as well as blazars, are displayed. The stacking limits are obtained by 
stacking the most luminous sources of the same class in the sky, indicated by 
sources with filled circles. Weaker sources (empty circles) are not included. 



FR-I 

w 

Q Q FR-II 



Blazar 




filled circles: stacking included 
empty circles: stacking excluded 



© 





Fig. 18. Schematic figure of the principle of source stacking. Objects are classi- 
fied into source classes (here: FR-I and FR-II galaxies as well as blazars). Under 
the signal hypothesis that photons in a certain bandpass are directly correlated to 
neutrinos, the sources are selected according to their strength in the correspond- 
ing bandpass. The most intense sources are then stacked for an analysis of a neu- 
trino signal from the source class. This is indicated for the case of FR-I galaxies 
as circles around the objects. Only the strongest sources are used in the stacking 
analy sis (fille d circles), while weaker objects are not used (empty circles). Figure 



after B+OTa 



The search for GRBs is almost background free, since the events are selected 
according to the time and the direction of the GRB. All bursts in the northern 
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hemisphere during the operation time of AMANDA have been analyzed by 
temporal and spatial selection in the data. This leaves a vanishing amount of 
atmospheric neutrinos in the sample. The background is estimated by choosing 
a 2 hour time window around the burst duration. A similar analysis was done 
for cascade events where also the southern hemisphere can be analyzed. Ad- 
ditionally, AMANDA searches for GRBs which are too weak to be detected 
in photons but have a potentially strong neutrino signal. Here, the data are 
analyzed according to density fluctuations in the temporal structure. No sig- 
nal above the atmospheric neutrino background was detected in any of the 
searches. Limits on the neutrino flux from GRBs could be set for each of the 
searches as indicated in Fig. 19. A Waxman-Bah call spectrum was assumed 



where each burst is parameterized the same way [WB97I . |WB99| . The theory 
of GRB neutrino emission is discussed in detail in Section 7. The search for 
up-going muons gives limits on the muon neutrino flux (Fig. 19(a)), while the 
cascade analysis gives all-flavor limits (Fig. 19(b)). Analyses are now underway 
which use individually parameterized spectra for each GR B. For m ore details 



on each of the GRB analysis methods of AMANDA, see |AI+05l |. While the 
AMANDA sensitivity is not high enough to reach the sensitivity of the flux 
predictions, IceCube is at the level of the diffuse signal prediction. Also, 
for individual parameterizations of GRB spectra, the strongest bursts would 
contribute with about ~ 1 event /km^. 

In the diffuse analysis, high-energy (HE) events from all directions are exam- 
ined with respect to the spectral energy behavior of the sample. A flattening 
of the total neutrino spectrum at E^, ~ 10^ — 10^ GeV is expected when a 
flat, astrophysical component {dNy/dEy ~ E^~'^) overcomes the steep atmo- 
spheric background [dNi,/dE^ ~ Ey~'^''^). The most restrictive limit = 
Ej^ dNjdE^ is given by the AMANDA experiment, from the diffuse analysis 
for the years 2000 to 2003, 

$DL = 7.4. 10-8 GeV cm-2s~Sr~^ (73) 



in the energy range of 4.2 < \og{E^/GeY) < 6.4, see (AI"'"07d . 

Diffuse neutrino flux limits are typically obtained by optimizing the analy- 
sis cuts on spectra. The dependency of the resp onse function of the 



detector to different spectra is considered in AI"'"07d| and limits could be 
set for different spectral shapes (e.g. E^~'^) or specific models. Varying the 
spectral index in the simulation shows that the event distribution simulated 
for AMANDA peaks at very different energies depending on the assumed 
spectral index. While, for an E^'"^ spectrum, 90% of the signal lie s between 



4.2 < log(£',,/GeV) < 6.4 for the analysis presented in |AI+07dl ]. an E^~-^ 
spectrum shows an event distribution located about an order of magnitude 
lower in energy while an Ey~^ spectrum shifts the sensitivity to higher en- 
ergies. This shows that it is useful and more accurate to model the spectra 
according to the predicted shape. 
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Fig. 19. Diffuse neutrino flux prediction by WB97 . WB99l | and limits to the flux 
prediction for different channels, (a) results based on the d etection of neutrino-in- 
duced muons [AljlOTq]; (b) cascade measurements [AI^07al |. 



4-5 Detection of extremely high-energy neutrinos 



Current neutrino experiments are able to measure the atmospheric neutrino 
spectrum up to ~ 10 TeV w i thout observiri p; a significant contribution from 
extragalactic sources AI+OTdl . iM+OTl . iMiinOTl ] . Successor experiments like ICE- 
CuBE and KM3NeT aim at the detection of neutrinos at 100 TeV up to 
100 PeV. The detection of a signal at even higher energies is restrictedly 
also possible with IceCube, but to achieve a good detection probability of 
the cosmogenic neutrino fiux, new methods are being developed which are 
complementary to optical detection. Acoustic and radio neutrino detection 
aims at the measurement of neutrinos at Extremely High Energies (EHEs), 
i.e. E-oiin ~ 10^ GeV. In addition, horizontal air-showers or air-showers from 
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rocks, the sea or the Earth are neutrino signatures, since other particles are 
not able to leave traces through such dense matter. Experiments like Auger, 
HiRes, magic, etc. can investigate the neutrino flux at energies around 
~ 10* GeV by searching for (sub-) horizontal events. 



4-5.1 Radio detection 

As it was discussed for the first time by Askaryan |Ask62 |. the Cherenkov effect 
is not only present at optical and UV wavelengths, but has a strong component 
in the radio. An electromagnetic shower is caused by a neutrino, produced by 
interactions of the secondary lepton with electrons from the molecules in the 
medium. Positrons from pair production annihilate with electrons in the atoms 
which leads to a negative charge excess of approximately 30%. The shower-size 
is around ~ 10 cm and the power P of Cherenkov radiation is proportional to 
the frequency u and the bandwidth Au: 



P (X u ■ Au . 



(74) 



The refraction index for light in a medium changes back to n ~ 1 at frequencies 
above the UV light and the differential power decreases. Therefore, Cherenkov 
radiation is optically seen as blue light. The power is proportional to the 
squared electrical field. 



E 



and for optical frequencies, E 
Thus, the power is proportiona 



(75) 



with N as the number of particles. 



to N: 



Popt oc N . 



(76) 



For lower frequencies the corresponding wavelength is bigger than the size 
of the shower (A > d) and the electromagnetic fields of the particles are in 
phase. Consequently, the amount of the electric field scales with the number 
of particles, E cc N. Therefore, 



E 



(77) 



The ratio between radio power P^ad and optical power Popt of the Cherenkov 
light is then 



rad 



P, 



N 



opt 



(7J 



Typical frequencies are Urad = (100, 500) MHz and Uopt = 75 ■ 10^ MHz and 
thus, 

^(1.8- 44.4) -lO-'^iV. (79) 



opt 



A particle excess of ~ 10^*^ — 10^*^ would be necessary to reach an intensity 
equivalent to the optical intensity. However, lower efficiencies are sufficient for 
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the observation of the ef fect: t he As karyan effect w as observ ed for cascade 



hke events in sand S"'"01 , salt G"'"05 



GA+06 . While cascade 



and also in ice 

signals result in a significant radio signal as detected, single muons do not ra- 
diate enough to be observed. It is used for the purpose of neutrino detection in 
the RICELiil experiments which aims at the detection of the radio signal from 
neutrino-induced cas cades a t energies of (10^^, 10^^) eV in the ice below the ge- 
ographic South Pole Kjl03|. A geometric volume of 0.008 km^ is instrumented 
with dipole receivers at a depth of 100 m to 300 m below the surface. The 
advantage of the high-energy threshold is that the atmospheric background of 
muons and neutrinos can be neglected. The challenge is the very small signal 
rate at such high energies. So far, no neutrino signal could be identified yet. 
Most recent studies concerned the search for indivi dual GR Bs, with limits to 
ORB afterglow neutrino emission models given in B"'"07bl ]. A second detec- 
tor trying to measure the radio signal from ne utrino-induced cascades is the 



ANITAl£! 

of (1018-5 



BA^06l |. sensitive in the energy range 



balloon-borne experiment 
]^g23.5^ eV. Radio emission arising from neutrino interactions with 
nucleons in the Antarctic ice sheet is measured from above. First limits come 
from test-flights of ANITA-LITE, the first full ANITA flight was successfully 
done in December 2006, with 3.5 circles around Antarctica, with no results 
published yet. ANITA is expected to investigate the model of GZK neutri- 
nos, either confirming a positive signal, or setting limits which will restrict 
source evolution scenarios for the production of ultra high-energy protons, see 
Section 9.1. 

Next generation detectors of radio emission from neutrino-ind uced casc ades 



are planned in salt and in ice. In ice, the AURAI I experiment |GME93l ] will 



be succeeding RICE. The final aim is the instrumentation of a surface of 
about ~ 5 — 10 km^ in Antarctica. A new project using radio ante nnas on 
top of the Ross Ice Shelf at the Antarctic coast is ARIANNAEI] |Bar06 



making use of neutrino-induced radio emission in the ice shelf. Salt is a good 
target material, since it has about 2.3 times the density of water and ice, 
and therefore also provides more than 2 times more nucleons within the same 
volume. The neutrino effective interaction volume is hence increased. The 
SalSA[2H experiment is planned, where parts of a salt dome is proposed to 



be instrumented for the p urpose o f neutrino detection. First tests are already 
being done in salt domes V"'"05bl . 



Another method of detecting neutrino-induced radio emission is to observe 
the moon with radio telescopes. The moon yields a good target for neutrino 
interaction and therefore, an enhanced radio signal is expected from the di- 



Radio Ice Cherenkov Experiment 

ANtarctic Impulsive Transient Antenna Experiment 

Askaryan Underice Radio Array 

Antarctic Ross Iceshelf ANtenna Neutrino Array 

Salt dome Shower Array 
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rection of the moon. Projects investigating this effect are planning to use the 
WSRiS or SKAED. T he latte r is planned to be used in the context of the 
lunaskaE^] project |LUN08l |. For a n overview of the status of neutrino 



radio detection methods, see e.g. jSalOG 



4-5.2 Acoustic detection 

Detectable acoustic signals are expected to accompany neutrino-induced cas- 
cades. While cascades are compact and have high energy densities, which 
produces the acoustic signal, single muons have too low energy density to be 
detected as acoustic pulses. An estimate of the thermal energy deposit and 



the energy densities is given in |Lea79l | . The SAUNDl I |VGL05l | experiment 



was the first array built for acoustic neutrino detection. It makes use of a 
military basis in the Bahamas, located between depths of 1570 m and 1600m, 
spread over a surface of ~ 250 km^. First results were presented in VGL05 



The acoustic technique in the ocean is also tested at the ANTARES site 
with the AMADEUS l^^l d etector, currently being deployed along with the 



ANTARES strings Njl06|. The main aim of these prototype detectors is the 
study of the background and general water properties. For the detection of an 
EeV neutrino signal, an extremely large array (~ 10 km^ surface area) needs 
to be deployed. In addition to water, ice and salt are good environments 
occurring naturally for the installat ion of a large v olume acoustic neutrino 



detector. In the SPATScU project B"'"05l . lAI"'"06a| . first hydrophones were 



deployed on three IceCube strings as of December 2007 in order to test the 
acoustic properties of the Antarctic ice. The final aim is again the construction 
of a ~ 10 km^ array, in a depth of 200 — 1500 m, for the detection of EeV 
neutrinos. For a su mmary of the current status of acoustic neutrino detection 
arrays, see 



Van06 



4-5.3 Detection of neutrino-induced air showers 

The hadronic and electromagnetic cascades which are produced in neutrino- 
nucleon interactions via the neutral and charged current can also be observed 
as air showers, in particular in the case of tau neutrinos. Neutrino-induced 
showers can be distinguished from regular Cosmic Ray showers by choosing 
events which originate from close to or below the horizon. While Cosmic Ray 
showers are absorbed by Earth, showers from Earth-skimming neutrinos are 
induced in the outer mantle of the Earth and propagate further through the 
atmosphere. Neutrinos traveling through the core of the Earth are attenuated 
and only nearly horizontal events can be observed as inclined air showers. 



Westerbork Synthesis Radio Telescope 
^"^ Square Kilometre Array 

Lunar Ultra high-energy Neutrino Astrophysics using SKA 
Study of Acoustic Ultra high-energy Neutrino Detection 
ANTARES Modules for Acoustic DEtection Under the Sea 
South Pole Acoustic Test Setup 
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see e.g. |Zas05l ]. Searches for neutrino-ind uced air s howers are p erformed by 
the Cosmic Ray detectors, e.g. by Auger |AA+07bl | and HiRes jMH+07| . A 
Cosmic Ray Tau Neutrino Telscope (CRTNT), dedicated to the detection 
of neutrino- induced air show ers, is proposed with possible locations in Tibet 



near the ARGO-YBJ array H"'"07l |. Cherenkov telescopes with the primary 
aim of TeV ph oton show ers are also performing neutrino-shower analyses, see 
MAGIC lGM"'"07bl l. The detection o f the r adio signal from air-showers 

also gives t he pos sibility for 



e.g 



F+06 



S+07a . At these 



as planned by the L O FA R[^ experiment 
neutrino detection with a threshold energy of > 10^^ eV 
highest energies, exotic particle decays, like topological defects, can be tested. 
The energy threshold is too high for neutrino fluxes as produced in AGN or 
GRBs, since the maximum neutrino energy lies below 10^° eV. The detec- 
tion of neutrino-induced radio showers from the moon has similar a detection 
threshold and is investigated by experiments like FORTE] I. GLUEcil and 
NuMoON. 



4-6 Hybrid arrays 

The combination of different experimental techniques in one experiment or at 
the same location has been proven to help the understanding of the fundamen- 
tal processes to be measured. The latest example is the Auger experiment 
which uses the complementary pieces of information about the air shower from 
both fluorescence and surface detectors. Events can be reconstructed with 
higher accuracy and systematic errors can be reduced. A sim ilar approach 



done with the combination of IceCube and IceTop [Gai06|. Surface sta- 
tions are deployed together with the IceCube strings in order to measure the 
electromagnetic component of air showers at the surface in combination with 
the high-energy muons in the air showers in IceCube. High-energy muons 
origina te from the fragmentation region of hadronic showers, which is not well 



known |Gai06l | . The combination of a surface array with an underground muon 



detector gives a good energy resolution, while the measurement of low-e nergy 



muons at the surface provides a better mass resolution, see e.g. |Gai06 . 

In the case of neutrino detection, the three techniques optical, acoustic and ra- 
dio detection of neutrinos provide complementary information on the neutrino- 
induced showers. Therefore, the best reconstruction of single neutrino events 
will be possible, if all three detection methods are combined in a hybrid array. 
First simulations of a com bined array of acoustic, radio and optical array have 



been done for ice [BI"'"05l ]. Here, the IceCube configuration is used together 
with a further outer ring of strings for optical detection, as well as an acous- 
tic and radio array with a diameter of 10 km around IceCube. A scheme 
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Table 6 
Limits to an E,, 



-shaped, diffuse neutrino flux. The results are organized accord- 



ing to the energy range. 

of the configuration and of the effe ctive volumes of the three arrays and the 



combined array is given in V+05a 



4-7 Current diffuse neutrino flux limits 



In order to optimize the signal to noise ratio, it is useful to assume a certain 
spectral shape for the potential extraterrestrial contribution. This way, specific 
cuts can be applied in the data in order to reduce the background as much 
as possible while keeping most of the signal. Extraterrestrial neutrinos are 
typically believed to have a spectral behavior close to E,^~'^, and thus, limits are 
usually derived assuming an E,^'"^ spectrum. Table 6 lists current neutrino flux 
limits, organized according to the lower value of the energy range. All limits 
are given for a 90% confidence level, except for the RICE results, where 95% 
confidence level values are quoted. It is important to note that the different 
experiments are sensitive to different neutrino flavors. It is typically assumed 
that the flavor ratio for extraterrestrial neutrinos is [Uf, : : i/^) = (1:1:1). 
Therefore, single-flavor limits can generally be scaled to all-flavor limits by 
applying a factor of 3. However, since the flavor ratio of 1 : 1 : 1 is not exact 



and can even change signiflcantly at energies above 1 TeV |KW05| , the limits 
are quoted as given by the different experiments, see references in the table. 

Speciflc models can also be tested by changing the signal expectation from 
an E'"^ spectrum to a speciflc spectral shap e. Such hmits have bee n set by 
e.g. AMAND A in the ca se of muon neutrinos AI"'"07d| and all flavors AI+07e 



and Baikal AB^OGbj . It is also possible to d erive di fferential, model 



independent limits with methods as described in [KR+06|. These typically 
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trace the energy- dependence of the detector's sensitivity, giving bes t result s 
at the m aximum performance e nergy, see e.g. results from GLUE G"'"04al |. 
FORTE Il+oJ and NuMoON |S+07J. 



During the past decade, neutrino limits could be improved from values of 
<|)DL ^ ^'jQ . ^^-2 g-i rp^Y gng].gigg (Frejus) to four orders 

of magnitude lower values, reaching limits as low as 
<|)DL ^ Q Y4 . 10^7 g-i g^-i pgY energies (AMANDA). With 

IceCube, it will be possible to further improve this value by one order of 
magnitude within only one year of full operation, and up to two orders of 
magnitude in less than five years, giving the opportunity to detect the first 
extraterrestrial neutrinos at >PeV energies. 



66 



5 Galactic sources 



The basic assumptions for neutrino production in galactic sources are the same 
as for extragalactic neutrino sources: in the case of optically thin sources, TeV 
photon emission can be directly correlated to neutrino emission in the case of 
a hadronic acceleration scenario. In addition, the emission at lower energies 
can be interpreted as an avalanched high-energy signal. One striking difference 
between galactic and extragalactic neutrino spectra is the maximum energy, 
since galactic sources produce maximum proton energies of -Ep < 3 ■ 10^^ eV 
which is mirrored in the neutrinos' relatively low maximum energy. Candidates 
for neutrino emission are the following source classes, which are also listed in 
table 7. 

• Supernova Explosions 

Supernova Explosions release a large amount of their energy via the emission 
of MeV neutrinos. This phenomenon was observed for the first time for 
SN1987A. Today's neutrino detectors have created an early warning system 
for such e xplosio ns in our galaxy, the SuperNova Early Warning System 
(SNEWS) |A+04 . Since neutrinos are emitted prior to the optical eruption. 



the SNEWS network can trigger optical, radio. X-ray and other telescopes 
which would probably miss the explosion otherwise. On average, a supernova 
explosion is believed to happen within our galaxy about every thirty years. 
Apart from the thermal neutrino emission prior to the optical outburst, a 
TeV-neutrino flux is expected to arise from proton -nuclei interactions in 



Type-II supernova explosions as described in [WLOlJ. Protons are believed 
to be accelerated in shock-fronts up to proton energies of Ep > 20 TeV, 
resulting in a neutrino flux up to neutrino energies of E^, < 10 TeV. 

Supernova remnants 

Young supernova remn ants m ake up of the dominant part of the charged 



Cosmic Ray spectrum |GS64| |. Neutrinos are e xpected to be produced in 



GA03 



the shocks of both shell-type SNRs, see e.g. |A1H02| ] . and plerions, see 
e.g. 



The observation of TeV emission from the shell can be interpreted as a 
signal from 7r°— decays. The spectral information from the shell-typ e SNR s 



and from plerions can be used to calculate the neutrino signal, see K"'"07 
The event rates for the single sources are of the order of 1 — 10 events in a 
cubic kilometer detector with a lower energy threshold of i?min = 1 TeV, at a 
rate of atmospheric neutrinos of the same order or even higher. This implies 
that an effective background rejection is necessary, possibly by performing 
a stacking analysis in order to observe a significant signal. 
Microquasars 

Ne utrino emission from jet-disk like microquasars was suggested 
in [Bed05|. Here, it is proposed that nuclei being accelerated in the jet 
of the system interact with the accretion disk of the compact object or di- 
rectly with the accompanying Wolf-Rayet star. While photons are absorbed 
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in the dense star or accretion disk, most neutrinos traverse the objects and 
escape the production region without significant energy loss. The same kind 
of processes can happen with B- and 0-type stars instead of the Wolf-Rayet 
star. 

The detection of a periodic TeV photon signal from the two microquasars 
LS 5039 and LS I 61+303 encourages models predicting high-energy neu- 
trino emission from massive star/compact companion- type binaries. If the 
TeV photons are produced in photohadronic interactions, the ne utrino sig- 
nal is directly correl ated to the TeV photon signal as described in THOTj. A 
model presented in BotOTj ] predicts an anti-correlation between very high- 
energy photons and neutrinos. While the region is optically thin in the high- 
state of high-energy photon emission, it is optically thick in the low-state, 
allowing for the production of high-energy neutrinos. 
Soft gamma repeaters and anomalous X-ray pulsars 

Most peculiar events in our galaxy are the so-called soft gamma repeaters or, 
if the detection is in X-rays rather than in soft gamma-rays, the anomalous 
X-ray pulsars. Their emission is irregular and does not follow a periodic 
structure. In addition, the intensity of the fiares varies immensely. Low- 
intensity fiares are typically thermal, while the most intense fiares have a 
non-thermal component. The latter are therefore good candidates for neu- 
trino emission. The n eutrino e mission from the giant fiare of SGR 1806- 



20 was calculated in |HLM05l |. This fiare had a total energy release of 
~ 10^^ erg, which is several orders of magnitude above typical fiaring states. 
Th e non-th ermal photon spectrum is extrapolated to TeV photon energies 
in 



HLM05 



assuming that the TeV signal is proportional to the neutrino 
luminosity. A detection rate of a few events is expected in a cubic kilometer 
array for such a giant fiare. 
Molecular clouds and diffuse TeV emission 



Molecular clouds have been observed by H.E.S.S. Aha^06g| |. These are 
good candidates for neutrino emission: TeV emission is observed, since an 
astrophysical accelerator hides behind the clouds. The protons are absorbed 
by the dense clouds {uh ~ 250/cm^) and produce TeV photons and high- 
energy neutrinos through proton-proton and proton-photon interactions. 
Recently, MIL AGROH detected diffuse TeV photon emission from the 
Cygnus region IGM+061 . This is a go od candidate for a neutrino production 
site as well as pointed out in HOOT . 



So far, only the outer part of the galactic plane can be observed thoroughly, 
since the most sensitive instruments, AMANDA and IceCube are located 
in the southern hemisphere. Due to the detector location, the northern hemi- 
sphere is visible (see Section 4). With KMSNeT to be built in the north, 
the southern hemisphere and therefore galactic center will be observable. The 
remaining challenge is the small maximum neutrino energy, E^, < 100 TeV. 



Multiple Institution Los Alamos Gamma Ray Observatory 
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Table 7 

Examples for galactic sources of neutrino emission. Source classes for which neutrino 
emission is likely are listed. The origin of the neutrino flux is listed together with the 
electromagnetic or hadronic flux which is used to normalize the neutrino spectrum. 
In addition, the maximum neutrino energy is listed. The last column gi ves refere nces 
of proposed models. A more detailed description can be found in e.g. 



For a dedicated summary of galactic neutrino sources, see BBM05 
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6 Neutrinos from Active Galactic Nuclei 



The most luminous, persistent objects in the sky are Active Galactic Nuclei. In 
hadronic acceleration models, it can be assumed that for each class of AGN, 
the electromagnetic emission is correlated to a neutrino signal. Apart from 
individual normalization factors, the corresponding cosmological integrations 
are basically mathematically identical. In this section, the correlation between 
the emission of neutrinos and photons at different wavelengths will be reviewed 
according to neutrino flux models which are currently being discussed in the 
literature. 



6. 1 Calculation of a diffuse neutrino flux 



The total diffuse neutrino spectrum from a certain source class at Earth is 
given by the single source spectrum d^y/dE^, convoluted with the redshift 
and luminosity dependent source number. A factor l/(47r(i|) is applied to 
account for the decrease of the source flux with the luminosity distance d^. 
The number of sources is given by the product of the luminosity function of 
AGN, dn/dL/dV, and the comoving volume dV/dz. Thus, the total diffuse 
neutrino flux dNj^/dEy can be written as 

dNy,r.^ f f , ,^d^u,^n ^ ^ dn ,^ . dV 1 , , 

^{El)= dzdL—^{ElL,z)-——iL,z)- — - . (80) 



dEy^ JzJl dEy' ' ' dLdV^ ' ' dz At: dii^zf 

A shift of the neutrino energy at the source Ey to lower energies at the detec- 
tion site due to the expanding Universe, E^ = E^/ (1 + ^) is taken into account. 
A simple cosmology of VLm = 1 and I^a = is still used for many models, since 
most spatial distributions of the source samples were given according to this 
cosmology at the time of the publication of the models. The changes in the 
normalization of the res ult whe n using the experimentally conflrmed r nodel o f 



~ 0.3 and I^a ~ 0.7 |S"'"07b| | should be negligible, as pointed out in |DP90 



6.1.1 Different classes, different normalization options 

The overall normalization of the neutrino spectra builds on the correlation be- 
tween the observed flux from photons or charged Cosmic Rays as discussed in 
Section 3. In the following section, the prevailing diffuse neutrino flux models 
are discussed. Table 8 summarizes the different models with their very basic 
assumptions. Figure 20 shows a selection of the models in comparison with 
current AMANDA data and the atmospheric prediction. 



6.1.2 TeV blazars 

Sources of electromagnetic TeV emission can be interpreted as opticall y thin to 
photon-neutron interactions, T^n 1, in hadronic acceleration models |M+03bl . 
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Fig. 20. AGN neutrino flux models: ® u-VleV correlatio n for blazars Ste05l ]: @ 
z^-radio correlation for FR-II galaxies & FSQRs [BBROSal ] . usi ng an opt ical depth 
of r = 0.2; ® maximum contribution from EGRET sources jMPROll ]: ® upper 
bound for us from optically thick sources (upper, s traight b ound) and from optically 
thin sources (lower, curved bound of shaded area) MPROl| : ® Optim um pred iction 
for High-peaked BL Lacs (HBL-opt) within the pro ton-bla zar model M"'"03bl |. The 
atmospheric prediction is shown as the dotted line [Voisd ]. data are for four years 
of AMANDA-II data [m+OTI . lMun07l |. as well as the limit given in jAI+07dt . 



MPR01| . In such a scenario, charged Cosmic Rays (CRs) are produced in the 
vicinity of the source through the decay of the escaping neutrons. In this case, 
the resulting neutrino energy density would be proportional to the extragalac- 
tic CR component measured at Earth. The neutrino flux from High frequency 
peaked BL Lacs (HBLs) was cal culated using the connection between Cosmic 
Rays and neutrino emission, see M"'"03bl | . The neutrino emission peaks at very 
high energies, i.e. ~ 10* GeV. Due to the low intensity of the signal, it will 
be difficult to detect such a contribution. A theoretical upper bound of such a 
contribution to the diffuse neut rino flux was derived by Mannheim, Protheroe 
and Rachen {MPR hound), see |MPR,ni| and Fig. 20. 

As discussed in Section 3, the high-energy photon tail is absorbed by the 
extragalactic background light. Therefore, only the local Universe can be ob- 
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Table 8 



Neutrino models for Active Galactic Nuclei. The models are developed for different 
AGN sub-classes, using different signal hypotheses. Each model uses either charged 
Cosmic Rays or non-thermal photon emission from the given source class to nor- 
malize the neutrino spectrum. The CR flux gives evidence for proton acceleration, 
protons interact with the photon fleld to produce neutrinos. The correlation be- 
tween neutrinos and MeV-GeV photons can be present if the photons arise from 
vr'^-decays which implies the coincident production of charged pions, the latter de- 
caying to leave neutrinos. The normalization to the radio signal from AGN can be 
used, since the radio power is connected to the total power of the AGN via the 
jet-disk model. A fraction of the total power goes into neutrinos. X-rays originate 
from the foot of AGN jets. In the cited models, these X-rays are assumed to be 
produced in vr'^ decays, the signal cascading in an optically thick environment from 
TeV-energies to X-ray wavelengths. 



served with TeV photons. The most distant source observed by lACTs so 
far is 3C 279 at a redshift oi z = 0.538, detected with the MAGIC tele- 



scope [PMiQTl], which is considered as the VHE photon horizon here. In the 
following, it is assumed that the sources can be observed to this maximum 
redshift. Neutrinos, being co-produced with TeV photons in hadronic scenar- 
ios, offer the possibility to observe the same source classes, but without any 
losses. The spectra to be observed at Earth correspond to those emitted at 
the sources, only shifted in energy by 1 + z. The neutrino flux, originating 
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from the local sources, which can be observed in coincidence with TeV photon 
experiments, is a factor of //(^max) smaller than the total neutrino flux. The 
absorption factor 77 is determined by calculating the ratio of the neutrino flux 
up to the first formation of galaxies and the total neutrino flux up to -Zmax, 



V 



dNy/dE^{total) 
dNy/dE^{zms.^) 



The total diffuse neutrino flux cannot exceed flux limits as determined from 
AMANDA data, the most stringent one being 



= 7.4 ■ 10^^ GeV cm~2 sr' 



f82) 



as presented in AI"'"07dl ]. Therefore, the flux which can be observed simulta- 
neously in TeV photons and neutrinos is restricted to 



ah/,. 



DL 



(83) 



The absorption factor 77 is presented in |B"'"07al | with values of ri{zraa.^ = 0.5) = 
17 and rj^Zmax = 1) = 5. Assuming a horizon for TeV photon detection of 
^max = 0.5 or Zjnax = 1 rcsults in a maximum flux of 



and 



E,, 



0.5) < 4.3- 10"^ GeV cm-^g-^ 



sr 



1) < 14.8 ■ 10"^ GeV cm-2 s"^ 



sr 



B4) 



5) 



dE, ' 

respectively. This result implies that the diffuse fluxes to be observed in both 
TeV photons and neutrinos are relatively small, even for IceCube prospects. 
It emphasized the importance of neutrino analyses inde pendent of phot on 
measurements, such as the general diffuse search, see e.g. 



AI+07d. Miin07 



6.1.3 EGRET and COMPTEL hlazars 



In the case of optically thick sources (r^„ > 1), the TeV photons from 7r°— decays 
interact with nucleons in the source before escaping at lower energies leading 
to the emission of sub- TeV photons. Therefore, the diffuse extragalactic back- 
ground measured by the EGRET0 experiment |SMR05| {E^ > 100 MeV) 
can be interpreted as an avalanched TeV signal from blazars and can thus 
be used to normalize the neutrino flux from EGRET or comptelEi]- 
type source s. Again, an upper bound to the contribution from such sources 



is given in [MPROl[ |. which is much less restrictive than the optically thin 



case. Apart from the bound, a calculation of the maximum contribution from 



Energetic Gamma Ray Experiment Telescope 
^'^ COMpton TELescope 
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bl azars is g iven in MPROll ]. In addition to the contribution from HBLs given 



in |M"'"03bl |. a contribution from the optically thick Low- frequency peaked BL 



Lacs (LBLs) can be calculated using the EGRET diffuse extragala ctic back- 
ground for a normalization of the neutrino spectrum, see M"'"03bl |. A model 
of proton-photon interactions in AGN and collisions of protons from the core 
with protons of the host galaxy is derived in Man95l | . For sources with an even 
higher optical thickness, only photons below 100 MeV escape. In this case, 
the neutrino signal can be normalized to the diffuse extragalactic contr ibution 



measured by COMPTEL at energies in the range of (0.8, 30) MeV K"'"96 



This would enhance the contribution of neutrinos fro m proto n-photon inter- 



actions by almost an order of magnitude as shown in |Man95 



A model by StS96l ] was originally using the diffuse cosmic X-ray background 
and has recently been modified in a way that it is usi ng the GOMPTEL 



diffuse background to normalize the neutrino spectrum [Ste05|, see Fig. 20. 
This reduces the formerly very high contribution by a factor of 10. In addition, 
oscillations have been taken into account which leads to a further reduction 
by a factor of 2. 



6.1.4 X-ray detected AGN 

The measurement of the diffuse extragalactic contribution in X-rays by RO- 
has raised the question whether it is produced by radio-weak AGN. 
Assuming that the X-ray emission comes from the foot of the j et, the X -ray 
signal would be accompanied by a neutrino flux. A model by NMB93I ] and 



an approach by |StS96[ | have been presented, predicting a correlation between 
X-ray emitting AGN and neutrino emission from the same sources. In these 
models, p 7 interactions produce high-energy photons, which cascade down to 
X-ray energies before escaping due to the optically thick environment. An al- 
ternative scenario for the explanation of the X-ray flux from AGN would be 
the up-scattering of thermal electrons via the Inverse Compton effect. In that 
cas e, the X-r ay component would not be accompanied by a neutrino signal, see 
Until today, about 75% of the diffuse X-ray signal has been re- 



MSR95 



e.g 

solved by ROSAT [BH05|, with the help of Cha ndra and XMM-Newtop>1^ 
data, this number can be updated to 90%, see |BH05l | and references therein. 
More than 70% of the diffuse background are possibly connected to the X-ray 
emission of Active Galactic Nuclei most of which are radio weak. The search 
for a diffuse, extraterrestrial neutrino flux with AMANDA did not yield a 
significant signal and restricts these models strongly. 



Roentgen SATellite 

X-ray Multi-Mirror Mission Newton 
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6.1.5 FR-II galaxies and FSRQs 

The diffuse neutrin o ffux fro m FR-II galaxies using average parameters has 



been calculated in [BBR05a|. For the derivation of a source spectrum, it is 
assumed that the neutrino spectrum follows the proton spectrum in first order 
approximation. Effects from multi-pion production are neglected. The changes 
to the spectrum due to this effect would be negligible with respect to other 
remaining uncertainties in the calculation, such as the lack of knowledge of the 
sources' optical depth. For the normalization of the spectrum, it is assumed 
that the neutrino output of a single AGN is proportio nal to the total power 



of the jet given by the disk luminosity L^isk |BBR05al |. The disk luminosity 
is linked to the ra dio ou t put Lmd i n via t he jet-disk symbiosis approach which 



was developed in |FB95l . IFMB95 
chosen as E^~'^, with a resulting spectrum of 



Fal96|. The spectral behavior is generically 



E^^-—^ = 1.43 ■ Teff ■ 10~^ GeV cm 
dE, 



^6) 



for one neutrino plus anti- neutrino flavor, including oscillations. Here, Tg// is 
the effective optical depth to p 7 interactions. 

With the currently most restrictive flux limit, see Equ. (82), the optical depth 
has to be smaller than Tg// < 0.6. The maximum flux for Te// = 0.6 is shown 
in Fig. 20. The limit, however, only applies if the spectrum does not devi- 
ate strongly from the E^~^— shape. It was pointed out by different authors. 



e.g. |Bar04l . lMBQ07bl |. that the spectral shape can vary significantly under 
certain conditions of relativistic shocks. While flatter spectra may still be ex- 
cluded by the limit, steeper spectra are more difficult to investigate. When 
interpreting the radio spectrum as synchrotron radiation from electrons, the 
radio spectral index s ~ 0.8 is correlated to the primary index as p = 2 ■ s -|- 1. 
The neutrino spectrum, following the proton spectrum, could therefore be as 
steep as E^'"^'^. 



6.2 Individual source spectra 



For the calculation of the diffuse spectra discussed above, mean parameters 
have been used. In this section, different predictions implementing the indi- 
vidual variations of the single sources are discussed. 



6. 2.1 TeV emitters 

A simple estimate of the neutrino flux from TeV sources is to use the photon 
flux at TeV energies under the assumption that it arises from photo-hadronic 
interactions. It is known that the TeV photon flux is diminished by interactions 
with the extragalactic background light. Thus, the TeV photon flux at Earth 
and its spectral index are used here as a lower limit for the photon emission 
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from the blazars. The neutrino flux is given as 



dE~^'^'^^ 



with the spectral index a'^, dTeV- The normalization factor is calculated by 
correlating the neutrino flux luminosity with the TeV photon flux as 

[^KdK = K. [ ^ E^^y dEr.y . (88) 

Here , K = 1 for interactions and K = 1/4 for p'j interactions, see 



e.g. |Hal06[ | . Assuming that the spectral indices for neutrino and TeV pho- 
ton emission and the energy threshold (£"111111 ~ 100 GeV) are about equal 
translates the above equation into a direct correlation between the flux nor- 
malizations, 

A'^^K-Frev (89) 
Table 9 lists the normalization for the neutrino spectrum of the 13 HBLs with 
given spectral information, the LBL BL Lacertae and M 87. The values have 
been calculated using the parameters in table 10 in Section 3. The results 
are listed for a 1 : 1 correlation between photons and {u^ + z7^)-neutrinos as 
it is the case for proton-proton interactions. Dividing the values by 4 yields 
the flux in case of p7 interactions. It needs to be considered that the TeV 
photon flux from blazars is highly variable. The intensity of the low state can 
be several orders of magnitude below the high-state flux. Both difficulties and 
opportunities arise from the strong variability. 

The estimated neutrino signal from hig h-state p hoton states in AGN and point 



source neutrino flux limits as given in AI"'"07bl | for 6 of the individual sources 



are listed in table 9. However, the limits cannot directly be compared to the 
estimated flux: the calculated values are only valid during the time of the flare, 
while the limits apply to the flux averaged over the five years of observation. 
The latter is significantly lower and for many of the TeV sources, the low- 
emission state is not known. Low-state observations require long observation 
times. The duty cycle of Cherenkov telescopes is quite low, since only dark 
nights can be used for measurements. Therefore, extensive long-term monitor- 
ing of the sources is not possible with the current telescopes. Many flares are 
missed and the average flux is difficult to determine. Small telescopes, dedi- 
cated to the long-term monitoring of Te V blazars, are planned to be rebuilt 



from HEGRA telescopes as described in [BMlLOTbJ. Currently, the Whipple 



telescope uses part of its tim e for bla zar monitoring in order to get a better es- 



timate of the TeV lightcurve SV"'"07l ] . These attempts may be an option to get 
a better understanding of the variability of the sources. The small telescopes 
are, however, typically not sensitive to low-state emission. This has to be ex- 
amined by the more sensitive telescopes, H.E.S.S., MAGIC, VERITAS & 

CO. 
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Table 9 



Neutrino flux (i^^ + u^) from AGN of TeV emission. Limits are taken from [AI+07b| . 
having been calculated from 5 years of AMANDA data. 

The flaring states also provide an advantage for high-energy neutrino tele- 
scopes: for known flares like the ones listed in table 9, temporal cuts in the 
data of neutrino telescopes can be used in addition to spatial cuts. If only the 
duration of the flare is used for an analysis, the estimated flux values can be 
used as reference values. This way, a much higher flux is expected, in addi- 
tion to a reduced background, which increases the probability of detecting a 
signal. Analysis methods concerning AGN flare searches are being developed 
as presented in jRI+06 . 



6.2.2 FR-II galaxies and FSRQs 

A sub-sam ple of 114 FR-II sources from the steep spectrum sample given 



m 



W"'"01| was used to analyze the differences between an average prediction 



and a source property based estimate, see jBBR05bl |. The radio flux S is 
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converted into the radio luminosity as 



-'radio 



indf-S (90) 



and the spectral radio index s is used to determine the proton spectral index 
for each case individually, ap = 2 ■ s + 1 . The single source normalization versus 
the primary spectral index is shown in Fig. 21. The normalization is given at 
Ey = 1 GeV. There is a strong correlation showing that flat neutrino spectra 
are stronger than steep ones. The coincidence spectrum is much flatter than 
the average prediction, 

Using the mean spectral index s = 0.8 results in a much steeper mean spec- 
trum due to the correlation between the normalization Ay and the spectral 
index ap. 
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Fig. 21. Correlation between the neutrino flux normalization and the parti- 
cle spectral index Op. While the majority of sources is concentrated around 
Ay ~ 10~^^ GeV/s/cm^ and ap ~ 2.6, a strong correlation between Ay and Up 
is visible - bright sources are typically flatter. 



6.3 Theoretical upper bounds on neutrino emission 



An upper bound on the diffuse neutrino flux from extragalactic sources was 
derived in [WB99|. Here, the observed spectrum of UHECRs was used in order 
to calculate the maximum neutrino flux possible to be produced together with 
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the highest energy Cosmic Rays. The bound is given as 



dE„ 



1.5 ■ 10"*^GeV cm-^s^^sr"^ 



(92) 



for muon and anti-muon neutrinos. It is pointed out in [MPROl[ | that this 
bound is, however, only vahd for sources optically thin to neutron-photon 
interactions. In those sources, neutrons escape, decay and produce the flux 
of UHECRs. High-energy photons from p roton-p hoton interactions escape at 
TeV-PeV energies. The bound derived in WB99| does therefore not apply to 
those classes of AGN which are opt ically th ick to neutron-photon interactions. 
In addition, it is pointed out in MPROl| that the calculations in WB99 
assume a strict spectrum. The neutrino spectrum, however, deviates 

from the primary spectrum of charged particles, which needs to be taken into 
account in the case of the derivation of an upper bound. Taking these facts 
into accoun t, a boun d for transparent sources and optically thick sources was 



derived in [MPROll ] as indicated in Fig. 20. The upper, straight line of the 
shaded area represents the bound for optically thick sources, while the lower 
line shows the upper bound concerning optically thin sources. 

It is generally very difficult to derive absolute upper bounds. There are still 
large uncertainties in many of the parameters. Taking for instance the flux 
of high energy Cosmic Rays, which is uncertain by at least a factor of two, 
when comparing the results from HiRes and AuGER, see Fig. 2. Further 
uncertainties come from the spectral index of the primary flux and the neutrino 
flux, the redshift distribution of the different source classes, as well as the factor 
which connects UHECRs or photon emission with neutrinos. 



It is shown in Hal06 



WB99 should rather be 



that the derived bound from 
interpreted as a flux estimate considering the connection between UHECRs 
and neutrinos, with an intensity between 



2 dN^, 
dK 



;i - 5) ■ 10"'' GeV cm-2 s"^ sr"^ 



(93) 



for muon and anti-muon neutrinos, using recent Cosmic Ray data. 



6.4 North and south view - next generation neutrino detectors 

Both for point source and diffuse searches, the orientation of the fleld of view 
of the detector is crucial. The sources are not distributed isotropically. Sky 
surveys can be biased by the location of the instruments and the brightest 
source in a sample sometimes makes up the dominant contribution from a 
source class. The latter is, for instance, the case for FR-I galaxies, where M 87 
gives most of the radio flux in the total sample. The reason is that M 87 is by 
far the closest of the sources, at a distance of only ~ 16 Mpc. 
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The measured diffuse ffux depends on the location of the detector. If the 
luminosity distribution is steep, implying that nearby sources make up most 
of the contribution, the signal is much stronger with the nearby sources in the 
field of view. 

IceCube is currently being built at the geographic South Pole with a full 
view on the northern hemisphere in the muon-neutrino channel. The com- 
plementary instrument, Km3NeT, planned to be built in the Mediterranean, 
will cover the southern hemisphere. Depending on the source class, optimal 
observation conditions are provided by either southern or northern telescopes. 
For instance, in the case of EGRET-detected sources, the three brightest ob- 
jects are in the southern sky. The largest fraction of identified FR-I and FR-II 
galaxies are located in the northern hemisphere, since the supergalactic plane 
is mostly located in the northern hemisphere and there may also be a bias in 
the radio selection of objec ts. A detailed description of the most important 



catalogs is given in B"'"07a . 



6.5 Physics conclusions from current neutrino flux limits 

The primary goal of large volume neutrino arrays is the detection of a signal 
from extraterrestrial sources. Different methods have been developed in order 
to reduce the atmospheric background, using the model predictions described 
above for an optimal selection of data. Here, energy cuts are applied for diffuse 
searches and spatial and temporal cuts are used for point source searches. 
The stacking of the most intense objects of a source class can further improve 
the signal to background ratio. This was done for 11 AGN classes with the 



AMANDA detector, as described in |AI+06bl . lAI+07bl |. Here, typically the ten 



brightest sources in the sample were used for the analysis. The point source 
limits $ can be interpreted as diffuse limits, as described in |B+07a 



2 7rsr- ^ ' 

The factor e includes identified sources contributing to the diffuse flux, which 



have not been considered in the stacking analysis I. The second factor S, 
accounts for the diffuse background of unidentified sources. This is only known 
for some of the source classes and needs to be estimated in other cases. The 
factor 1/(2 vr sr) accounts for the field of view of the AMANDA detector. 

The best results are achieved for identified EGRET blazars. The stacking 
diffuse limit is given as 

•^loREr = 7.25 ■ 10-^ GeV cm'^ s'^ sr'^ (95) 



In many samples, most of the sources are too faint to be considered in the stacking 
analysis - including them would reduce the sensitivity. They still have a significant 
contribution to the diffuse flux. 
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in the energy range of 1 TeV< E^, < 1000 TeV. This is not a general diffuse 
hmit, but it only applies for diffuse flux predictions building on a correlation 
between the photon flux as measured by EGRET and neutrino emission. 
While the limit cannot restrict current predictions, second generation neutrino 
telescopes like IceCube and KM3NeT will be sensiti ve to the prediction 
of neutrinos from pp interactions described in Man95l ] at energies around 



1 TeV. The advantage of the stacking method is that it is sensitive to neutrino 
fluxes at relatively low energies, while the general diffuse search is restricted 
to energies of above 10 TeV due to the high atmospheric background. 

The diffuse search, on the other hand, is more sensitive in most of the cases 
at higher energies, since the stacking method loses sensitivity due to the high 
diffusive factors ^ which can reach values up to ~ 10^ in the worst case. With 
current neutrino flux limits , it is po ssible to exclude the correla tion between X- 



rays an d neutrino emission B^07al ]. The three predictions, i.e. StS96l . lNMB93 



AM04j ] show a violation of the limit by about an order of magnitude. It is 
therefore excluded that neutrinos are produced in coincid ence w ith X-rays 



from radio-quite AGN. The modified model of Stecker et al. [Ste05|], assuming 



a correlation between neutrinos and MeV photons rather than X-rays is a 
factor of 20 lower than the original prediction. The factor of 20 arises from 
the fact that the diffuse MeV flux is lower than the X-ray flux, which translates 
directly into a lower neutrino flux. 



A search for ultra high-energy neutrinos was performed in AI"'"07el |. At en- 
ergies around 1 PeV, neutrinos are effectively absorbed by the Earth 
and only reach the detector when coming from close to the horizon. Events 
from above are difficult to distinguish from atmospheric muons, which leaves 
horizontal neutrinos for the analysis. The limit for three year of data taking 
(year 2000-2003) lies at 



= 2.7 ■ 10"^ GeV cm-2 s"^ sr"^ . 



(96) 



Limits have been derived acco rding to different models. The limit lies a factor 



of 1.25 below the prediction of |MPROll |. presenting the maximum contribution 



from optically thin sources. The most optimistic parametrization can there- 
fore be excluded. With IceGube, a sensitivity gain of about two orders of 
magnitude in one year of observation is expected, providing more information 
on the parameters of different flux models at high energies. 
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6. 6 The Olbers ' paradox for neutrinos 



First predictions of high-energy neutrino emission were made in the 1990s 
with maximum parameter settings in order to estimate what is possible. These 
models could now be revised using more realistic parameters due to current 
neutrino flux limits. While it is obviously the aim to detect high-energy neu- 
trinos from extraterrestrial sources, the limits which have so far resulted from 
the searches are very valuable in the context of multimessenger astrophysics as 
discussed above. These limits are in analogy to the Olbers ' paradox. The para- 
dox arose when Olbers calculated that - with a homogeneous distribution of 
astrophysical objects in the sky, and with these objects radiati ng isotr opically 
- that the sky should be isotropically bright as the sun itself 01b26l |. There 
was, however, a very obvious, optical limit on the radiation during night-time. 
The sky was, and still is, dark! The paradox could not be solved before the 
1950s, when different assumptions from Olbers' calculations fell: the sources 
are not isotropically distributed in the sky, the Universe is not inflnite and it 
is not in thermal equilibrium. 

Similar considerations can be done for every other wavelength, where a pre- 
diction is made which violates experimental limits. The limits give the oppor- 
tunity to constrain the parameters of the model. With neutrino-astrophysics, 
we are still at the beginning of exploring these predictions. Still, current neu- 
trino flux limits start to become sensitive to prevailing neutrino flux models 
and start to constrain them. Within a short time, the next generation of large 
volume detectors will be able to reach sensitivities more than 1 — 2 orders of 
magnitude better than current instruments and it can be expected, that the 
flrst extraterrestrial neutrinos will be detected soon. 
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7 Gamma Ray Bursts as neutrino sources 

A summary of the photon and neutrino emission from GRBs is shown in 
Fig. 22. Three phases o f non-the rmal emission are expected, i.e. pr ecursors 
hours prior to the GRB IRMW03I I . em ission d uring the prompt phase WB97 



jprior to tne LrKo [KM W UdlJ . em ission d uring tne prompt pnase |VVoyi, 
9| as well as afterglow emission |WB00| . The lower row shows the neu- 



trino energy spectra for the three phases while the upper row displays the 
corresponding electromagnetic output schematically. 

The basic idea is that 



A precursor model has been developed in [RMWOS 



a shock forms when the pre-GRB matter collides with the wind of the 
central pulsar or the SNR shell. At this point, the burst is still opaque 
to photon emission and the shock environment yields a good target for 
neutrino production by shock-accelerated protons interacting with thermal 
X-rays in the sub-stellar jet cavity. The shocks happen at smaller radii than 
the prompt emission and at lower boost factors F. The neutrino signal could 
also be accompanied by a signal in the far infrared. The low energy part of 
the neutrino spectrum arises from neutrino production in pp interactions 
and is -E^~^-shaped. The flux at energies > 10^ GeV originates from 
proton-photon interactions. The two lines in the lower left panel of Fig. 22 
represent different shock/jet radii and envelope mass settings. 
The prompt photon emission from GRBs can be correlated to the observed 
flux of UHECRs, since protons are belie ved to be accelerated in the highly 
relativistic shocks Vie95l . IWB97l . IWB99l ] . The acceleration of protons in turn 
implies the production of neutrinos through photon-hadronic interactions. 
The resulting flux is shown in the lower middle graph of Fig. 22. 
Afterglow neutrinos are produced when the internal shocks from the original 
fireball hit the interstellar medium and produce external 
shocks 



WBOO 



As for the prompt emission, the synchrotron emission of 
electrons gives evidence for the existence of relativistic charged particles 
which in turn implies neutrino production by the baryonic component of 
the jet and t he pho ton field. In the case of afterglow emission, Waxman 
and Bahcall |WBOO | predict the acceleration of ultr a high- energy protons 



{Ep > 10^° eV) in reverse, mildly relativistic shocks MR97l |. The accelera- 



tion of protons to such high energies implies the production of neutrinos in 
proton-photon interactions in environments optically thick to proton -photon 



or proton-proton interactions. Waxman and Bahcall conclude in |WBOO 



that a significant neutrino flux during the afterglow phase is directly corre- 
lated to the electromagnetic afterglow emission. 

The three different emission phases are discussed in more detail in the next 
paragraphs. All variables are given in the observer's frame. 

A general question is the timing of photon and neutrino emission. Depending 
on the environment of the astrophysical shock, neutrinos are not exactly emit- 
ted at the same time as the photons. This is important in the case of temporal 
selection of data in high-energy neutrino experiments. It is typically still as- 
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7.1 Precursor neutrino emission 



Since the observation of GRB980425 and GRB030329, it is known that at 
least a fract ion of long GRBs originates in heavy supernova explosions, see 



e.g. [MiQSaj]. The emission is beamed, so that not every supernova explosion 
can be observed as a GRB as well. The Gamma Ray Burst can only be detected 
when the jet is directed towards Earth. The pulsar, remaining in the center of 
the SNR, is believed to form prior to the GRB. The pre-GRB is propagating 
outwards, and a shock is formed when it meets the magneto-hydrodynamic 
pulsar wind. In this shock, protons can be accelerated. The fast protons in 
turn interact 

(1) with the supernova remnant's photon field, 

(2) with the photon field of the MHD wind itself, and 

(3) with cold protons from the slow plasma. 



In the these interactions, neutrinos are produced |RMW03l | . The total spec- 
trum is calculated by assuming that a fraction of the fast proton spectra, 
dNp/dEp, are converted into neutrinos. 



dN, III fpp ■ M,{Ep)^ for Ep < Ep 



th 

(97) 



dE^ Andl ] UdN^ > ^ *^ 

' i dEp '■^^ ^ -^p ■ 

Here, di is the luminosity distance and Ep^^ 10^ GeV is the threshold energy 
of the protons for the production of the A-resonance. Below this threshold, 
pp interactions dominate, where fpp is the proton-proton efficiency. Above the 
threshold, the total spectrum is given by neutrinos generated in p7 interac- 
tions. For the case of photo-hadronic interactions, the proton spectrum needs 
to be weighted by the fraction of protons converted into pions, «i 1, and 
it needs to be considered that the flux distributes equally to the four decay 
products (e"*", Ueyi'^^T'fj). For the case of proton-proton interactions, it needs 
to be considered that the neutrino multiplicity for an incident proton, M^, is 
energy dependent as 



M,{Ep) = - N^/K ■ {{E/GeV) ■ [l/2 ■ In (lO^^ GeV/Ep)] } 



/ GeV 



(9J 



Here, G is a step-function, which is O = 1 for the given range and 6 = 
elsewhere. The product m^//^ ■7cm is the lower energy threshold for pion/kaon 
production, with m^//^ as the mass of the particle and as the Lorentz 
boost factor at the center of mass in the lab-frame. The normalization factors 
A^^ or Nx are given by integrating the distribution of pions, resp. kaons. The 
total numbers are given as N„ ^ 7 and Nk ~ 0.6. The proton-proton efficiency 
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is assumed to be fpp = 1, see |RMW03 



The proton spectra dNp/dEp are calculated in [RMW03|, assuming an Ep~ — 
shaped spectrum. The resulting flux is shown in the lower left graph of Fig. 22. 
The precursor is not observable in photons, since the environment is optically 
thick. Thus, the detection of neutrinos from precursor GRBs in the future 
provides th e possi bility to alert photon experiments before the actual GRB 



occurs. In [KMOTj, it is proposed to use neutrino events in IceCube as a 
possible trigger for optical telescopes. If the neutrinos arrive in multiples within 
a short time window, they might originate from a core-collapse supernova. 
Such a trigger would allow for the observation of SNe in the early phase of the 
collapse. 



7.2 Prompt GRB spectra & neutrinos 



Th e prorn pt GRB photon spectrum, dN^/dE.y{E^), is usually given by a Band 
fit IB+93L 



dK 



A 



A 



100 keV 



100 keV 



^ ^{-E-,/Eo) 



100 keV 



E'y < {a-y — (3y ) Eq 
,{-K-/3^)) > -/3y)Eo 



(99) 



Here, -E^ is the photon energy, Eq is the reference energy with (a^ — ■ Eq = 
as the break energy in the photon spectrum. The function was designed 
empirically in order to match the GRB spectra. It is the best mathematical 
description for broadband GRB spectra. The energy range of detection may 
not cover both parts of the spectrum and spectral fits with simpler powerlaw 
approximations are often applied. The break energy for the photon spectrum 



IS given as e 



j,MeV 



(3y) ■ Eo/MeV. 



The spectrum is presumably produced by sy nchrotro n radiation of electrons 
in the internal shock fronts of the jet, see e.g. HaH02[ | and references therein. 
There are two approaches to explain the break in the spectrum at a break 
energy of typically e!^ ~ 250 keV. The most common explanation is the steep- 
ening of the spectrum b y one power d ue to the cooling of electrons at high 
energies, see for example ZM04I . IPirOq as a review. The break can also be ex- 



plained by assuming an Inverse Compton scattering scenario, see e.g. |DR04 
and references therein. Throughout this review, all energies concerning GRB 
spectra are given in the observer's frame at Earth unless declared otherwise. 
For a regular GRB, the spectral indices are usually distributed around aver- 



age values of 



-1 and Py 



-2. These values scatter over a wide range. 



Short GRBs tend to have harder spectra with ~ and ~ — 1. Assuming 
hadronic acceleration in the jet, a prompt neutrino flux that is correlated to 
the photon spectrum results from photo-hadronic interactions in the source. 

The neutrino spectrum, dN^^/dE^, can be derived assuming that the proton 
spectrum follows the electron spectrum of the source. Furthermore, it needs 
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to be assumed that electron losses can be neglected. The neutrino flux in 
turn follows the proton spectrum in a flrst order approximation, since it can 
be connected to the observed synchrotron spectrum of the sources. An im- 
prove ment of the calculations is generally possible by calculating the electron 



losses |You06| . Here, this effect is neglected and the result is taken as a lower 



limit for the neutrino flux. 

Protons accelerated in astrophysical shocks can interact with the source's pho- 
ton fleld and produce a A— resonance which in turn leads to neutrino produc- 
tion. Proton and neutrino energy are directly proportional, E^, = Ep/20. The 
product of proton and photon energy remains constant, since the energy of 
the A mass has to be produced in the center of mass system, Ep- = const. 
Hence, photon and neutrino energy are inversely proportional, E^, oc E^~^ . 
The neutrino flux is then given as 

(100) 

[ {EjelY^^ for el<E,<el. 

The photon spectral indices can be used to describe = [3^ + 1 and (3y = 
+ 1. It should be noted that the usage of a broken powerlaw instead of the 
original Band function can lead to increased uncertainties in the calculation as 



pointed out in [TD07|. The use of the original Band function is more accurate, 
since it is able to flt the GRB spectrum very precisely over the whole energy 
band, (20, 2000) keV. A complete Band flt, however, cannot be extracted 
for many bursts. In the case of BATSE, it was still possible to reconstruct 
the entire Band function, while with Swift, only single powerlaw flts are 
possible in most of the cases. Also, in combination with the uncertainties in 
the measurements, the difference between using a broken powerlaw compared 
to a Band function may be negligible in many cases. Therefore, for simplicity, 
a broken powerlaw is used in the following paragraphs for cases where the 
entire bandwidth is detected. 

A second break at E,y = appears in the neutrino spectrum which is caused by 
synchrotron radiation of the neutrino-producing pions. Pions with sufficiently 
high energies suffer synchrotron losses before decaying and do not produce 
neutrinos. Thus, the high-energy tail of the neutrino spectrum is steepened by 
one power, 

^ El oc (E^/el)-^^ (KK)-' for E^ > et . (101) 

The spectrum is normalized to the 7-ray fluence which is assumed to be 
proportional to the neutrino luminosity, 

x-F^= -p^dE,- In(lO) ■ A, . (102) 

All parameters occurring in the following calculations are listed in table 10. 
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The factor x is given by the product of the energy transferred to the pions, 
a factor 1/8 since half of the photo-hadronic interactions result in four 
neutrinos and a factor 1/fe to account for th e fraction of total energy in 



electrons compared to protons in the jet G^04bl |. The normalization constant 
Ai^ is therefore given as 

In the following, the normalization of a single burst will be modified to a 
quasi-diffuse normalization by multiplying A^, with the number of bursts per 
year (1000 bursts per year, of which 77 = 2/3 are long and 77 = 1/3 are short) 
and dividing the result by 4 tt sr, 

K.,.^A.^ (104) 

yr 4 7r sr 

This number is comparable to the normalization which is achieved when using 
the flux of UHECRs to normalize the spectrum, since the power in UHECRs is 
of the same order as the electromagnetic emission of GRBs. The total number 
of 1000 bursts per year is chosen t o be c onsistent with calculations given 



by Waxman and Bahcall in |WB97l . |WB99| ] . A estimate of ~ 700 bursts per 



year, however, is a more reasonable number when looking at the BATSE 
catalog |P+99l ]. In that case, all results could simply be weighted by another 
factor of 2/3. 

The first break energy in the spectrum, e'^, is related to the break energy in 
the photon spectrum by 

, (ml - ml) ■ , ^ 

It is determined through the minimal energy necessary to produce a A-reson ance 



in the shock fronts of the bursts. The numerical values given in YP^06| for 
the proton mass, = 0.94 GeV, and the A mass, m^ = 1.23 GeV, lead to 

el = 7 ■ 10^ ■ (1 + z)-^ GeV (106) 

^■y,MeV 

in the observer's frame. 

The second break energy is connected to the pion's synchrotron loss time. It 
depends on the neutrino flavor and for muon neutrinos, it is given as 



= C ty 4 

^ ^Ar^areBL, {1 + z) ■ m, ' ^ ' 

with the Thomson cross section ax = 0.665 ■ 10~^^ cm^ and r° = 2.6 ■ 10"^ s 
as the pion life time at rest. For electron and anti-muon neutrinos, the break 
energy is about an order of magnitude lower, since these neutrinos originate 
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from muon decays. The muon lifetime is about a factor of 100 higher than the 
pion hfetime, which gives a lower energy threshold for synchrotron losses. Here, 
and are the fractions of the burst's internal energy going into the magnetic 
field, respectively into electrons. The equipartition fractions have been set to 
ee = 0.1 and = 0.1. There is no good way of determining the equipartition 
fractions theoretically, yet. However, afterglow observations indicate values on 



the order of 0.1 for prompt emission |WKF98| . The remaining parameters in 



Equ. (107) are listed in table 10, with the values as used in the following 
calculations. Inserting all numerical values gives 



10^ 

l + z 



(108) 



A detaile d derivation of the neutrino spectrum as presented above is given in 
WKF98I. 



7.2.1 The average spectrum 

Th e avera ge spectrum of prompt neutrinos from regular GRBs was calculated 
in WB99| . hereafter referred to as WB. The following parameters were used. 



A. 



3 • 10^^ GeV cm-^s-^r- 



= -1 , 
/5. = 0, 
el = 10^ GeV 
ef, = 10^ GeV 



(109) 
(110) 
(111) 
(112) 
(113) 



Though the sensitivity of current neutrino detectors is not sufficient yet to 
detect this flux (see also Section 4), they come close to the predicted flux. 
With IceCube, the model can therefore be tested in detail. 

A numerica l calculation of the prompt neutrino spectrum from GRBs is given 
in MN06al | . Multi-pion production is included in the calculations by using the 
GEANT4 generator, and proton cooling was calculated. Thi s more p recise ap- 
proach results in a spectrum similar to the one calculated in WB99l | , but with 
smooth transitions at the break energies and a slightly different normalization. 



7.2.2 The parameters 

The parameters in the neutrino flux calculations have been adjusted accord- 
ing to the different classes of GRBs. A summary of the mean values of the 
used parameters is given in table 10, for Gamma Ray Bursts, X-Ray Flashes 
(XRFs), Short Hard Bursts (SHBs) and Short High- Energy (SHE) bursts. The 
different source classes are discussed in section 7.2.3. 

The boost factor F is constrained to 100 < F < 1000 for regular GRBs, 
since for boost factors less than 100, the medium would be optically thick to 
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photons a nd for T > 1000, protons lose most of their energy to synchrotron 



radiation |HaH02l ]. The possibihty of fluct uating; T using the photon break 



energy is given as demonstrated in G"'"04bl ]. but there are several arguments 



for using a constant value: bursts can be misaligned which would lead to 
a misinterpretation of the boost factor. Also, varying the break energy for 
each single burst might implicitly already include boost factor fluctuations. 
Therefore, a constant boost factor of F = 300 is used in following calculations. 



In previous publications, e.g. G^04bl ]. a variation of the energy going into 



pions was discussed. Such a variation would further increase the width of 
the distribution of the neutrino spectrum normalization. The burst luminosity, 
the boost factor, the photon break energy and the variability time influence 
as they are correlated as 

A - 0.2 . ^ . (114) 

2.5 t^^MeV 

In the following calculations, this fraction will be kept at a constant value of 
~ 0.2 for classical GRBs for the following reasons: 

(1) strongly depends on the boost factor F which will be used as a constant 
as discussed before. The dependence on the other three parameters is only 
linear. The variation of F is more striking, since it is amplified by the F'^ 
behavior. 

(2) The main uncertainties in the current calculations result from the lack of 
knowledge of the parameters Cg and as well as from uncertainties in the 
redshift relation which leads to uncertainties in L^. These three parame- 
ters are all important in the determination of the spectral normalization 
and thus a constant value is favorable. 

In the case of SHBs, the variability time is smaller, t^,_2 ~ 0.1 and also the 
luminosity is about an order of magnitude lower, ~ 0.1. Thus, the pion 
energy constant stays the same and is assumed to be /^r = 0.2 for SHBs as 
well. 

X-Ray Flashes (XRFs) will produce very high pion efficiencies due to their 
low boost factors and break energies, /tt ^ 1 indicates an optically thick 
source in comparison to a classical GRB and yields a good target for neutrino 
production. 

The pred iction of the di ffuse neutrino flux derived from the model above is 



given in |WB97l . [WB99|, where the authors use average parameters to de- 
termine the shape of the spectrum. Cosmological evolution of the sources 
was considered in that model by taking into account the redshift evolution 
of GRBs. It is assumed that GRBs follow the star formation rate, since they 
appear to be connected to supernova-Ic explosions. Thus, the diffuse Waxman- 
Bahcall flux has to be weighted by a factor ?7 = 2/3or?7 = l/3 for long and 
short burst samples respectively. 
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1.2.3 Conventional versus peculiar bursts 

As discussed in Section 3, GRBs can be divided into different sub-classes, such 
as long and short events and a further division of long bursts into XRR bursts, 
XRFs and regular GRBs. While short and long bursts differ in duration and 
in the hardness of the spectrum, XRR/XRFs and regular GRBs differ in the 
energy range of prompt emission. A further class of short bursts is predicted 
with a high-energy component, which has not be observed yet, but which 
would be observable with GLAST. 

The photon spectral properties are reflected in the neutrino spectra. A schematic 
representation of the neutrino spectra of the four different classes is shown in 
Fig. 23. Due to a low boost factor of F ~ 30, XRR and XRFs have a lower 
break energy in the photon spectrum, which leads to higher neutrino break 
energy, ~ 10^ GeV for XRR/XRFs (b) compared to 10^ GeV for regular 
GRBs (a). The second break is at approximately the same energy, 10^ GeV 
due to the combination of a lower boost factor and a higher variability time, 
tv,-2 ~ 10"^ due to large emission regions. XRRs and XRFs can actually be 
very efficient neutrino emitters, given that for the low boost factor, the pion 
parameter can become as large as /^r = 10'^. Such a high value of would, 
however, violate recent AMANDA limits, see Section 3. Therefore, the flux 
is normalized to the current limit in Fig. 23, which implies that the pion pa- 
rameter must be < 4. Other possibilities to get around the limit would be 
to use smaller scales for the emission region and therefore smaller variabil- 
ity time scales. This reduces the synchrotron break energy significantly. Also, 
there could be much less XRR and XRF events than regular GRBs, which 
would affect the normalization. The variation of the total luminosity plays an 
additional role, affecting both the normalization and the second break energy. 

The photon spectra of short bursts are typically harder than long ones which 
leads to a much softer neutrino spectrum of ~ E'"^ below the break energy and 
~ above the break energy (see Fig. 23 (c)). The class of high-energy short 
bursts would produce a neutrino spectrum as displayed in Fig. 23 (d). There 
is, however, no experimental evidence yet that such a class exists. While the 
photon spectrum would be observed at > 10 MeV energies, the corresponding 
neutrino spectrum is supposedly visible at lower energies as indicated in the 
figure. 

The reflections above show that the assumption of a (2:1) ratio of long to 
short bursts and a (1:1:1) ratio of (GRB:XRR:XRF) would result in a modified 
spectrum which can look significantly different from the spectrum which has 
so far been classified as the typical neutrino GRB spectrum. 

In addition to the classification scheme as presented above, there are two 
types of bursts which are of low photon luminosity, but can possibly be 
sources of a significant flux of high-energy neutrinos: low-luminosity GRBs like 
GRB060218 a nd chok ed bursts. The neutrino flux from low-luminosity GRBs 
is estimated in |M+06| to be at a level oidN^/dE^ EJ^ ~ 10"^ GeV cm"^ s^^ sr^^ 
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at 10^ GeV, which would be observable by IceCube after a fe w years o f obser- 



vation. Neutrino emission from choked GRBs is predicted in [MWOll ]: While 



the electromagnetic emission is not observed due to the optical thickness of 
the event, e.g. a supernova explosion into a dense SN-wind, neutrinos from 
proton proton or proton photon interactions can still escape. Further acceler- 
ation of secondary charged particles like muons, pions and kaons can lead to 
an enhanced neutrino signal KWOTj. Proton proton interaction for instance 
would only produce high-energy neutrinos up to ~ 0.1 — 1 TeV. Assuming 
further acceleration, the maximum neutrino energy lies above 10^ TeV and a 
flux of dNy/dEy ~ 5 ■ 10~® GeV cm"^ s^^ sr~^ is expected for an type 
neutrino spectrum. 

Apart from spectral fluctuations between different burst classes, there are 
relatively large variations within the classes themselves. 
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Fig. 23. Prompt emission neutrino spectra from (a) a classical, long GRB as opposed 
to an XRR/XRF event (b) and from (c) short, high-energy bursts which are not 
observed yet compared to short, regular bursts as observed by KONUS, Swift, etc 
(d). 



7.2.4 Average versus single source variation 

A different approach to predicting the neutrino flux from GRBs is to look at 
each burst individually and add up the individual spectra to make a predic- 
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tion of the t otal flux from t hese so urces. The study of different catalogs are 



presented in G"'"04bl . IB^OGI . iBecOTj ] and are summarized here. 



For each burst i in a sample of n GRBs, a prediction of the prompt neutrino 
flux, $j, from this source can be made as described above. With a large volume 
neutrino telescope like AMANDA, a large number of bursts can be analyzed 
simultaneously in order to increase the detection signiflcance. By stacking 
many bursts, the signal is increased more than the background. This makes 
it interesting to look at a coincident spectrum of burst samples, i.e. a quasi- 
diffuse flux. The total flux $ is given as 



n 



(115) 



Three different burst samples are used here for a systematic investigation: 
® BATSE sample: the energy spectrum of 568 bursts from the BATSE 



G+04b 



G+04b 



The ne utrino spectra for the individual 
and in B"'"06l ] using different parameter 



B+06 



FROO for a 



catalog are given in 
bursts was derived in 

settings. Calculations here are following the latter publica tion 
Redshifts are determined using the variability method, see 
description. 

KoNUS sam ple of sh ort bursts: a catalog of 130 short bursts has been 
presented in |M"'"05a| . of which 99 have fitted energy spectra and 34 are 
localized. Redshifts are set to ^ = 0.1, if unknown . 



@ Swift sample of long bursts: mea sureme nts from Swift-BAT have been 
published on the Swift web-page |Swi08[ | . In this sample used here, those 
188 long bursts with published spectral information until March 12, 2007 
are included. For the Swift bursts without measured redshifts, 2; = 2 is 
assumed. 

To compare the different samples to each other, the coincidence fluxes of the 
samples are shown in Fig. 24. Three major conclusions can be drawn from the 
comparison of coincidence spectra. 

(1) While the spectral shape for the three samples of long bursts are rela- 
tively comparable, even small differences in the spectral indices result in 
order of magnitude deviations in certain parts of the spectrum. The main, 
systematic difference between the samples is the low normalization of the 
Swift sample. Here, the reason lies in the properties of the BATSE and 
Swift experiments. While BATSE had a wide fleld of view and could 
observe about 2 vr sr, its sensitivity was not as good as the one of Swift. 
This fact leads to the observation of a large number of strong GRBs for 
BATSE (~ 1 per day) compared with a lower number (~ 100 per year) 



Short bursts are believed to occur when neutron-star /black hole binaries merge. 
These events typically happen in old regions, and the mean redshift is therefore 
expected to be as low as ^ = 0.1. 
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Fig. 24. Summary of the coincidence spectra of three GRB samples. Bursts in the 
BATSE sample ® are typically quite strong. SwiFT-detected bursts (sample (D) 
have a flat spectrum and a high second break energy, while the short KONUS bursts 
(sample ®) are generally very weak and produce a steep spectrum. All calculated 
mean spectra differ significantly from the average diffuse prediction made by Wax- 
man and Bahcall (dashed line). 



of bursts with lower fluence for Swift. Consequently, the neutrino flux 
which is proportional to the photon fluence is also lower. The differences 
in the field of view are corrected by assuming a GRB rate of 1000/yr for 
all of the samples. The coincidence spectra of the two BATSE-samples 
are more than two orders of magnitude more intense than the Swift 
sample only due to the fact that the increased number of bursts in the 
Swift sample is not due to a large field of view but due to a good sensi- 
tivity. There are fewer strong bursts but a large number of low-intensity 
ones which reduces the total contribution of the signal. Such a bias from 
experimental data will be present in all of the experiments, since the re- 
sults in each case will be dependent on properties such as energy range, 
sensitivity, field of view etc. 
(2) The energy spectrum for short bursts is about one power steeper than 
the general flux for long bursts. This general conclusion is observable in 
the coincidence spectra. The neutrino spectrum from KoNUS bursts is 
much steeper than the ones from long duration bursts. The reason lies in 
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the hardness of the photon spectra. Since the short keV photon spectra 
are typically much harder than the long burst spectra, neutrino spectra 
for short bursts tend to be softer (7 and v energies are inversely propor- 
tional). This underlines that long and short bursts have to be treated 
separately. 

(3) For all samples both the average and the coincidence spectrum differ from 
the spectrum that is considered as the standard diffuse spectrum. 

Due to the variation in the single source spectra, the detection rate for single 
bursts in a neutrino detector like AMANDA or IceCube scatters over almost 
six orders of magnitude. There are a few bursts with close to one event per 
km^, while the mean of the distribution rather lies around i? ~ 0.01 yr~^. The 
values of the mean for each sample are given in table 11. 

Only 34 of the 99 sources with spectral information in the KoNUS satellite are 
well-localized by the IPN3. Consequently, the detection rates are only given 
for these 34 events. The single source parameterization causes a much wider 
scattering than the WB rates. Also, the mean detection rate is about half an 
order of magnitude lower in the case of single source parameterization (see 
table 11). The main reason lies in the steepness of the actual GRB spectra. 
These are typically about one power steeper than assumed in the WB spectrum 
and have a relatively low flux at the main detection energy range in neutrino 
Cherenkov telescopes [Ey > lO'* — 10^ GeV). This demonstrates again that at 
least systematic differences in the spectra need to be taken into account. For 
the 165 long bursts with given localization, the results from previous samples 
can be confirmed. Here, the systematic shift between WB rate and single 
source parameterization comes from the lower normalization of SwiFT bursts. 

While the detection of a single burst will be difficult even with instrumented 
arrays of one cubic kilometer like the ones planned for IceCube or KMSNeT 
(the most intense bursts in the detector have a detection probability of < 
1 neutrino event per km^ effective area), the stacking of photon-detected GRBs 
increases the significance immensely, so that the model can be tested in the 
very near future. GRB analyses usually make use of the temporal and spatial 
information from the GRB satellites. GRB030329 was analyzed and it turns 
out that the single source sen sitivity in this case is several orders of magni- 



tude above the predicted flux |S"'"05d |. The burst was analyzed adjusting the 
spectrum to the observed properties of the GRB, similar to how it is described 
above. An alternative approach like the one done by AMANDA, is the stack- 
ing of all bursts in a sample using the WB spectrum for al l of the bursts. 
By sta cking da ta from the sample of neutrino-induced muons AI"'"07cl | or cas- 



cades |AI^07al | according to the arrival direction and their arrival time, the 
sensitivity to a diffuse flux from GRBs according to WB can be determined. 
An analysis looking for the temporal cluste ring of e vents within small time 



bins (100 s) is done for cascade-like events AI"'"07a . Such a search is inde 



pendent of the photon-detection of GRBs and therefore also includes choked 
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bursts and bursts which are missed by GRB sateUites. 



7.3 Afterglow neutrinos 



For the afterglow of a Gamma Ray Burst, it is typically assumed that the 
emission arises when the burst hits the interstellar medium. The reverse shocks 
accelerate charged particles and electrons emit synchrotron radiation at X-ray 
energies. The spectrum is observed to decay as 

^ oc Ex-'^- . (116) 

The characteristic synchrotron frequency at which the synchrotron cooling 
time becomes comparable to the ejecta expansion time, e\ lies below keV- 
energies, 

e\ ^ 0.3 e-^'Ll n^' E^^'^ t7,^'^ keV . (117) 
Here, E =: £"53 ■ 10^^ erg is the isotropic, kinetic energy of the fireball and 
65 ^2 := ■ 10~^. At energies Ex > ^x, the spectrum steepens 



dN 



dE 



- oc Ey-^"" (118) 



X 



and Px = Cix + 1/2. On average, the spectral indices scatter around ax ~ 3/2 
and /^x ~ 2. Just as for the prompt emission, individual bursts can have spec- 
tral indices deviating significantly from the average value as it is investigated 
for the Swift sample later this section. 



Provided that proton and pion synchrotron losses can be neglected I. the 
afterglow neutrino spectrum is then given as 

E,'-^ = A:-{ " " (119) 

The normalization A^, can be determined by assuming that UHECRs are pro- 
duced by GRBs, implying that the power in UHECRs is comparable to the 
electromagnetic output from Gamma Ray Bursts, 



/ EcR ^ dEcR ^ 10'' erg MpC ^r"! . (120) 



The neutrino break energy is inversely correlated to the photon critical energy 
e\ as described in Equ. (106) with 0.1 — 1 keV. Therefore, 



ef =0.7-7- 10^(1 + ^)-' GeV (121) 



Reverse shocks happen at radii of ~ 10^"^ cm, which is large enough to avoid 



synchrotron losses WBOO(]. 
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for a boost factor of F ~ 300. The spectral indices are given as a^^a = —f^x + 1 
and Pi^^a = — /3x + l-5. With an average value of Px ~ 2, the neutrino spectrum 
behaves as 



E,^ ioxE, < el 
E,^'^ forE, >e^ 



(122) 



The break energy is fixed in WBOOl ] to = 10^ GeV, and the normalization 
factor is calculated to be 



A„ = 10~^° GeV cm-2 s"^ sr"^ . 



(123) 



The maximum energy is again given by the maximum proton energy with 

£;max _ ^q21 ^^^^ ^max ^ £;max/20 ^ 5 ■ 10^ GeV. 



7.4 Single source spectra 

The result as presented above relies on the averaged parameters as observed 
for GRBs. Single source spectra can in some cases deviate significantly from 
the m ean sp ectrum. SwiFT-XRT-detected bursts are given on the Swift web- 
page jSwiOSf (as of August 2007). Since XRT is sensitive in the energy range 
of (0.2, 10) keV. Single powerlaw fits are presented and for the break energy, 
the lower energy threshold ~ 0.2 keV can be used, resulting in a neutrino 
break energy of 

^,^sn,^ft ^ 3 5 . . (1 + z)-2 GcV . (124) 

For those bursts without given redshift, z = 2 is used. The distribution of 
spectral indices is shown in Fig. 25. The mean value lies near f3x ~ 2. Figure 26 
shows the spectrum for each individual GRB in the sample. The normalization 
calculates as it was done for the prompt emission, assuming that a single burst 
produces a fraction of the total flux of charged Cosmic Rays, resulting in the 
neutrino normalization factor of 

11 F 

Again, ^ 0.2, but this should be considered as an upper limit, since 
the efficie ncy for pion pro duction in afterglow emission can be reduced sig- 



nificantly [WBOOl . lB+07bl |. The maximum energy for neutrino emission is 
^^max ^ lo^-^ GeV, given that protons are accelerated up to 10^^ GeV, and 
neutrinos carry l/20th of the proton energy. Figure 26 shows that the single 
source spectra (dotted lines) deviate from the mean spectrum (solid, thick 
line) by up to two orders of magnitude. Furthermore, the spectral index varies 
from burst to burst by more than one power. 
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7.5 Neutrinos and early afterglow 



As discussed in Section 3, the rapid slew of XRT after GRB alerts from 
BAT, both on board of the Swift satellite, made early afterglow detection 



possible. It is discussed in |Mur07l | that in the case of late internal shocks, 
neutrino production via proton-photon interactions can be very efficient, while 
reverse shocks cannot produce high-energy neutrinos. In particular, short X- 
ray flares as observed by Swift can be connected to neutrino emission in case 
of had ronic par ticle acceleration within the frame of the late internal shock 



model |MN06bl ]. For neutrino telescopes like IceCube, a diffuse flux from X- 



ray flares can be tested, the flux level is near the IceCube sensitivity level for 
one year of detection {dN^/dE^ ■ Ej^ ~ lO"^-^ GeV cm'^ s'^ sr-^ at 10^ GeV). 
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Fig. 25. Histogram of the afterglow spectral index ax for long Swift bursts. 



7. 6 Neutrinos from GRBs: presence and future 



Neutrinos can originate in the three phases of a GRB. Neutrino emission is 
expected prior to the G RB, when the pre-GRB interacts with the central 



pulsar or the SNR shell RMW03I ] . The prompt emission is accompanied by 
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Fig. 26. Single source spectra from SwiFT-detected GRB afterglow s (dott ed lines). 
The solid, thick line indicates the mean spectrum as calculated in WB00l |. 



a neutr ino flux if the o bserved spectrum of charged Cosmic Rays comes from 



GRBs |WB97l IWBQQj . The after glow em ission can for the same reason be 



accompanied by a neutrino signal 



WBOO 



Searches for neutrinos from GRBs have been performed by the AMANDA 
experiment. In the case of muon neutrinos, directional information provided by 
IPN3 was used in order to reduce the atmospheric background. No signiflcant 
GRB signal was dete cted. The neutrino flux limit from AM ANDA data as 



presented in AI+OTcl ] comes close to the flux as predicted in WB99 



$DL ^ g . ^^^-2 g-i 



(126) 



for a WB-like spectrum. With IceCube, it will be possible to constrain the 
model. In the case of cascade-like events, two analyses have been performed. 
The flrst one is also based on the selection of events according to IPN-3- 
detected bursts. The second one does not use photon-based measurements, 
but simply looks for unexpected temporal clustering of events. The cascade- 
limits are still far above the predictions and no general conclusions can be 
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drawn yet. The RICE experiment has investigated five bur sts with respect to 
a possible flux connected to the afterglow emission B"'"07bl |. with limits a few 
orders of r nagnitu de above the prediction, derived for each burst individually, 
following IWBOd . With the stacking of more bursts, it can be possible to 
improve such limits. Larger detection arrays like AURA or SalSA will have 
increased effective areas and will hopefully be possible to detect neutrinos 
from GRB afterglows. 

The investigation of individual burst fluctuations has revealed that the detec- 
tion rate^jactua^es within each class between ~ 10~^ — 1 events/km~^/burst, 
see ™ 



G+04bl . IB+Od lBec07 



The systematic variations between the flux pre- 
dictions for different burst samples come from the bias connected to each of 
the satellites. The energy band pass and the fleld of view of the detector are 
important features which lead to a selective sampling of events. For instance, 
when comparing BATSE and Swift bursts, the latter are on average less lu- 
minous, since SwiFT has a smaller fleld of view and a better sensitivity than 
BATSE. Short bursts detected by Swift typically do not belong to the class 
of short-hard bursts. This is a new su b-GRB-class detected by Swift, which 
has not been observed before jSaSOd ]. Swift can detect these events due to 
its low energy band pass of (15, 150) keV. 

In the future, GLAST will open the possibility of observing the high-energy 
component of GRBs, probably revealing new classes of GRBs. TeV photon 
detectors like MILAGRO and MAGIC can observe the component of very 
high-energy photons in coincidence with GLAST in orde r to iinprove the 
unders tanding of high-energy emission from GRBs, see e.g. [VMOTl . lGM+07al . 



Con07l |. It is likely that the differences between these classes are intrinsic 
as they are for long and short bursts. The detection of neutrinos in the three 
different phases can contribute to distinguishing the different sub-classes and it 
can also to explain the underlying physics. With the completion of IceCube in 
2011, flrst results can be expected, constraining current neutrino flux models. 
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Table 10 

Parameters for the prompt photon and neutrino emission in GRBs. The numbers 
quoted as typical values for GRBs, XRFs, Short Hard Bursts (SHBs) and Short 
High-Energy (SHE) bursts are to be taken as rough bench marks, since all of these 
parameters fluctuate strongly as emphasized in the text. The calculation of the 
luminosity distance at a redshift between 2 = 1 — 2 for long bursts and z = 0.1 
for short events is done using cosmological parameters of Vtm = 0.3, VL\ = 0.7 and 
h = 0.72. 
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mean rate/GRB [km ^] 


total number of i^^j [km ^] 


BATSE 


2^Q-2.05±0.67 


14 


WB [BATSE] 


-[^Q-1.86±0.23 


9.1 


KONUS 


]^Q-2.76±0.69 


0.28 


WB [KONUS] 


■|^Q-2.36±0.17 


0.17 


Swift 


2^Q-2.72±1.35 


3.4 


WB [Swift] 


]^Q-1.95±0.17 


2.0 



Table 11 

Mean neutrino spectra parameters for the BATSE, KONUS and Swift samples. 
The standard deviation to the mean values was calculated as an error estimate. 
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8 Neutrinos from starburst galaxies 



The emiss ion of neutrinos from starburst galaxies was suggested by Loeb & 



Waxman |LW06| (LW). The observation of synchrotron radiation at radio 
wavelengths implies the existence of relativistic electrons, which have a charac- 
teristic energy of ~ 1 GeV, assuming that magnetic fields in starburst galaxies 
exceed the galactic magnetic field by two orders of magnitude. If protons are 
co-accelerated along with the electrons, neutrinos can be produced in proton- 
proton interactions. For high densities of tt, > 100 cm~^, starburst galaxies 
serve as a "proton calorimeter", in which all protons are absorbed in pp in- 
teractions to produce neutrinos. There are three crucial parameters for the 
calculation of the neutrino spectrum from starburst galaxies: 

(1) The spectral behavior - It is assumed by LW that the neutrino spectrum 
follows the induced proton spectrum, leading to a spectral index of a,^ ~ 
ftp ~ 2.00 — 2.25. The maximum energy for the neutrino spectrum is 
given by EJ^^^ = 0.1 PeV, derived from the fact that starbursts serve as 
proton calorimeters up to a proton energy of E'p < 3 PeV. Here, the ratio 
between proton and neutrino energy is Ep/Ei, = 20. 

(2) The normalization of the total spectrum - In order to estimate the inten- 
sity of the neutrino fiux, the synchrotron radiation from electrons is used. 
This can be done using the following connection between the different lu- 
minosities: The observed synchrotron emission is directly connected to 
the energy flux of the primary electrons, 

2 dNe 

Eg ~ 2 h' synch Ly^gynch • (-^^ ' ) 

ahie 

Here, dNg/dEg is the energy spectrum of the primary electrons with an 
energy Ee, v synch is the characteristic frequency of the synchrotron pho- 
tons and Ly^synch is the synchrotron luminosity per frequency interval. 
The factor 2 arises from the fact that v synch c< -Eg, so that the energy in- 
terval AE'e corresponds to twice the width of the synchrotron frequency 
interval, see Section 3.1.3. Assuming a fixed ratio of protons to electrons, 
^p/e = 6 for starbursts, the proton luminosity is connected to the electron 
luminosity as 

„ dLrt „o dNp , . 

If all protons are absorbed in pp— interactions, the proton luminosity 
corresponds to the neutrino luminosity, with a reduction factor 1/3, ac- 
counting for the fact that 2/3 of the A— resonance processes goes into 
charged pions and 1/2 of the charged pion's energy goes into neutrinos. 
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Here, / accounts for the fraction of starbursts serving as a calorimeter. It 
is assumed by LW that all starbursts are calorimeters, resulting in / = 1. 
Combining Equations (127), (128) and (129) gives a connection between 
synchrotron radiation and neutrino emission as 

The normalization is performed by using the observed synchrotron spec- 
trum at V synch = 1.4 GHz, which corresponds to an electron energy of 
Ee = 1 GeV. The value is calculated by using the correlation between ra- 
dio and FIR emission. Here, it is assumed that the entire FIR background 
arises from the emission from starburst galaxies. 
(3) The source evolution function with redshift- The distribution of starburst 
galaxies in the sky is taken i nto acc ount, using recent models of star 



formation rate as discussed in LW06 and T"'"06 



The resulting neutrino flux is given as 

^<'°'=^) 'I^-^-'"- 1^ ''^"l..,....=i.4GHz (131) 

assuming that the characteristic electron energy corresponds to the neutrino 
energy, E^ ^ E^ ^ 1 GeV. The factor ^ 0.5 takes into account the source 
evolution which follows the star formation rate and tn is the age of the Uni- 
verse. This flux gives the total muon-neutrino signal for a flavor ratio at the 
source. If, on the one hand, a mixing of neutrinos from (z/g, z/^, u^-) = (1:2:0) 
at the source to (z/g, u^, z/^-) = (1 : 1 : 1) is assumed, the muon-neutrino flux 
must be divided by a factor of two. On the other hand, if muon-neutrinos 
and anti-muon neutrinos are considered, the flux is enhanced by a factor of 
two. Thus, the flux which can be measured with large volume neutrino de- 
tectors in the muon-neutrino channel corresponds to the result of Equ. (131). 
Numerically, the neutrino flux is given as 



„ r, dN,, 7 / E,, 

eJ — ^ = 10" 



dE, VlGeVy 



X 2-ap 

:) GeV cm-2 s-^r-^ (132) 



An signal is already close to recent AMANDA limits AI"'"07d| . If the 

induced proton flux is steeper, i.e. ap ~ 2.15 as suggested by Cosmic Ray data, 
the flux is more difficult to observe with AMANDA, and next generation 
experiments like IceCube and KM3NeT are needed to test the model. 

The question on the fraction of the diffuse far infrared (FIR) flux coming from 
starbursts is still unclear. While LW a ssume that 100% of the detected signal 



comes from starbursts, Stecker |Ste07l | presents that the fraction is only 23% 



on average. On the whole, the diffuse flux from starbursts is probably lower 
than predicted: LW assume that starbursts are loss dominated, which means 
that most primaries interact and do not escape the source. This enhances the 
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neutrino flux, since basically all protons lose their energy in proton-proton 
interactions and produce neutrinos. This fact is reflected in the variable / 
which, for the loss-dominated case of a proton calorimeter is / = 1. Observa- 
tions of the spectral radio index of the sources {S ~ u^^-^) indicate, however, 
that starbursts are in the diffusion limit, indicating that a negligible fraction 
of protons interact and only few neutrinos are produced. With IceCube soon 
reaching the sensitivity level of the prediction of LW, it will be possible to 
solve this question. 

There is another possibility to expect enhanced neutrino emission from star- 
bursts. In the past few years, it could be shown that long GRBs are typically 
connected to the explosion of Wolf-Rayet stars into a supernova Ic. These oc- 
cur preferably in star forming reg ions. Thus, a diffuse flux of GRBs similar to 



the prediction of |WB97l . IWB99l | should originate from the direction of these 
galaxies. There are two different ways to normalize the diffuse ORB spectrum. 
One method is to assume that the observed keV-photon flux is proportional to 
the neutrino flux. In that case, the normalization is dependent on the number 
of observed GRBs per year. This number is strongly dependent on the instru- 
ment and the number is not very exact. Under the assumption that GRBs 
accelerate protons up to the highest energies, Ep ~ 10^^ eV, the neutrino 
spectrum can also be normalized to the flux of ultra high-energy cosmic rays 
(UHECRs). In this case, the normalization is independent of GRB observa- 
tions. It should be kept in mind that the spectral index of the spectrum still 
varies from burst to burst - in the model of Waxman&Bahcall, an average 
spectral index was used. 

It is possible to look for a neutrino signal from GRBs by stacking starburst 
galaxies. This method has one advantage over a triggered-GRB search: it is 
a systematic search, since independent of GRB data. A disadvantage is that 
only nearby events can be included, since the sample of starbursts only reaches 
out to redshifts of z = 0.07. It should, however, be possible to use IRAS data 
to identify starburst galaxies at higher redshift. The search for a GRB signal 
from starbursts should be considered as a systematic search for choked and 
undetected GRBs. The sources can be selected according to their FIR-flux, 
since this is a measure of the SN rate in a starbursts. A higher FIR flux 
indicates a high star formation rate, thus more SNe and therefore, also more 
GRBs. 
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9 Cosmogenic Neutrinos 



The prediction of the GZK cutoff at energies in the spectrum of UHECRs due 
to proton interactions with the CMB implies the existence of a flux of ultra 
high-energy neutrinos. The production scenario and the expected event rates 
of such a flux of cosmogenic neutrinos are reviewed here. 

9. 1 Production of cosmogenic neutrinos 



High-energy protons interact with the CMB on their way to Earth via Bethe- 
Heitler pair production, 

P7cA/iJ ^ pe+e" (133) 

and via the A— resonance, 

P7CA/B^A. (134) 
The energy threshold for Bethe-Heitler pair production is = 5 ■ 10^^ eV, 
while it is higher for the A— resonance, Ef^ ^ 5 • 10^^ eV. The attenuation 
length A = Ep/ {—dEp/dx) determines the propagation radius of the protons 
from a given source, with Ep as the energy of the proton and dEp/dx as the 
mean energy loss rate. The latter depends on the photohadronic cross sections 
as well as on the CMB photon field density, and is determined by Monte 
Ca rlo meth ods. Figure 27 shows the simulated attenuation length as presented 



m 



HTS05|. Below 10^^ eV, Bethe-Heitler pair production is the dominant 



process, before the production of the A— resonance becomes dominant above 
its threshold energy, Ep > 5-10^^ eV. Pair production results in an attenuation 
length of A ~ 1000 Mpc above 10^^ eV, while the A— resonance allows for 
the undisturbed propagation over a range of only ~ 10 Mpc. The latter is 
therefore the dominant process leading to the prediction of the so-called GZK 
cutoff at the highest energies, i.e. Ep > 5 ■ 10^^ eV. 

Neutrinos are produced in the decay of the A— resonance via the channel of 
charged pion production, see Section 3. The resulting diffuse neutrino flux is 
usually referred to as the flux of cosmogenic neutrinos or GZK neutrinos. The 
neutrino flux ranges from approximately E^, ~ 10^^ eV up to E^, ~ 10^^ eV. 
The shape and strength of the neutrino flux depends on various parameters 
which are described in the next paragraph. 



9.1.1 The parameter space 

The strength and shape of the cosmogenic neutrino flux depends on seven 
factors: 

(i) the energy spectrum of primary protons - spectral shape, normalization and 
energy cutoff, 

(ii) the proton-photon cross section ap^, 

(iii) the CMB photon density n^'^'^, 
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Fig. 27. Attenuation length as a function of energy. The straight line represents 
adiabatic energy losses due to the expansion of the Universe. The upper (green) 
hue shows energy losses due to Bethe-Heitler pair production, while the lower 
curve is the a ttenuation length resulting from pion-production processes. Reprinted 
from [HTS07|, with permission from Elsevier. 



(iv) the evolution of proton-emitting sources with redshift, 

(v) the composition of UHECRs, 

(vi) magnetic field deflections, 

(vii) neutrino oscillations . 



9.1.1.1 Primary spectrum The spectral shape and normalization of the 
sources determines the cosmogenic neutrino flux. Typically, it is assumed that 
all sources follow a powerlaw with an exponential cutoff at a critical energy 

Er, 



dNp 



exp 



E.. 



(135) 



Such a spectrum can be produced in shock fronts of AGN or GRBs, see Sec- 
tion 3. Monte Carlo simulations indicate a spectral index close to Up ~ 2. 
A generic spectrum with an exact value of a = 2 is therefore used for all 
sources in calculations of the cosmogenic neutrino flux. However, it must be 
noted that recent results from test particle shock acceleration indicate that 
this value can deviate from = 2 dep ending on the exact shock configu- 



ration, see e.g. [Bar04l . lMBQ07bl . |SBS07|. The spectral powerlaw behavior is 



in accordance with the powerlaw structure of the UHECRs, which are used 
for the normalization of the total proton spectrum at Earth. This way, no 
concrete normalization for the single sources needs to be given. 
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9.1.1.2 Cross section Cross section measurements from accelerator data 
are used in order to estimate the fraction of energy going into neutrinos and the 
number of neutrinos produced in one interaction. This is done for the complete 
energy spectrum at the given redshift z, and the spectra are propagated from 
the source to Earth, using a large number of discrete steps of Az. The value 
of Az is chosen to be much smaller than the average mean free path. 



9.1.1.3 CMB photon density The photon density of the cosmic mi- 
crowave backgrou nd, n^'^ ^^, is well known since its observation with the 

cobeT^ 
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satellite jMi90|]- The spectrum can be fitted with a blackbody spec- 
trum of T ~ 2.7 K, 



n 



CMB ( 



7 



(z/) OC 



exp 



hv 



(136) 



9.1.1.4 Evolution of source population with redshift Typically, it 
is assumed that ultra high-energy protons are produced in AGN or GRBs, 
see Section 3. These objects are believed to follow the evolution of star for- 
mation and the redshift behavior i s often assumed as p ~ (1 + z)*" up to a 
redshift z = 2max, with m > 0, e.g. |Sch72l | . At high redshifts, 2 > 3 to 4, the 
source number counts become small and the exact behavior of for example the 
distance distribution of AGN is difficult to determine. 



9.1.1.5 Composition of UHECRs In most neutrino flux predictions, 
it is assumed that a dominant fraction to the spectrum of UHECRs are pro- 
tons. While protons are very efficient neutrino gener ators v ia p 7— interactions. 



heavier nuclei lose energy via photo-disintegration [Ste69|. For instance, iron 



nuclei interacting with CMB photons induce the following reaction: 

''Fe7cAfB^''Mnp (137) 
^^Mn 7c7MiJ Mnn (138) 
n pe^ Ue . (139) 

In these processes, only very few neutrinos are produced, i.e. via neutron 
decays. Most models neglect the fraction of heavy nuclei in the spectrum of 
UHECRs. However, if heavy nuclei dominate the Cosmic Ray spec trum at the 



highest energies, the neutrino flux can be small as pointed out in (HTS05 . 
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COsmic Background Explorer 
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9.1.1.6 Magnetic field defiections Protons are deflected by intergalac- 
tic magnetic fields and by tlie magnetic field of the Milky Way on their way to 
Earth. While this effect could play a role, it is usually neglected in the calcu- 
lation of cosmogenic neutrino flux predictions. The modelling of extragalactic 
fields is very difficult, since they are not well- determined. 



9.1.1.7 Neutrino oscillations Neutrino oscillations are usually not con- 
sidered in the calculation of the GZK neutrino flux. For a precise determination 
the changes in the spectrum, oscillation effects would need to be applied af- 
ter each single propagation effect. However, a good approximation is add up 
the z/^ and spectra at Earth and assume that each of the three flavors are 
present in equal amounts, as it is the case if the neutrinos are produced far 
away from Earth. 



9.1.2 On the evolution of the GZK neutrino flux model 

The exist ence o f the cosmogenic neutrino flux was proposed by Berezi nsky fc 



Zatsep in |BZ69j shortly after the predict ion of t he GZK cutoff was made |Gre66 



KZ68| . The first prediction presented in |BZ69l | was revised during the years by 



applying the gained knowledge on the different parameters described above. 
Table 12 summarizes the development of the model. Each row includes at least 
one new development in physics or astronomy which was implemented by one 
or more authors: 



1969 
1972 

1985 



When the first prediction was made in 1969, the mass of the py^— bosons, 
M^yi, was not yet determined experimentally. 

The implementation of the H^— boson propagator with the correct mass term 
in the calculation of the c osmogenic neutrino flux lead to the improvement 



of the prediction, see e.g. |WTW72l . ISte73l . IBS75 



Hill & Schramm normalized the spectrum according to recent measurements 
of UHECRs. At the same time, the cognition that sources are not homoge- 
neously distribute d over the the s ky, but evolve with redshift was considered 



for the first time HS85. HSW86 



1991 The evolution scenario, cross sections and normalization factors were further 
improve d using recent resul ts from astrophysical and particle physics data. 



S+91, YT93. PJ96 



2001 



2005 



2007 



see e.g. 

The observation of UHECRs above 10^° eV lead to the prediction that 
UHECRs are accelerated up to at least -Emax ~ 10^^ eV. The change in 
the maximum energy of the primary sp ectrum also changes the cosmogenic 



neutrino flux wh ich was pointed out in [ESSOl 



It was shown in |HTS05l | that the cosmogenic neutrino flux is reduced sig- 
nificantly if a large fraction of heavy nuclei is present in the observed flux 
of UHECRs. 

It is not clear yet which are the sources of UHECRs. The two primary can- 
didates are AGN and GRBs. These source classes follow different evolution 
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ref 




Ap 
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B- 


photo- 
















dn. 


dis. 


1969 


[BZ69] 


unknown M|y± 


estimate 


iO 


eV 


no 


no 


no 


'72-75 


e fWTW72l 
[Ste73, BS751 


inrl A'fiTz-\- 

±X±\^Xa JIJ. M/ HZ 


pstimatp 


1020 


eV 


no 


no 


no 


'85-'86 


[HS85], 
1 no vv ou 1 


incl. Myj/i 


CRs 




eV 


yes 


no 


no 




P o- fS+911 
e.g. 1 J 

[YT93 RB931 


je& 


iixipr • 


1 n20 
10 


eV 


imp)i*a 




n n 


2001 


e.e. fESSOll 


yes 


CRs 


1021 


eV 


yes 


no 


no 


2005 


fHTS05l 


yes 


CRs 


1021 


-22.5 eV 


yes 


no 


yes 


2007 


[YK071 


yes 


CRs 


1021 


eV 


impr. 


no 


no 



Table 12 

Evolution of GZK models. Boldface: improvement compared to previous models. 
Impr.=Improved; dfl.=deflections; dis.=disintegration. 



scenarios, which is pointed out in YK07l | . GRBs evolve stronger than AGN 
which leads to an enhanced signal of cosmogenic neutrinos. 

In Fig. 28, the energy spectrum of cosmogenic neutrinos is shown for three 
different models. The uppermost, solid line shows the resulting spectrum if 
GRBs are responsible for the flux of UHECRs. GRBs are observed to evolve 
even stronger than dif ferent A GN classes, which leads to the enhancement of 



the flux as shown in jYK07|. The current AuGER limit comes close to the 
strong evolution model and with data from the completed Auger array, it 
will be possible to test this model in detail. The middle, dashed lin e show s 



the spectrum for an AGN-type evolution of sources as calculated in RB93 



The lower, dot-dashed line represents the hyp othetical spectru m using that 
UHECRs consist of iron only as calculated in HTS05I . iTayO?! ] to show that 
the cosmogenic neutrino flux is decreased significantly with increasin g contri- 



butions from heavy nuclei in the Cosmic Ray spectrum. It is shown in |HTS05 



that, the higher the mass number of the nucleus, the less neutrinos are pro- 
duced during propagation. A maximum energy of Ep = 10^2-4 QeV is used 
for the primary iron nuclei, taking into account the increase of the maximum 
energy with the charge Z, see Equ. (12) in Section 3. For lower maximum 
energies, the flux is further diminished. 
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Fig. 28. Comparison of different GZK neutrino flux models, for all neutrino fla- 
vors, neutrinos and anti-neutrinos. Uppermost, solid line: prediction using an evo- 
lution functi on mat ching the one from GRBs YKOTJ. Middle, dashed line: model 
presented in RB93I ]. using an evolution f unction matchi ng those of AGN. Lower, 
dot-dashed line: pur e Fe-based spectrum jHTSOSl . lTav07l | . The model-independent 
limits from and ANITA-LITE are shown tog ether wi th the sensitivity for 

a 45 day long flight with the full ANITA experiment 



9.2 Comparison of event rates in underground detectors 



In this section, the expected neutrino rates is calculated for underground de- 
tectors in order to get an estimate of the necessary array size for the detection 
of ultra high-energy neutrinos. Four detector properties have to be considered: 

(1) The threshold energy -Emin of the detector determines the total integral 
neutrino rate per year and km^. The threshold energy is varied between 
Emin = 10^ GeV and -Emin = 10^° GeV in the following calculations. 

(2) The depth of the detector below the Earth's surface d determines the 
fraction of neutrinos absorbed before reaching the detector. 

(3) The effective range of neutrino- induced cascade, r^ff. It is determined 
by the size of the detector, since the cascade has to be induced within 
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the detection array in order to be observed. All following calculations are 
performed using r^ff = 1 km. For a constant effective range, the results 
scale linearly with the range a simple factor can be used to account for 
different detection array sizes. 
(4) The density of the detector material, p' = p/ PH2O has to be considered 
as well. The density will be given in terms of p' = 1, assuming water or 
ice as the detector medium. For detection arrays in salt, the result can 
be multiplied by a factor 2.3 in order to receive the corresponding event 
rates. 

The angle integration was performed over the upper hemisphere (2 vr) only, 
since at these energies, no contribution from the lower hemisphere due to 
absorption. The result is given in units of yr~^ km^'^ GeV^^. Three models 
are compared here. 



[YK07] GRB evolution scenario YKOTj]. With a source evolution matching the one 
of GRBs, a high neutrino flux is expected. This model is used as the maxi- 
mum flux possible, since it is close to, but still below, most recent neutrino 
flux limits. 

AGN evolution scenario 



[RB93] 
[HTS05] 



RB93 



The less strong evolution of AGN leads to 
a reduced neutrino flux co mpared to the G RB evolution scenario. 
Fe iron nuclei as primaries HTSOSl . iTayOTi ] . It is assumed that the dominant 
fraction of Cosmic Rays consists of heavy nuclei. A pure iron spectrum is 
an absolute lower limit - the neutrino flux increases with the decrease of the 
mass of the primary nucleus, and it is likely that the spectrum of UHECRs 
is a composition of light and heavy nuclei. The question of the fraction of 
heavy nuclei in UHECRs is not answered yet. Therefore, this model is used 
as a lower limit here. 

Since the cross sections for nucleon neutrino interactions do not change sig- 
nificantly with neutrino flavor, the results are basically the same for both 
flavors. The calculations presented here are performed for electron neutrinos. 
For muon neutrinos, the results are comparable. In the case of tau-neutrinos, 
the regeneration of tau leptons needs to be considered as well. Figure 29 shows 
the integrated electron neutrino detection rates depending on the threshold 
energy of the detector for the four different GZK neutrino flux predictions. 
A depth oi d = 0.5 km is assumed in these calculations. With a detection 
threshold below -Emin ~ 10^ GeV, a rate of i? ~ 1 events/yr/km^ can be ex- 
pected for the maximum model [YK07]. The AGN evolution model [RB93] 
yields a factor of ~ 5 less neutrinos and the minimal model using Fe nuclei as 
primaries results in ~ 0.03 events/yr/km'^. This shows that the instrumenta- 
tion of volumes significantly larger than 1 km^ is eligible in order to achieve 
a significant rate of GZK neutrinos per year. This is possible with radio and 
acoustic detection arrays, since the cascades have a large extension, so that a 
sparsely instrumented area suffices. Simulations for the extension of IceCube 
to IceCube Hybrid have been performed, assuming a radio/acoustic array of 
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10 X 10 km^ in the ice around the IceCube detecto r, together with additional 



IceCube strings around the original array [VjjOSaJ. The spacing between the 
acoustic/radio strings is assumed to be 1 km. While the size of the cascades in 
the radio and in the acoustic channel is an advantage of these new techniques, 
the high energy threshold at -Emin = 10^'^ GeV for radio measurements, and 
at -Emin = 10^'^ for acoustic detection of neutrino-induced showers diminishes 
the rate by up to an order of magnitude as it can be seen from Fig. 29. The 
total event rate per year is still a few up to 10 events per year for an array of 
> 100 km'^ and an AGN evolution scenario. 

Table 13 shows how the event rate for electron neutrinos changes with the 
depth of the detector below the Earth's surface. Results for five different 
depths, d = 0.4, 0.6, 1.0, 1.5, 2.5 km, are shown for a threshold energy of 
Emin = 10^ GeV. The rate decreases with the depth. It is reduced by a factor 
of two when the detector is built at 1.5 km depth instead of 0.5 km. For thresh- 
old energies above -Emin ~ 10^ GeV as they are expected for radio and acoustic 
detection, the background of atmospheric muons and neutrinos is negligible. 
In the case of optical detection, energy cuts can be applied in order to remove 
lower-energy atmospheric events. Consequently, it is not necessary to deploy 
the detector deep into the Earth to reduce any background and a detector has 
the best detection potential at shallow depths of ~ 0.5 km. 
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Fig. 29. Integrated neutrino rate, representing the number of detected electron neu- 
trino induced cascades above a certain energy Emm- Three different models are 



shown. The solid line represents the model from [YK07| | . using G RBs as the sources 
of UHECRs, the dashed line shows an AGN-based model fr om RB93|. Th e dot- 
dashed line represents a pure iron spectrum as calculated in jHTSOSl . lTav07l | . 
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Table 13 

Electron neutrino event rates for different detector depths, d = 
0.4,0.6,1.0,1.6,2.5 km. The detector threshold energy is Emm = 10^ GeV. 
The rate decreases with the depth, since neutrinos are absorbed in the Earth. 

9.3 Conclusions and outlook 



Shortly after the prediction of the GZK cutoff, the first estimates of the cos- 



mogenic neutrino flux were made [BZ69|. During the past 40 years, advances 



in particle physics, e.g. the determination of the W mass, as well as in astro- 
and astroparticle physics, e.g. improving evolution functions, have lead to the 
improvement of predicted GZK neutrino flux. If UHECRs are protons origi- 
nating from AGN, an detection rate of about 1 — 2 events can be expected in 
IceCube per year. This event rate is relatively low, since the GZK neutrino 
flux becomes important at energies above -Emin = 10^ GeV, while the sensi- 
tivity of IceCube peaks at lower energies around Ei, ~ PeV. Still, the flux 
should become visible within only a few years of operation due to the extremely 
low background at these high energies. The event rate for the GZK neutrino 
flux is expected to be significantly higher with future detection methods, like 
the acoustic and radio detection. Here, much larger arrays of the order of 
> 100 km^ are planned to be instrumented. For acoustic detection, the lower 
energy threshold lies at around ~ 10^'^ GeV, for radio detection, it is close to 
~ 10^'^ eV. Apart from underground arrays, neutrino-induced electromagnetic 
and hadronic showers can be detected with balloon experiments, e.g. ANITA, 
airshower experiments like Auger and Cherenkov telescopes like MAGIC. All 
these methods have relatively high detection thresholds. Therefore, they are 
not suited for the observation of point source fluxes from e.g. AGN, but are 
well-suited for the detection of the diffuse GZK neutrino flux. 

Within the next years, it will be possible to further improve the model of the 
GZK neutrino flux either by the detection of a signal or by setting limits. 
The sensitivity of a 45 days-flight by ANITA is indicated in Fig. 28. Here, it 
is already possible to constrict strong evolution models, like those assuming 
that UHECRs come from GRBs. Together with the detailed examination of 
the GZK cutoff observed with charged Cosmic Rays by the Auger observatory 
and other airshower experiments, the search for the GZK neutrino flux will 
render possible the constraint or the determination of the evolution function 
of the sources of UHECRs. 
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10 Summary and conclusions 



The concept of multimessenger physics becomes more and more important 
these days in order to get a complete picture of the sources of investigation. 
In particular, the non-thermal Universe can only be explained by combining 
pieces of information from all three messengers, i.e. photons, charged Cos- 
mic Rays and neutrinos. The connection between the observation of charged 
Cosmic Rays and high-energy photons can only be established by means of 
high-energy neutrino astrophysics. Hence, the deployment of 1 km^ neutrino 
detectors is crucial in order to further increase the sensitivity to the flux from 
extraterrestrial neutrinos. 

The modeling of neutrino emission from different sources is essential for the se- 
lection of data in large volume neutrino detectors. The different source classes 
reviewed here, i.e. galactic sources, AGN, GRBs, starbursts and GZK neutri- 
nos, are all analyzed with different methods. These include the search for a 
diffuse signal at the highest energies, single point source searches, transient 
source searches and stacking approaches. The different methods help to reduce 
the background towards the signal expectation. 

The fact that neutrino flux limits can be used for physics constraints can be 
formulated as the Olbers paradox for neutrinos. Olbers pointed out in the 
19th century that the fact of a dark night sky contradicted the prevailing 
theory of the Universe: a homogeneous distribution of stars in an spatially 
and temporally infinite Universe in thermal equilibrium implied an infinitely 
bright sky during day and night time. The simple limit of the dark night sky 
was a first indication that the world view at the time was incomplete. Today, 
it is known that the Universe is not in thermal equilibrium, which, among 
other things, leads to the dark night sky. In analogy to this optical limit, the 
so far dark neutrino sky provides first indications of how to improve current 
neutrino flux models. The current neutrino flux limits can already be used 
to restrict different models. In particular, a correlation between X-rays from 
AGN and neutrinos can be excluded, since models based on this assumption 
violate recent neutrino flux limits. 

Next generation of neutrino telescopes like IceCube and KM3NeT enables a 
deeper view into the neutrino sky by more than an order of magnitude. With 
new detection techniques, like the observation of neutrino-induced cascades 
at radio wavelengths, or the detection of an acoustic cascade in ice or water, 
enables the investigation of neutrino emission at the highest energies. The 
combination of different observational techniques renders possible the detailed 
investigation of a neutrino signal and thereby further advancement in the 
understanding of multimessenger astrophysics. 
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